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PREFACE. 



The present volume contains all that Professor Boole 
wrote for the purpose of enlarging his Treatise on Differential 
Equations. Had he lived to publish the second edition he 
would doubtless have incorporated his more recent investi- 
gations with the original work, and it is therefore necessary 
to explain why another plan has been adopted. 

In some cases Professor Boole had indicated that certain 
portions of the original work were to be omitted and their 
places supplied from the manuscripts; but on examination 
it appeared that in subsequent passages of the work there 
were references and allusions to the portions thus marked to 
be omitted which would not apply to the substituted matter. 
Thus in attempting to carry out the directions it would 
have been necessary to accept the responsibility of making 
many alterations, and consequently to incur the risk of fail- 
ing in the attempt to improve the original form. 

Moreover the Treatise had been for some time out of 
print, and the long delay which must have been caused by 
the labour of reconstruction would have produced serious 
inconvenience to students at Cambridge and elsewhere. Pro- 
fessor Boole himself was always especially anxious to consult 



the advantage of students, and those who had the cliarge 
of hia manuscripts were naturally inclined to adopt a course 
of which they believed he would himself have approved. 

The design of reconstructing the Treatise was therefore 
ahandoned; and it was resolved that the original volume 
should he reprinted, and that the manuscripts should be 
collected and published separately. This plan has the ob- 
vious recommendation of enabling those who are already 
familiar vrith lie original work to turn their attention 
readily to the new investigations, It will be seen that 
many of the Chapters of the present volume may be re- 
garded as independent essays or memoirs which lose nothing 
by being separated from the other volume ; and indeed no 
indications had been left by Professor Boole of the place 
which such Chapters were to occupy in the enlarged edition. 

I have printed all the unpublished matter relating to 
Differential Equations which I fomid among Professor Boole's 
pa|>ers. In a few cases it will be seen that an investigation 
is incomplete ; such investigations have however been in- 
cluded in the volume, because I was unwilling that anything 
should be lost which so gi'eat a mathematician had written 
on a subject he had long and carefully studied. 

I trust that no serious error will be found in the volume, 
and that any faults which may be detected will be excused 
on account of the nature and diiSculty of the task that had 
to be performed. Many of the manuscripts had not been 
finally revised ; some of them were very obscure and had 
to be carefully and laboriously copied for the press. In 
general the equations were not numbered, and thus only 
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blanks occurred in place of references; this drcumstance 
often caused great trouble and peiplexitj: I bope however 
that a satisfactoij result has been finally attained. 

1 maj state for the benefit of those who are conversant 
with the first edition of the original work that the theo- 
rem which in the present volume is cited as contained in 
Chap. n. Art. 1 will be found in Chap. IV. Art. 2 of 
the first edition: the change was made hj the direction of 
Professor Boole's interleaved copy. It was judged conve- 
nient to number the Chapters in the present volume in con- 
tinuation of those in the original work. 

All additions of mj own are enclosed within square 
brackets. The sheets have been read by the Eev. J. Sephton, 
Fellow of St John's College, as well as by myself, and the 
volume is much indebted to his care and accuracy. Obvious 
mistakes in the manuscripts were of course corrected; thus, 
for example, the table at the end of the volume was calcu- 
lated by Mr Sephton, because the table in the manuscript 
was rendered erroneous by the use of a wrong sign in a 
formula. 
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CHAPTER XIX. 



ADDITIONS TO CHAPTER II. 



1. [In Chapter ii. Art. 9, two methods are given for 
solving the differential equation 

{ax + Jy + c) tZa? + {ax + Vy + c') dy = 0.] 

But there exists another transformation by which the equa- 
tion may be reduced to, (because it may be constructed from), 
an equation in which the variables are separated. 

Assume as this equation 

(^y+C)(?aj'+(^V+(7Vy' = (I) 

and let a?' = a? + w^y, y' = a? + m^. 

It will be seen that in these equations united we have as 
many constants as in the original equation. Now on substi- 
tuting in the assumed equation the values of cc' and y', and 
comparing with the equation given, we deduce a system of 
relations equivalent to the following, viz.: 

The quantities tw^, m^ are roots of the quadratic 

am' — (6 + a') w + J' = 0. 

The c^uantities -4, A^ (7, C are determined by the system 
of equations 

A^A^a, 0+C'=^c, 

Am^ + Am^ = a , Cin^ + C'm^ = c, 
B.D.E. II. \ 

2 c 
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from which we find 

. am^ — a' ^ __ cm^ — c 

A. = — ) O • 

Now (1) gives on dividing by {Aoi + C) {^Ay + (7) and 
integratin 



g 



^ log (^V +0^0+2 ^^S (^2^' + ^) = ^^^^^-^ 

or {A^ + C) ^' (^y' + CJ^ = const, 

which on substitution and reduction gives 



{(am, — a') (a? + m, v) + cm. — c>"»i^' ^ . , 

■i^^ — ^ — ^^ ^ ^— i— = const (2) 

{{am^ — a) {x + m^y) + cm^ — c'}«^V^' 

2. Under certain circumstances the general solutions of 
differential equations of the first order fail. This happens in 
the above example if m, = m^, the solution then reducing to 

X.= const. 

The theory of the deduction of the true limiting form of 
the solution in such cases requires a distinct statement. 

Let the supposed general solution be represented by 

G being the arbitrary constant and u a function of a?, y, and 
constants which are not arbitrary. Suppose too that when 
one of these constants k assumes a particular value /c, the 
function u reduces to a constant v. Then we have 

w — V G — v 

fC "~ tc fc ~~ /c 

Now the second member being a function of an arbitrary 
constant is equivalent to an arbitrary constant and may be 
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replaced by C. The first member is a vanishing fraction, the 

limiting value of which is i-^A, the brackets being used to 

denote that after the differentiation k is to be made equal to k. 
Hence the solution becomes 



(: 



S)-« 



In applying this theory to the reduction of the general 
solution (2) in the case in which m^ = m^^ it must be 
observed that the numerator of the first member is the same 
function of m , a?, y, as the denominator is oi m^, x, y ; or 
attending solely to their functional character with respect to 
m^, 7n,, we may affirm that the numerator is the same function 
of m^ as the denominator is of m^. Representing these func- 
tions by <^(mj, <l>{m^ respectively, we have 

But wij, w, being roots of a quadratic equation may be 
represented in the form 

m^^m + k, m^^m — ky 

the roots becoming equal when k = 0. Hence 

__ ^(m + k) 
(j>{m — k)* 

Therefore since 

d<f) (m + k) _^d(p(m + k) d^ (m — k) __ d^ {m — k) 

dk ~" dm ' dk dm ' 

we have 

^^ <}>{m-k) ^^^^ ^ -.<l>{m^k) ^^ 

dk ld)(m-A:)f ' 
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therefore ^^j |^(«,)j^ = -J(^ 

c 

Thus the solution becomes on putting (7 for — , 

J 1 

3. [The next Article seems to have been intended to ap- 
pear in the enlarged form of Chap, ii.; but 1 cannot discover 
what precise position it would have occupied. I conjecture 
that " the above demonstration" refers to Chap. il. Arts. 2, 3; 
and I have accordingly supplied a reference to equation (3) of 
Chap. II. 

I had myself drawn Professor Boole's attention to Chap. ii. 
Arts. 2, 3. The geometrical process of Chap, ii. Art. 3, ap- 
])ears to have been first given by D'Alembert in his Opus- 
cules, Vol. IV. p. 255. D'Alembert calls it a demonstration; it 
seems to me only an illustration, at least in the brief form of 
the text : and that such was Cauchy's opinion may perhaps 
be inferred from the elaborate investigation given by Moigno, 
to which Professor Boole refers in Art. 5 of the present 
Chapter. 

I had also drawn Professor Boole's attention to the state- 
ment at the end of Chap. il. Art. 12, that only one arbitrary 
constant was involved. Accordingly Article 5 of the present 
Chapter developes this statement, and Article 4 seems intended 
to bear on the same subject.] 

4. In the above demonstration the relation between y and 
X is regarded as one of pure magnitude, and the interpreta- 
tion of the differential equation becomes a limiting case of 
that of the equation oi fiiitte differences (Eq. (3), Chap. ii.). 
But if we represent x and y by the rectangular co-ordinates 
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of a moving point on a plane the differential equation may be 
interpreted directly. For supposing it reduced to the fona 



|-=/(^,y). 



we see that the direction of motion is constantly assigned as 
a function of the co-ordinates of position. The entire motion 
is therefore determinate as soon as the initial point is fixed. 
The result of the motion is a line or curve wholly continuous 
or subject to in'egularities according to the nature of tlie func- 
tion f(x^ y). That the arbitrariness of origin is geometri- 
cally equivalent to the appearance of a single arbitrary con- 
stant in the relation connecting x and y may be shewn thus. 

Let y = ^(^o»yo»^) 

be the relation between x and y indicated by the supposed 
motion, a;^, y^ being the initial point of departure. Then this 
point being on the line of motion, a;^, y^ are particular values 
of X and y, so that we have from the above equation 

which establishes a relation between x^ and y^, and shews 
that there exists virtually but one arbitrary constant. 

5. It is proved in Art. 3, Chap, ii,, that the constants 
a?^, y^, initial values of the variables a?, y in the solution of 
the differential equation of the first order, are necessarily 
equivalent to one arbitrary constant. I shall shew from the 
form of the above solution that this a priori condition is 
actually satisfied. 

Developing the expression for y [see Eq. (30) of Chap, ii.] 
in ascending powers of a?, we have 

y = A + Aa' + ^+l^ + &c., (32) 

in which A, = ^f^^^^JiAtL<!ll 

1,2 ... (» — r) 

the summation extending from n — r to n = oo. Forming 
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hence the differential coefficients of A^ with respect to x^ and 
y^^ and reducing by (28), we shall find 

dA^ - , . dA^ ^ 





dx^ 


-fJiK^t, 


''' dy. 


-V/, 


whence 


in particular 












+/i(«.. 


^''' dy. 


= 0. 



Eliminate between these equations^ (a?^, y^^ and we have 

dA^ dA^ ^ dA^ dA^ __ ^ 
d^o d.y^" dy^ dx^" 

Therefore, by Prop, i., A^ is a function of -4^^, so that tlie 
solution reduced to the form (32) contains but the single 
arbitrary constant A^. 

It remains to notice that the solution must be applied 
only under the conditions of convergency, i.e. under the con- 
dition that the ratio of the w*^ to the (n — 1)*^ term tends to a 
limit less than unity as n tends to infinity. For a discus- 
sion of the failing cases of this test see ' Finite Differences,' 
Chap. V. Generally it is desirable, in order to secure rapid 
convergency, to divide the interval x — x^ into separate equal 
portions, to each of which the general theorem of solution 
may be applied. If a? — a?^ be very small the theorem may 
be approximately represented by 

2'-yo=/(^o,yo)(-^~^o)- 

On these principles Cauchy has founded remarkable methods 
of solution, which deserve attention from the commentary on 
the limits of error on their application by which they are 
accompanied (Moigno, Vol. li. pp* 385—434). 
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CHAPTER XX. 



ADDITIONS TO CHAPTER VII. 



1. [This Article relates to Art. 2 of Chap, vii.] 

The sense In which (9) may be said to constitute tlie 
general solution of the differential equation is this. We 
obtain from it 

giving any particular value to G this will geometrically 
represent a curve consisting of two branches, and giving to 
C every possible value we obtain an infinite system of such 
curves, each consisting of two branches. The aggregate of 
branches thus obtained is evidently the same as tlie aggre- 
gate of curves given by the two primitives (5) and (6), un- 
restricted by any connexion between c^ and c^. In this sense 
then the solution (9) is general, that it includes all the parti- 
cular relations between y and x which are deducible from 
the original primitives (5) and (6). And it is only in this 
sense not general that it groups these relations together in a 
particular manner. 

To the expression of the complete primitive a certain 
variety of form may be given without affecting its generality 
in the sense above affirmed. Thus, if to the solutions of the 
component differential equations we give the forms 

ye^-c^ = 0, logy + aa? - Cj= 0, 

we should have, by the same procedure, as the expression of 
the complete primitive, 

{]/€^ - c) (logy + oa; - c) = 0, 
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an equation which may equally with (9) be regarded as the 
complete primitive of the differential equation given, and 
which in geometry represents the same totality of branches of 
curves as (9), with this difference only, that they are differ- 
ently paired together. 

2. [This Article relates to Art. 3 of Chap, vii.] 

The question will here naturally arise, Since if F= c be 

a solution of one of the component differential equations, 

/{V) = c, in which /(F) denotes any function of F, is also a 

solution, by Chap. iv. Art. 3, why not give to the complete 

primitive the form 

{/.(^J-cH/,(F,)-c} {/.(F.)-c]=0, 

or the stricter form 

A{y^)fAy.) /«(F.)=o , {F), 

in which /(F,), f^{V^^ •••/n(l^n) denote arbitrary functions 
of Fj, Fg,..., V^ respectively — stricter because the presence of 
arbitrary constants and functions in the previous form is a 
superfluous generality ? It is replied that though the form 
just given is analytically more general than (15), it is not 
more general than (15) with such freedom as is permitted 
in the interpretation of the arbitrary constants. In a physi- 
cal or geometrical application we should not only be per- 
mitted to assign a particular value to the arbitrary constant 
in (15), so deducing what in reference to its source would 
then be termed a particular primitive, but to combine the re- 
sults of different determinations of c together, so as to obtain 
every form of solution which is implied either in the func- 
tional equation (-F), or in its component primitives 

The same considerations justify us in speaking of (15) as 
the complete primitive, and not as a complete primitive. 
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CHAPTER XXI. 



ADDITIONS TO CHAPTER VIII. 



1. [The Singular Solutions of Differential Equations of 
the First Order received great attention from Professor Boole, 
and the Chapter devoted to that subject is one of the most 
valuable and important in his work. He continued his re- 
searches after the publication of his first edition, and intended 
to reconstruct the Chapter with great improvements in the 
second edition. After carefully examining the manuscripts I 
came to the conclusion that it would be very difficult to re- 
write this portion of the work so as to connect the old matter 
with the new ; and thus it seemed best to reprint the original 
Chapter Vlll. with corrections of obvious misprints, and to 
print the matter intended for the revised form in the present 
volume. The plan gives rise to some repetition ; but this 
seems unimportant, compared with the advantage of preserv- 
ing in the author's own language all that he left on an in- 
teresting and important point which he had carefully studied. 

2. It may be of service to the student to reproduce the 
substance of some remarks on his Chapter viii. which were 
sent to Professor Boole soon after the publication of his first 
edition ; for there is evidence in his manuscripts that he paid 
great attention to sucli remarks while engaged in the revision 
of his work, and thus the reason and the meaning of some of 
his additions and changes may be made more obvious. These 
remarks will occupy the next Article. 

3. The two pages beginning with " And these conditions 
are sufficient," and ending with " do not lead to conflicting 
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results,''^ forming part of Arts. 3 and 4 of Cliapter viii., seem 
obscure and difficult. The following may perhaps be substi- 
tuted with advantage. 

The only ways in which 

dy ^ df{x, c) ^^^ dy ^ df{x, c) df{x, c) dc 
dx dx dx dx dc dx 

can be equivalent when c is variable, are 

(1) when^ = 0, 

(2) when^^)=oo; 

in the latter case -^^ = go , and therefore -r- = 0, and this 

ax dy 

implies that the singular solution is of the form x = constant. 

Thus there can be no singular solutions except such as 

df (x c) 
are found from "^^ - — - = 0, and such as are found from 

dc 

X = constant. 

Similarly, if the complete primitive be expressed in the 
form x — F{y^ c), there can be no singular solutions except 

such as are found from -j^ — = 0, and such as are found 

dc 

from y = constant. 

In Art. 8 of Chapter viii. we read, " We may pass over 
the case in which the above equation is satisfied independ- 
ently of c, because the relation obtained would involve x 

only, while it is a condition accompanying the use of -r- = 

that it leads to solutions involving y at least." It is ob- 
jected, Why may we pass over this case? Such a case might 
occur and furnish a solution, and then we should want to 
know the character of that solution. Take for example 

jp = cc*y ; here if n is negative, -~ is infinite when a; = 0, and 
this is a singular solution. For the general solution isy =C€»+i, 



QO 
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and 80 a?=0 is not a case of it. The words — lohih it is 

a condition,,. at least — seem very diflBcult, for by supposition 

dt) 
we are now investigating what is furnished hy ~ =cc . 

Professor Boole met the objection in substance thus : 

" It will be found that the rules in the book are correct in 
this case. What is implied in the Chapter, though not stated 

with sufficient clearness, is that if ^ = co leads to a solution 

ay 

which does not involve y in its expression, nothing is to be 
inferred whether it is singular or not. Then the proper test is 

dx \p) 

• In tliis example we have 

-^ = Qo gives cc" = CO ; no inference ; 



dx \pj 



GO gives a?"^*"^*^^*^ = co . 



Hence a; = 0, provided n is between and — 1, or y = 0. 

Consider these separately : 

First. Let w be between and— 1, and x = 0. This is 
by the test a singular solution. Substituting it in the com- 
plete primitive we get y = c, which confirms this. 

Second. Let y = 0. This satisfies the differential equa- 

J /I \ 
tion; but from the fact that it comes from ^ f-j = Go we 

have no inference ; from the fact that it does not come from 

-^ = 00 we have the inference that it is a particular integral: it 

corresponds to c = 0. 
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There remains the case of a: = when n is between — 1 

and — QO . As this does not satisfy -7- ( - j = 00 , we infer that 

it is a particular integral. To prove this we have 

^+1 



When a; = this gives, since 1 + w is negative, 

c = co orc = — CO, 

according as y is positive or negative. This is like Ex. 2 of 
Chap. VIII. Art. 8." 

The remark made hy Professor Boole in the above reply, 

that «y -p = CO leads to a solution which does not involve y 

nothing is to he inferred,,. is important. It corrects the state- 
ment put too strongly in Chap. viii. Art. 7, " All we can aflSrm 

is that if -^ = CO gives a solution at all it will be a singular 

solution." 

From Art. 8 onwards it seems assumed that a solution for 

which -^ = is always to count as a singular solution, even if 

it should coincide with a particular integral. This does not 
seem to have been quite tne view of the former part of Chap- 
ter VIIT. : see Arts. 5 and 6 of the Chapter. 

In Ex. 3 of Art. 9 we read, " the second is obviously a 
singular solution." This means that since we have a solu- 
tion which makes -^ infinite, we conclude that it is a singular 
solution. 

So in Ex. 5 of Art. 11 we read, " is evidently a singular 
solution," when it seems better to say, ** and is therefore a 
singular solution," 

4. The additional matter relating to Chapter vili. begins 
with another example which was to be placed at the close of 
Art. 3 of that Chapter.] 
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Ex. The differential equation 

has for its complete primitive 

Vaj* + y — 7/1* — y — c = 0. 
Here ^^ , ^ -1, ^ ^ 



^ = -1. 
Hence -—- = ' ^ =, -r- = =^^ . 

« 

Both -T^ and -j- vanish then if 
dc dc 

This therefore is the singular solution and it satisfies both 
the tests, as both x and y are contained in its expression. . 

Of the partial tests 

c?c ' da? ^ dy * 
the first is not satisfied, the last two are satisfied. 

The determination of c as a function of x by the solution 

dfix c) . . 

of the equation -^ ^' ^ = is equivalent to determining 

what particular primitive has contact with the envelope at 
that point of the latter which corresponds to a given value 
of x. 

One important remark yet remains. The elimination of c 
between a primitive y =/ (a?, c) and the derived equation 

-T^ = 0, does not necessarily lead to a singular solution in the 
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sense above explained. For it is possible that the derived 
equation 

dc 

may neither on the one }iand enable us to determlue c as a 
function of at, so leading to a singular solution ; nor, on the 
other hand, as an absolute constant, so leading to a particular 
primitive. Thus the particular primitive 



being given, the condition ->- gives 



-0, 



^lm^ 



vrhence c is + oo if ic be negative, and -co if cb be posil 
It is a dependent constant. The resulting solution y = 
does not then represent an envelope of the curves of particu- 
lar primitives, nor strictly one of those curves. It represents 
a curve formed of branches from two of them. It ia most 
fitly characterized aa a particular primitive marked by a sin- 
gularity in the mode ol its derivation from the complete pri- 
mitive. 



All the foregoing observations and conclusions may be 
e.-itended to the case of solutions derived from the condition 

5. We Iiave seen that the equation zr =0 may be satisfied 

by an absolutely constant value of c, so leading to a particu- 
lar primitive and not a singular solution. In this case 
■^ {a:+ h, c) as well as i^ (x, c) would vanish, and the nume- 
rator of (9), instead of being.the difference of a finite and an 
infinite quantity, would be the difference of two infinite and 
equal quantities. [See Chap. Viii, Art. 8.] It would not there- 
fore be infinite. Hence we conclude that -,- would not become 
infinite for a particular primitive In the strict sense of that 
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term, i. e. for a solution derived from the complete primitive 
by giving to c an absolutely constant value. 

This is one point of contrast between the conditions 

^^ = ^ = 00. 
dc ^ dy 

There is another not less important. As the numerator 
of (9) may become infinite not only when -^ [x, c) = 0, but 
also when yfr (a?, c) = infinite, we see that a relation between 

y and x which makes ^ infinite will not necessarily satisfy 

the differential equation. On the other hand, it is not a par- 
ticular primitive in the strict sense of that term. 

. . dx 
Exactly in the same way the condition -i- = 0, as relating 

to the complete primitive, leads to the condition 

as relating to the differential equation, with the same points of 
difference in the respective applications. 



dy__ 



m-l 



Ex. Let -p- = my ^ , and suppose m a positive constant 
greater than 1. 

Here ^ = (jti - 1) y^, 

which becomes infinite when y = 0. As this involves y and 
satisfies the differential equation it is a singular solution. 

To confirm this conclusion we may refer to the complete 
primitive 

which does not give y = for any particular value of c. 

Now let m be a positive constant less than 1. We have 
still -T- — ^ when y = ; but this value of y no longer satis- 
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fies the differential equation. It is not a solution at all, nor 

would it result from the application of the condition ^ = 

to the complete primitive. The distinction of character of 
the two tests is here made manifest. 

6. We may express the most important results of the 
foregoing investigations in the following theorem. 

Theorem. Every solution of a differential equation of the 
first order which is derived from the complete primitive by 
giving to c a variable value will, if it involve y in its expres- 
sion, satisfy the condition 

dp 
and if it involve a?, the relation 



dx \p) 



But relations satisfying these conditions will not neces- 
sarily be solutions of the differential equation. 

In applying this theorem the following points must be 
carefully attended to. 

1st. No conclusion can be drawn from the satisfying of 
the condition -^ = oo when the relation in question does not 
contain y in its expression, nor from the satisfying of 

— f-l = oo 
dx \pj 

when the relation in question does not involve x in its ex- 
pression. For these conditions being respectively derived 

d'U dx 

from -f-=0 and ;t— =0 are subject to the same limitations 

in their application, 

2ndly. It may be that ^ ot -r- ( -J assumes for a particu- 
lar relation between x and y the indefinite form - • In this 



HIbt. 6,] ADDIi'lOSS TO CHAPTER VIII. 17 

case we must seek by the development of its terras or by 
other known modes its true limiting value or values. Finite 
values will indicate particular primitives, infinite values sin- 
gular solutions, and when such values emerge together out of 
the same relation between the variables, the solution wiU ha 
a particular primitive poaaeasing the geometrical propertiea of 
a singular solution. Ifa locus will be a particular eurve en- 
■Veloping other curves of the same family. 



See Examples 2 and 3 of Chap. viii. Art. 11. 
We have seen that the conditions 



dy 



■Velo 

I 

^Bidicate in general the existence of a relation between c and 
^Hb or c and y. And when that relation is such as to eoable ua 
™lto determine c as a continuoua function of one of the vari- 
ables, the corresponding aolution of the differential equation 
is singular, and is geometrically represented by an envelope 
of the curves of primitives. But it may be, as we have seen 
in a particular example, that the relation does not determine c 
aa a function of a: or ^ ; but according to the language already 
used, c is a dependent constant, or in some other way different 
_&om the constant of an ordinary particular primitive. Let 
|a examine in particular instances the kind of singularity 
vhich may hence arise. 



Ex.1. Giv 
Here 






This becomes infinite if a; = 0; but this not involving y 
must be rejected. Again, it becomes infinite if y = 0, and 
biiis proves to be a solution of the differential equation, the 

jlimitrng value of the indeterminate function in the second 

■member l>eing (Todhunter's Differential Calculus, Chap, x.). 

■ "ow the complete primitive is y = e°, diacusscd in Art. 4. 

Il'he constant c is there shewn to be dependent, the solu- 
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tion y = emerging from the complete primitive by making 
c = — 00 if a? be positive, and c = oo if a? be negative. 

Ex. 2. Given (^ ^ a^ ^ + y' log y = 0. 
Here ^^xy ±y{a?^^\o^y)^ ^^ 



therefore 



2 



dp __x ±(x^ — 4: logy) ^ _ 1 

dy 2 (a?'-4logy)i' 



and this is made infinite by y = and by cc* — 4logy=0, 
i.e. by 

y = o> y=€*. 

Both satisfy the difierential equation. 
Now the complete primitive is 

We see at once therefore that the second of the above solu- 
tions is singular. The first however is deducible from the 
complete primitive by making c = co or c = — oo, irrespec- 
tively of the sign or value of a?, provided only that x he 
finite ; not so however if x be infinite. The value of c is not 
therefore in the most absolute sense independent of that of x. 
If from the complete primitive we seek the singular solution 

by the condition --^ = 0, we get the two equations 

€«^=0, aj-2c=0. 

The second of these determines c as a function of a?, and 
leads to the second of the solutions obtained above. The first, 
though it does not determine c as a function of x, still ex- 
presses a relation between c and a?, which is the ground of 
the fulfilment of the condition 

dp 
dy 
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► 



B>RT 

■ We may further notice a peculiarity arising from this rela- 
W Son. Supposing x finite and the solution j = a particular 
H. integral, it presents the singularity that it is the only case in 
I which two particular integrals agree. We might in any com- ■ 
plete primitive, by changing c into c", get two values of c for 
the same particular integral, bnt then it would be for every 
particular integral. 

One negative character seems indeed to mart all the cases 
1 which a solution involving y in its expression satisfies the 

condition hC=ix. It is that such solutions do not emerge 
from the complete primitive by the attributing of a single and 
absolutely constant value to c. The relation which makes -f 
infinite satisfies the differential ecjuation only because it satis- 
fies the condition -f^ = 0, and this implies a connexion be- 

ac '^ 

tween c and x, which is the ground of a real though it may 
be unimportant singularity in the solution itself. 

At this point, then, the question arises, whether the term 
singular solution shall be confined to that class of solutions, 
the loci of which represent the envelopes of curves of primi- 
tives, or shall be extended to all solutions which, satisfying the 

condition -7- = oo , indicate the existence of a relation be- 

tween c and x, and possess an actual singularity arising from 
this source. While the all but universal consent of mathe- 
maticians is in favour of the former course, it is to be remem- 
bered that the question is solely one of definition. Not such 
is the question how singular solutions of the envelope species, 
or as would more generally be said true singular solutions, 
are to be distinguished from all other solutions. This we 
now propose to consider. The question is not an isolated one. 
It stands in close relation to a series of properties of singular 
Bolutions which admit of an orderly development. 



1 
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DiscrimtTtation of singular solutions of the envelope ^ecies. 



7. A negative test, which in the ^eat majority of casea 
suffices for the present object, ia suggested by the followiDg 
consideration . 

The differential equation determining -^ aa a function of 

X and y determines also y-^ , -t4 , ad inf., and the know 

iedge of these enablea us to construct in a developed form 
the complete primitive. See Chap. ii. Art. 12. 

The values of —■ , -r^, &c. ad inf., as derived from the 

differential equation, are the same as tlioae derived from the 
complete primitive. 

= CO is only 

constructed so as to yield the same value of ,- as the given 

differential equation does. If it be of the envelope species, 
llio curve it represents has in general no continuous contact 
with tiie curve of any particular primitive. It will not there- 
fore generally yield the same values for j-\, -j^, &c. as 

tlie differential equation docs. It will not therefore generally 
aatisty the differential equations of an order higher than the 
tirst, which would be derived from the given equation by dif- 
ferentiatioQ. Hence we have the following Proposition. 

Proposition. If a relation which makes -y- infinite satufy 

the given differential equation of the first order, but do not 
satisfy all the higher differential equations obtained from ii^ 
such solution will be singular and oj the envel<^e species. 
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Ex. 1. By comparison with its complete primitive we saw 
in Art. 5 that -^ = m^ "• has for a singular solution y = 
when nt is a constant greater than 1. 

We will first suppose m a fractional quantity greater than 
1, and en(leavour to deduce the character of the solution with- 
out making use of the complete primitive. 

From the solution we have 

^ = 0, ^ = 0, &c, adinf. 

But from the differential equation 



g = rm^lU--^ = 






and generally 



m—r 



g=»(»-x,..,»-,.,,^. 

Hence, when r is less than m, the substitution of y = gives 
as before. But if r is greater than m, it gires 

We conclude that the solution is of the envelope species. 

Secondly, suppose m a positive integer greater than 1. 

In this case we find, when r is less than m, the same series 
of values as before ; but for r = m we have 

^ = m(w-l)...l, 

and this also shews the solution to be of the envelope species. 
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Ex. 2. The differential equation 

1+y 4 

is satisfied hj 

Is this a singular solution or a particular integral ? 

From the solution we find 

dy ^ X d^y __ 4 
dx^ y ' doi^ y^ * 

From the differential equation we shall have 

^^_(l±£!)V±i(!I 
doi? 2(y-a?p) * 

substituting in which the value of -J^ , obtained from the 
proposed solution, we find 

da? 2f /• 

Now this differing from the value before obtained, we con- 
clude that the solution is singular and of the envelope species. 

And this result is verified by comparing the solution with 
the complete primitive 

As the test above exemplified is merely negative, it is in- 
sufficient. For it is conceivable that an enveloping curve 
should have an infinite order of contact with each oi the curves 
which it envelopes, and this is also possible. Anv test found- 
ed upon a comparison of the values of differential coefficients, 
any test therefore famished by the Differential Calculus, would 
be insufficient for the discrimination of such cases. 



Ex. 3. Given -j- =?= y (log i/)*. 



dx 



\ 



t 
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Here 
r tial equation. 

From this solution we find 

From the differential equation we have 

_ Wliich consists of y multiplied by a rational and entire func- 
tion of log^. It is easy to see that all the higher differential 
coefficients of y lietice derived will possess the same character. 
And all 3uch vanish with y. 

We can therefore neither affirm nor deny that the proposed 
solution is of the envelope species, 

8. Before demonstrating a general Rule for the discrimi- 
nation of solutions of this character, we shall notice certain of 
their properties which serve to indicate in what direction the 
Rule IS to be sought [See Chap. viii. Art. 14.] 

Ab the exact differential equation differs from the sup- 
posed given differential equation by having acquired a factor 
which the singular solution makes infinite, so the given dif- 
ferential equation may fee said to differ from the correspond- 
ing exact one by containing a factor which the singular solu- 
■tion makes to vanish. If we knew that factor, we could "by 
rejecting it reduce the given differential equation to a form in 
wnich it would no longer be satisfied by the singular solution. 
Now PoisBon has shewn on a particular assumption, which 
does not however affect the principle of the demonstration, 

ItiiaX this factor can be found when the singular solution is 
Ildowu. His demonstration is in substance as follows. 
^ Let US represent l!ie given singular solution of the dif- 
ferential equation by 

w = 0, 
» }Klng a given function of x and y. Then introducing u and 



^ 



J 
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X instead of y and x as variables, the differential equation 
after transformation will assume the form 

g =/(«', «). 

Now this equation being satisfied by w = and the first 
member vanishing, the second must also, Poisson now 
assumes, and the assumption must be carefully noted, this 
second member to be capable of being developed in ascending 
positive powers of u. Supposing it so developed, the diffe- 
rential equation becomes 

in which A, B,... are functions of a?, and a, /8,... ascending 
positive indices. 

Hence if w = be a singular solution we have, putting j> 

^^ = Jaw-i + ^)9tt^i+ &c. = Q0 . 

But this demands that there should be at least one nega- 
tive power of ti in the development in the second member. 
Therefore a — 1, the lowest index, must be negative. There- 
fore a being already positive must lie between and 1. 

We may give therefore to the transformed differential 
equation the form 

a being a positive fraction, and Q not vanishing with w. 
Hence, dividing by w", 
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a differential equation wliicli is not satisfied by m = 0, since 
w = gives M^"" — 0, and the first member vanishes while the 
second member does not vanish. In its present form then 
the equation is not satisfied by ii = 0. We see also that the 
property of being satisfied by a = has been lost not in 
reality through a transformation, but through the rejection of 
an algebraic factor u" from the transformed equation. It has 
been shewn in the treatment of Clairaut's equation, how in 
the ascent by differentiation to an equation of^ a higher order 
a somewhat analogous effect is produced, the singular solu- 
tion emerging out of a factor of that higher equation. 

If we inquire what ia essential in Pwsson's demonstration, 
we shall find it to consist in that the transformed equation ia 
cf the form 



I in which while Q neither vanishes nor becomes infinite when 
r« = 0, the functions 



UslbA 



f'du. 



Klx)th TBniflh with w. The question whether U ia of the form 
■ tC as Poisson supposes, or ia not, ia wholly immaterial, 
ft'Thie will fiiily appear from the demonstration of the foUow- 
ving theorem, which ia in effect Poisson's freed from arbitrary 
lesumptions. 

Proposition, If u={i he a solution of a differential 
tguation ofthejirst order between y and x, and 
du ^ 
dx 



-/(=',«) 



weaenl the form which that equation assumes when u and x 

■ara assumed as variables instead of y and x, then if f{xi ") ^^ 

resolved into two factors Q, U, of which Q neither vanishes 

nor becomes infinite when u = 0, while the functions U and 

m \ -fi l>ol^ vanish when m = 0, then the differential equation can 

' J reduced to a form in which it shall cease to he 



J 
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In the statement of this proposition x is supposed to be 
constant in the integration relative to u. 

The differential equation after the transformation which 
introduces u and x as variables becomes 

1= «^- 



Let 






so that V is in general a function of x and w, the form of 
which is known by integration when that of U is given. 
And again, transform the differential equation by making v 
and x the variables instead of u and a;. We have 

dv\ dv dv du 



(S- 



dx du dx* 



dv 
in which ->— is the differential coeflScient of v with respect to 

a?, on the above hypothesis as to the constitution of v as a 

function of x and w, while {-^j is the differential coefficient 

on the hypothesis that v is reduced to a function of x alone 
by the conversion of u into a function of x. 

c*> dv 1 du ^Tr 

the above equation becomes 

'dv\ dv 



\dx) dx 



Now if t* = give v = for all values of Xy it will there- 
fore give 

^dv^ 



(S-- 



and further, 

* du 



dx dx] a 



AET. 9.] ADDITI0S3 TO CHAPTER VIII. 

I^iiice we are permitted to make m = before effecting the J 
differentiation with respect to x. Hence the equation re* 1 
daces to 



And thia is not satisfied, since by hypothesis Q does not 
vanish with u. 

Hence if m = make I -,, = 0, the ti'ansformed differen- 
tial equation will no longer Iiave m = for a solution. 

CoE. Assuming Q — \, which does not violate the hypo- 
esia respecting Q, and gives 

e see that if 

! satisfied by m = 0, and if at the same time k = gives 



he differential equation can be transformed so as to cease 1 
admit of the solution u= 0, 

It is obvious however that it ia best to assume Q so as to 
make tlie subsequent integration for determining v the sim- 
plest possible. 

It is manifest that a solution which can thus be made to 
ease to satisfy the differential equation cannot be a particular 
nimitive. I'or the complete primitive of the transformed 
ifferential equation which it does not aatiafy is convertible 
)nto the complete primitive of the original differential equa- 
tion which it does satisfy, merely by writing therein for v its 
expression as a function of x and v. It cannot therefore be a 
ise of the complete primitive m any sense. It must be 
singular solution of the envelope species. 

The converse proposition still remains to be proved. 



28 ADDITIONS TO CHAPTER VIII. [CH. XXI. 

10. Proposition. If u^Ohe a singular solution of the 
envelope species of a differential equation of the first order ^ and 
if by assuming u and x as the variables^ the differential equa- 
tion is reduced to the form 



then will 



du 



Jo 



become when w = 0. 

Let the complete primitive be represented by 

F{x, u) = (7, 
then, since 

dF(xy u) dF{x, u) ^ _ /. 
dx du dx ' 

we have if for brevity we represent F{x, u) by F^ 

dF 

du dx ^ 

du 



« 



dF 



therefore / ^.^ , =- / -y^rfw- 

Jof{x,u) I dF 

^Q dx 

Now w = being a singular solution, F{x, 0) is not a con- 
stant ; for if it were, the complete primitive would, on giving 
to G the constant value in question, yield w = as a particu- 
lar primitive. And this would equally be the case whether 
that constant were finite or infinite in value. We see then 

that F{x, 0) must be a function of x, and therefore — \ ' 

must either be a function of a;, or a finite constant difiering 
from ; the latter if F{x, 0) be of the form ax-\-b, the former 
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if it be not of that form. Therefore the value of 7-^ — - 

ax 

Avhen w = 0, since in this we are permitted to make w = 

before differentiating with respect to a?, will be a function of 

Xy or a finite constant differing from 0. 

Now if is manifest that in general 

r* 1 dF . ^rdF 



J dx 




where H is some value intermediate between the greatest and 

least values which -jrv assumes within the limits of integra- 

dx 
tion. When these limits are, as in the above case, infinitesi- 
mal, we have 

dx 

TT I du 1 f^dF , 

Hence I -tt r= \n. ... I -^ du 



JoA^yU)'' dF(x,0) l du 

dx J. 



dx 

dx 

But we have seen that \ ' does not vanish. Hence 

ax 

its reciprocal, the first factor of the right-hand member of the 

above equation, does not become infinite. Again, 

F{x,u)-F{x,0) 
vanishing when tt = 0, we have 

'* du 



1. 



when u is made infinitesimal as was to be shewn* 
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It will be observed that the previous general expression 

for I 7- -. becomes infinite if M = ia a particular integral. 

For then, F{x, 0) being a constant, — ' ' Tanishea, wbilc 

^(a;, «} —F{x, 0) does not vanish bo long as w differs by 
however small a quantity from 0. 

These propositions form the ground of the following Rule 
for the discrimination of singular solutions of the envelope 
species from all others. 

11. Hole. The proposed solution being represented by 
M = 0, let the differential equation, transformed by makj 
and ic the variables, be 



dx 
Determine as a function of x and h the integral 






lakin^^l 



in which U is either equal to f{x, u), or to f{x, u) deprived J 
of any factor which neitlier vanishes nor becomes infinite 
when K = 0. If that integral tend to with u the solution la 
singular. 

Ex. 1. Determine whether y = is a singular aolutio^or 
particular integral of the differential equation 



io^j^J 



Here, since u = y, no preliminaiy transformation is needed. 



We have 



^y 



if (logy)" logy' 
which tends to with y. Hence the solution is singularJ 
To verify this we obaei-ve that the complete primitive 
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wd tills cannot be reduced to ^ = by giving any constant ' 
^alue to c. 

We have seen in Art, 7 that the test which ia founded upon 
the compariaon of differential coefficients does not suffice to | 
xsharacterize the above solution, 

Ex.2. The equation 4' = ^ - - is satisfied by if=0. 
Js this solution singular or particular? 

Here also no transformation is required. We have, reject- 
ing the factor — which neither vanishes nor becomes infinite 
|rlieny = 0, 

I — i-^^ — = loe los V ~ log loET 

id this being infinite, however small y may be, may properly I 
1 said to tend to iutinity as y tends to 0. The solution ia 1 

therefore particular. 

It will perhaps appear at first sight as if in the above es- 

iimiple we ought to write 



log 



iogy „ 



logo 

rhen y Is made equal to 0. But the course of the demonsfra- 
ion shews that the value of the definite integral must be first 
ibtained on the hypothesis that u (in this case replaced by y) 
B finite, and then the limiting value which its expression 
approaches to, as u approaches to 0, he sought. And in this 
case, since for all finite values of u however small the integral 
is infinite, its limiting value is infinite. 
The complete primitive in the above case is 

and the nature of the solution y = has already been dis- 
cussed in Art, 4, 
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12. It is remarkable that wliile the theory of enveloping 
curves and surfaces was at once founded and developed by 
Leibnitz in 1692 — 4*, the corresponding theory of the singular 
solutions of differential equations has been of very slow growth. 
1'he existence of these solutions was first recognised in 1715 
by Brook Taylor; it was scarcely more than recognised by 
Clairaat in 1734. Euler, in a apecia! memoir, entitled Expo- 
sition de quelques Paradoxes dans le Calcul Integral, publifiiied 
in the Memoirs of the Academy of Berlin for 175C, tirst made 
them a direct obje6t of investigation ; but the foundations of 
their true theory were only laid in 1768 in his Institutional 
Calculi Integralis. Laplace, Lagrange, Legendre, Foisson, 
Cauchy, and De Morgan have in various ways developed and 
extended that theory; but tliere has been so remarkable a 
want of unity and connexion in this long series of researches, 
that important portions of the theory appearing in a too 
isolated form have been neglected, forgotten, and rediscovered. 
I purpose here to give a brief account of what seems most cha- 
racteristic, rather than of what is most original in their several 
researches ; for the germs of nearly ail subseqnent discoveries 
on the subject are to be found in the great work of Euler. 

Taylor and Clairaut appear to have been led by accident to 
the noticing of singular solutions; the former while directly 
occupied on the solution of differential equations, the latter 
while discussing a remarkable class of problems relating to 
the connecting properties of different branches of the same 
curve. Taylor gave them the name singular, white Clairaut, 
and Euler too in bis memoir, regarded them as a species of 
paradox, not merely from their non-inclusion in the general 
integral, but from the mode of their discovery through a 
jfrocess of differentiation. The memoir of Euler, though it 
sheds no Light on the real nature of tliese solutions, contains 

* Ada Emditonm, 16S2, p. 168; 16S4, p. 311. Opera, Tom. in. pp. 284, 
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an interesting theorem concerning their connexion with the 
form of the differential equation, viz. If this equation can 
be brought to the form 

Vdz=^Z{Pdx'VQdy), 

in which 2 is a function of x and y, and Z of z, then will 

be a singular solution. In his Institutiones Calculi Integralis^ 
Tom. I. p. 393, however, Euler gives a rule which is the 
counterpart of that of Cauchy. [See Chap. Vlil. Art. 12.] 
He shews that if w = be a particular integral, and if the 
differential equation be reduced to the form 

f* du 
then -PT r = 60 . 

-'o <^(^> w) 

The limits of integration are here supplied. The reasoning, 
which is not fully developed, is the following. From the 
transformed equation we have 

, du 

ax = 



Hence a; = (7 + 1 



du 



(f> {x, u) ' 
X ^ 1 f du 



; ^ c] 



G ^ CJ<I> {x, u) • 

If this be satisfied by a solution involving x and y, and if 
that solution be a particular integral, then on putting for x 
its value in terms of u and integrating, the above equation 
will be satisfied by giving some particular constant value to 
C, But if the supposed particular integral be t^ = 0, then x 
and u being independent, we may perform the integration 
with respect to t^ as if a; were constant. The resulting equa- 
tion cannot be free from x unless C be infinite, and then it 

B. D.E. II. • S 
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will evidently- not be satisfied unless ) -j-r ; be infinite, 

) 0(a;, u) 
We infer then that this is a neceasary conditioQ in order that 
« = may be a particular integral. 

This is Euler's fundamental tbeorem, and from this, by 
means of an hypothesis agreeing with that of Poiason con- 
cerning the forta of the transformed differential e<iuation, he 
arrives at the condition 



[In the passage to which Professor Boole refers, Euler 
does not undertake to discuss the nature of any solution, 
but only of a solution of the form at = constant. On his 
page 408 Euier proceeds to discuss the nature of any solu- 
tion. Professor Boole seems to me to attribute too much 
to Euler. For the convenience of those who wish to ex- 
amine the original, I will give the reference to the passages 
in the later editions of Euler s Institutiones Calculi Integralis .• 
Vol. I. pages 343 and 355 of the edition of 1792; Vol, 1. 
pages 342 and 354 of the edition of 1824.] 

Laplace in the Memoirs of the French Academy for 1772, 



d^y 



p. 343, established the t 

dp_ 

and shewed their respective uses. He established also the 
test wliich consists in tbe comparison of differentia! coefficients, 
and he supposes it universal. He adopts the hypothesis of 
his predecessors as to the forma of expansion, but with some 
recognition of its insufficiency. 

Lagrange in the Memoirs of the Academy of Berlia for 
1774, p. 197, and 1779, p. 121, appears first to have developed 
the theory of singular solutions in its two forms of derivation 
from the complete primitive and derivation from the difieren- 
lial equation, and to have established the essential connexion 
of these. But supposing the difi'erential equation to be ex- 
pressible in the rational Ibrm 

Fix,y,p)=0, 
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and employing the differential coefficients of F[x,y,p) in- 1 
fltead oftboae of p he was led to sacrifice rigour to symme- 
try. One of his results has often since been adopted as a ■ 
■ test of singular solutions. It may be thus stated. J 

H Prop, a Bmgular aolutloa makes the general value of I 

-j^ , deduced from the differential equation in its rational and 1 

integral expression, to assume the form - . 

[The demonstration is given in Chap. vill. Art. 14.] 

This ambiguity of value of --rr, is evidently but au expres- 

flion of the fact that the contact of a curve of the complete 
primitive and that of the singular solution is not in general i 

of the second order. 

The result given in equation (5) of Chap. yiii. Art, 14 has 
also been adopted as the test of singular solutions. 

The researches of Poisson and Cauchy havo already been 
noticed. It is certainly remarkable that the final test to 
which Cauchy's analysis led should be essentially the same as 
that which had been discovered by Euler bo long before. 

Professor De Morgan has thrown an important light upon 

»the natm-e of the conditions 
^P = ^ '-^ = 00 
dy ' dj; ' 



which are fulfilled by all singular solutions in the expression 
of whicli X and y are both involved. He has shewn that any 
relation between x and y which satisfies these conditions will 



from the differential equation, infinite ; that it may satisfy the 
differential equation even if it make -r-^ infinite ; lastly, that 

3—2 



1 
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if it do not satisfy the differential equation, the curve it 
represents is a locus of points of infinite curvature, usually 
cusps, in the curves of complete primitives. 



The proof is as follows : 
Let p = (l> {x, y) 

be the differential equation. Then the proposed conditions 
are 

#(^>y) ^^ # (a?> y) ^ CO 
dy ^ dx ' 

therefore by differentiation, 

dxdy dy^ dx ' dx*^ dxdy dx ' 

whence we have 

d"^ d^ 

dy __ dxdy ^ daf 

dx d^<^ "" d'^<i> 

dy^ dx dy 

These are two equivalent expressions for the same value of 
-,- . The question now is, under what circumstances this- 

value of -^ will satisfy the differential equation. 

Now from that equation we have by differentiation 

d^y _ d<j> d<j} dy 
da? dx dy dx^ 

whence 



dy _ 


d^y d<f> 
dx^ dx 


dx~ 


d(l> 
dy 
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If then -j-^ be finite we have, since -^ and -^ are both 
(IX ax ay 

infinite, 

^ ' 

dy __ dx 

dx d<f)\ 

dy 

and this by the rule for the evaluation of fractions of the 
form ~ is equivalent to the value in either of its forms before 

obtained for -~. Hence, any relation which satisfies the 

d\ 
given conditions and makes -y^ finite, will satisfy the difie- 

rential equation. 

And the same result holds even if -j4 be infinite, provided 

that ^ -f- -7^ vanish. 
dxr ay 

Lastly, as when this result does not hold, the failure is due 

d^v 
to the infinite value of -t4 , we see that the line in which the 

locus of the proposed relation intersects the curves of primi- 
tives will be a locus of their points of infinite curvature. 

[Transactions of the Cambridge Philosophical Society, 
Vol. IX. Part II.] 

Legendre's Memoir of 1790 throws but little light upon 
the subject of this Chapter. But it exhibits the theory of 
the singular solutions of diflferential equations of the higher 
orders, both ordinary and partial, in a form of great beauty, 
and will be noticed in the proper places. 
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CHAPTEE XXIL 



ADDITIONS TO CHAPTEE IX. 



1. By successiye application of the second theorem of 
Chap. IX. Art. 13, a linear equation of the n^ order may be 
reduced to one of the (ti — r)"* order, if r distinct integrals of 
what the given equation depriyed of its second term would 
be are known. 

The reduction may howeyer be effected immediately by 
the method of the variation of parameters. In this and in 
most general investigations connected with differential equa- 
tions great advantages in point of brevity and of the power of 
expression are gained by the employment of the symbol of 
summation 2, and of the language of determinants. I shall 
exemplify this here. 

Suppose the given equation to be 

d^^'^^'d^''^^'^^ +^n--5C (1), 

and let y^, y^,.*>yr ^^ ^ particular values of y, satisfying the 
equation 

d*i/ J dr^y J dJ^'^y /i — a 
^■*'^^^"^^»^ +^«-0 (2). 

Thus y = c^y^ + c^, ... + c^r 

is a solution of the latter equation including these particular 
solutions. We shall represent this by 

y^tc^i (3), 

and regarding the quantities c,, c„... <v, represented here by 
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Cj as variable parameters, shall seek to determine them so that 
the above value of y may satisfy the equation given. 

These r parameters, enabling us to satisfy r — 1 arbitrary 
conditions, besides satisfying the differential equation, we may 
choose these so that 

may be the same inform as if Cj, c,, . . . c, were constant. Now 
from (3) 

whence 

Tx-^'di* 

provided that the condition 

be satisfied. Differentiating the first of these equations, we 
find in the same way that 

provided that the condition 

^^^^^ 
dx dx 

be satisfied. And thus continuing we see that the system of 
r equations 

^ dy ^ dyi d^^y ^^ dr% , ,n 

will hold true provided that the r — 1 conditions 
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be satisfied. In each of these equations the symbol X is to 
be interpreted by giving to i the successive values 1, 2,... r, 
and taking the sum of the results. 

Differentiating the last of the equations (4), we have 

^ = tc-^'-ht ^'^'^^ — * . 
rfa?*" * rfic** dx""'^ dx * 

As we cannot impose the condition that the last term of 
this equation shall vanish, let z represent its unknown value, 
then 

d^^-^'^dx^^'^"' W- 

Now the system of equations (5), together with 

^d^dci_ 
^ dx'-' dx ' 

constitute a system of r simple algebraic equations deter- 
mining by solution the r quantities 

dc^ dc^ dCf 
dx^ dx '"' dx 

in terms of their coefficients and of z, and therefore in teniis 
of x and z^ since the coefficients are known as functions of 
X. It is evident also that as the second members of all the 
equations but one vanish, and the second member of that is 
Zy the values so determined will be of the form 

XjjXj,... X^ being known functions of x. Thus the r un- 

dc dc 
known quantities ^~>"-3^ ^^^ made to depend upon only 

one unknown quantity, viz. z. It remains then to deter- 
mine z. 

For this purpose we must complete the expression of the 
differential coefficients of y, and substitute in the given dif- 
ferential equation, and then seek to satisfy that equation. 
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Now differentiating (6) we have 

dx^^'^'^^'Ux^^^^ djJ dx^dx 

^^""'d^^^^W^j'^dx 

dc 
on substituting for -j^ the value XiZ as above determined. 

We observe that the coefficient of z is here a known function 

of X, If we differentiate this equation and in the result sub- 

dc 
stitute as above for -j^ , we shall have a result of the form 

d^'v ^ d^yi , ^ . ^,dz d\ 

L and M being known functions of a?. Ultimately then 
we have 






ft— r„ 
Z 



dx"" 'dx'' ^ dx dx''-'" 

Thus, while y and the differential coefficients of y up to the 
(|. — l)th j^^g Qf ^jjQ same form as if c, , Cg,... c^ were constant, 
tlie succeeding ones differ in containing an additional portion 
consisting of z, and differential coefficients of z multiplied by 
known functions of x. The result of substitution of these 
values in the given differential equation will therefore consist 
also of two classes of terms, viz. terms under the sign of 
summation, which will be the same in form as if c^, c^^,,,c^ 
were constant, and terms involving the differential coefficients 
of z up to the {n — r)***, with multipliers which are known 
functions of x. We shall in fact have 

^'^^"^^^^^"^^-^^ ^^\ 

-^l^^-^d^r^' + ^=x 
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Now yi being by hypothesis an integral of (2), the first 
line of the above equation vanishes, and there remains the 
linear equation of the (n — r)"* order 

^^ + ^S^^ + ^ = X 

Supposing z hence determined, we have in general 

Ci=^ iXizdxy 
and hence 

y^yAX^zdx + yAx^zdx +i/JX^dx, 

and as z will have n — r distinct values, each involving an ar- 
bitrary constant, the above equation will furnish n — r distinct 
values of y, each involving an arbitrary constant. It is to be 
observed that no arbitrary constant need be added in the inte- 
gration of the terms Xizdx, for the effect of such addition 
would only be to reproduce the known integrals c»yi. In 
this way, however, the equation would represent the general 
integral of the differential equation given. 

2. Let us examine the form of the result in the particular 
case in which r = w — 1 . 

Here we have 

dx"^ ^ * dx"^ 
from 7n = Oto7n = w— 2, then 

d^^-^'^-^^^K^'d^^J'^di' 
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Accordingly the diflferential equation for z will be 

|+s(x,p).+^..-x.. („. 



Now the equations for determining 

dc^ dc„ dc, 



• • • 



dx ' dx '*** dx 

dci 

become on putting X^ for -^ , and writing for brevity y', for 

^ W'for^' &c 
^,yi lor ^ , etc., 

Vi^i. + y»-^» • • • + Vn-t^i = o> 






Whence, by the theory of determinants, 






■» Jf dy.t"^ 



if standing for the determinant 



'-^~Mdy, (■")' 



'*-! 



Viy y^r 



y 



«-l 
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Now the determinant is ultimately a function of x ; and 
such indeed that 



^ / d«-V,\ 1 dM 



For 

dM^ ^ dM dy, ^dMdyl_ ^ dM dyj""^^ 
dx dyi dx dyl dx '" dyi"""^^ dx 

^ dM , , ^ dM „ , ^ dM ,„_!. . . 

Now M being homogeneous and of the first degree with 
respect to the quantities y^, y^, y^j, we have 

^dM ,-. 
2 -J- yi = M. 
dyi ^* 

Hence 2 -j— yl is what M becomes when in its expression 

.Vi, 3^a, ...yn_i are changed into y/, y^', ... y^^, therefore it is 
what M becomes when two of its rows of elements become 
identical ; therefore it vanishes. In like manner all the other 
sums in (8) vanish excepting the last, iox y^''~^\,..,y^^'^^^ is 
not a row of elements of the determinant M. Thus we have 

^ d^ ,,(n-i)^dM 
^dyf"^^' '^ dx' 

Hence 

^X,yl )=2^^p,2^/ '^Mdx' 

Thus the equation (7) becomes 

dz . ' / 1 dM . . N ^. 

therefore 2 = ^e'-^"'^ JMe-^^^" Xdx. 



ART. 2.] ADDITIONS TO CHAPTER IX, 45 

Hence, since 

^_ Y — -L ^^ 

whence c,=jj^^^zdx, 

we have . y = %,|1 ^^^ ^e&, 

2 being given above. 

In the case of X= 0, we have 






whence 



y-M^w^'"'-''- 
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1. The theory of singular solutions of differential equations 
of the higher orders has been presented in the most complete 
form which it has yet received by Legendre. (Mimoires de 
r AcadSmie Boyale des Sciences^ 1790, p. 218.) He determines 
first the possible forms of these solutions considered as emerg- 
ing from the complete primitive by the variations of its arbi- 
trary constants, and secondly the theory of their derivation 
from the differential equation itself. 1 shall follow the same 
order, and shall in the end endeavour to point out in what 
respect Legendre's theory may be regarded as complete, and 
in what respect it is imperfect. 

Suppose the differential equation to be of the n^ order, 
and let it when solved with respect to the highest differential 
coefficient of y be represented by 

y«=<^(aj,y,yi,y„...y*_J (i), 

in which, for brevity, 

_dy _d^y _ d^ 

Let also its complete primitive, solved with respect to y, 
be represented by 

y =/(iC, a^, o„ ... a^ (2), 

ttj, a,, ... a,, being the arbitrary constants of the solution. If 
we differentiate (2) with respect to a?, regarding «i, a,, ... o^ 
no longer as constants but as functions of a?, so to be deter- 
mined as to leave the expressions for yj^, y^, ... y» as functions 
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of a-, a^, ... a« the same as before, we shall have, on repre- 
senting the second member of (2) by/, 



whence 



^^ dx da^ dx da^ dx da^ dx ' 



^^ dx' 



provided that 



df_da^^d£da^^ ,+_^^=0. 
da^ dx da^ dx da^ dx 

Differentiating on the same hypothesis the first of these two 
equations, we find in the same way 






provided that 



^ ^i + _i^ ^3 + _iZ_-- = o. 



da. 



dxda^ dx dxda^ dx 



dxda^ dx 



And continuing thus, it results that the system 

_df _^ _dy 

y^^dx' y^^ dx^'-'-^y^^ dx'' 



(3), 



-will be satisfied, i.e., that ^1,^3,... y^ ^^^ have the same ex- 
pressions when a^, a^^...a^ are variable as they have when 
these are constant, provided that the law of their variation be 
determined by the conditions 



df_da^^ da^ 
da^ dx da^ dx 



+ 



^f ^^n_n1 



da., dx 



= 



d'f da, ^ d^f da, ^ d^f da^^^ 

da^ dx dx da^ dx dx da^ dx dx 



^...(4). 



d^f da, d^f da^ 

da^ dx^'^ dx da^ da?*"^ dx 



^f da^ ^ 
da^dx"" ^ dx 
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In this system the coefficients of 

da^ da^ da^ 

dx ^ dx^ *" dx 

are known functions of a;, a^, a^, . . . a^ when the form of/ is 
known. 

Eliminating 

da^ da^ do„ 



dx^ dx^ "* dx ^ 



we have a relation between x, a^, a^, . • . a^; and this relation, 
with the given complete primitive and the first n — 1 of the 
derived and reduced equations, viz., with 

_/ __rf/ _dy _ d^'\f 

will enable us to eliminate a^^ a^^ , . , a^, and to obtain a rela- 
tion of the form 

^v^y^ dx' rfaj^'-'Vx^'v ^^^• 

This is a differential equation of the (n — 1)*** order. It dif- 
fers in its origin from the given differential equation, in that 
a new relation between ar, a^, eZg, . . . a^ has been employed in 
place of the n^ equation, derived by difi'erentiation from the 
complete primitive, for the elimination of the constants. 

• 

The differential equation of the {n — 1)**^ order thus obtained 
has an integral expressing y in terms of a*, and w — 1 arbi- 
trary constants. This is the most general form of a singular 
solution of the differential equation. 

It is possible that the elimination of a^, a,,, . .• a. may 
lead to a resulting difierential equation which, instead of 
being of the order n — 1, is of the order w — 2, w — 3, &c. 
The complete integral of such equation would be a singular 
solution of the differential equation. These possible types of 
solutions are distinguished by Legendre according to the 
number of arbitrary constants which they contain. A solu- 



H^-j 
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tion containing m — 1 arbitrary constants is called hy Lim a 
Blngnlai solution of the first order ; one containing » — 2 ar- 
bitrary constants a singular solution of the second order, and 
80 on. 

Adopting this language we might term the complete priml- 
tire a singular solution of the oi-der 0. 

Lastly, any relation between x and y, which satisfies the 
given differential equation, will constitute a particular case, 
■ either of the complete primitive or of one of the general 
forms of singular solutions above defined. In the case of 
diiFerential equations of the first order it is seen that no arbi- 
trary constant can appear in the expression of the singular 
solution. 

■ Ex. The equation 

has for its complete primitive 

y=~ +hx+a^-i-F (6), 

Liequired its singular solution. 

L Proceeding as above, we find on the hypothesis of a and b 
wng variable parameters, the same foi-mal expressions for 

I jj ^ if those parameters were constant, viz. 



■■('), 



provided that the variation of a and h be such as to satisfy 
the conditions 



I 



(?"«)!-(- ^')|-1 


(3) 


f B. D.E. ir. 
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Eliminating hence -j- and — , we liave 

a? 
2a - - - 2&a; = 0. 

And from this, the complete primitive, and the first of the 
derived equations (7) eliminating a and 6, we find 



(I)"- 



This 13 the differential equation of the first order, by the 
solution of which the most general form of the singular aoltt- 
tions of the given differential equation will be determined. 

Reducing it to the form 

and integrating, we find 

{16!/ + ix^ + x*)^ =x[l + x')^ + log \x + V(l + aJ")! + G, 

This then is the general expression for the singular aofa- 
tions of the given differential equation. We see tliat it in- 
volves in its expression one arbitrary constant. 

The differential equation (9) may properly be termed a aiib- 
gular first integi-al of the given differential equation. The 
singular first integral (9) has itself abo a singular eolutioD, 



hut this 13 not a solution of the original differential 
Nor have we any right to expect that it should be so. A 
singular solution of a dift'erential equation of the first order 
does not necessarily.satisfy the differential equations of higher 
orders derived from that equation, Chapter xxii. Art. 7. 

2. It remains to establish the theory of the derivation of 
the singular solution from the differential equation withont 
the mediation of the complete primitive. 



« 
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Resuming the differential equation in its reduced form (1), 
and representing its second member by tf>, suppose an infini- 
tesimal variation given to the arbitrary constants of its com- 
plete primitive, and let the symbol 8 be used to denote the 
corresponding derived variations of y, y^, ...y^. Then we 
have 



'•-1 



and so on. Hence, substituting and transposing, 

'd^'^d^.l^^d^l^ '^''•"^ ^^"^• 

> Let us consider the real nature of this equation. 

If a value of y, suppose y — '^ix), satisfy the given differ- 
ential equation, that value substituted in the coefficients 

dy^^' dy, ' 



'f»-a 



of the above equation will convert them into functions of a:, 
and the equation itself will become a linear differential equa- 
tion, the solution of which will determine hy as a function of 
35. If the differential equation (10) be really, as it is appa- 
rently, of the w*^ degree, hy will have n arbitrary constants, 
Qj, • • • Om, and will be of the form 

■Pj, jFJ,, . . . P^ being functions of x. Hence 

y + Sy = f(aj) + ajP,... + a,P,. 

We gee thus 'that the given solution y = -^ (a?) will, l^' a 
Piirticular case of this general integral involving n coi:^tants. 
h will therefore be a particular integral of the proposed. 

4—^ 
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If, owing to the constitution of its coefficients, the difierential 
equation (10) be of the degree n—\, we shall have 

and ?/ = 1^ (x) will then be a particular case of a eolution 
involving n — 1 arbitrary conatanta. It will therefore be a 
singular solution of t!ie first order. Even ao, if the differen- 
tial equation (10) be of the degree n — 2, y — ^{x) will be a 
singular solution of the second order. And generally, if 
the differential equation be of the r"" degree, y = '^{x) wUl he 
a singular solution of the order n — r. 

Resuming the equation (10) it is evident that it cannot 

be of the degree n—l, unless j— ^ he infinite. For, dividing 

by ■ , , we have 



j_.^i'_^::%_&e.=o 

1^0 dx" dx' ' ' 



d^ 



in which the first term does not vanish unless -r~- be in&- 

nite. This then is the necessary condition for a singular 
solution of the first order. For one of the second order we 
must have in like manner 



'^y^x 



It follows hence that to find the singular solutions of a 
differential equation of the r^ order, we ought to differentiate 

the eqnation, regarding y, -J , -r^ , &c. as varying through 

the variation of the arbitrary constants, to form in this way 
a linear differential equation for Sy, to examine the conditions 
under which this equation reduces to the (m — 1 )*, or to a lower 
degree, and to examine whether the most general relation be-. 
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tween x and y which aatiafies such condition, satisfies also the 
given differential equation. If bo it may be regarded aa a 
singular solution. 



Besiiming the last Kxample, viz. 

«-,^ + l™>^_('W_('^_,W_n 



and operating with S we have 



2^ 'd^^^'W "dx'/jda? 



-hKI-3)}° 



which reduces to a linear differential equation of the first 
order for determining Sy, provided that we have 

Elimmating -j^, from the given equation by means of this 
there results 

(l)"M^#I-(' + »^^-fe = »' 

and we find on differentiating this that it does constitute a 
solution of the given equation. It is therefore a singular 
first integral of that equation. We see that it agrees with 
the result obtained under the same name in the previous 
Article, and the rest of the solution need not be repeated. 

3. Upon Legendre's theory, and upon its results, the fol- 
lowing observations may be made. 

Ist. We leaVn from it that there may exist different 
general forms of the solution of a differential equation of the 
f** order, viz. the complete primitive invplving n arbitrary 
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conslant8j and general forms of singular solutions containing 
fewer than n arbitrary constants. A solution y = V'(^) °^ 
unknown origin being given, we construct a differential equa- 
tion for determining Sy, and, solving it, form tlie expresaion 
for y+iy, and from the number of infinitesimal arbitrary 
constants it contains, determine tbe nature of that general 
value of y of which tlie given value is a particular case. 
Now we are not to infer from this that the form o{y+ By will 
be the same as the general value of y in question. But we 
may infer that it will be a form to which that general value 
is reducible. And the actual reduction will be effected by 
expressing the general solution (as is always possible) in a 
form permitting its expansion in ascending powers of the 
arbitrary constants, and in the expansion making these con- 
stants infinitesimal, and rejecting all powers of them above 
the first. In fact, if 

J, -/(«;, .„<.., ...»,) 

be any general form of solution which, when we assign to 
a,, o,,...a, particular values (e.g. make them vanish) re- 
duces to 

then we shall have 

the brackets denoting that after differentiation we make 
a,, a,, ...a, vanish. 

This is that limiting form of the solution which Legendre'a 
method enables us to construct by the solution of a linear dif- 
I'erential eqtiation ; and the ground of the sufficiency of hia 
method consists in this, that the infinitesimal q^oantities 

Sa,, Ba^, ... Ba,, 

which are io fact the arbitrary constants of that solution, are 
e'jual in number to the arbitrary constants of the general 
unlimited Bolation, the nature of wliidi is thus made known. 
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2ndly, Legendre's tests for differential equations of the 
higher orders are in kind and effect analogous to the tests 

dp ^ d 1 ^ 

eZy"" ' dxp" 

for differential equations of the first order. They enable us 
to decide whether a solution possesses singularity, not whether 
it possesses the envelope species of singularity. The comple- 
tion of Legendre's theory would consist in the discovery of 
those further tests dependent upon integration which corre- 
spond to the test of JSuIer and Cauchy for differential equar- 
tions of the first order. 
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CHAPTEK XXIV, 



ADDITIONS TO CHAPTEE XIV. 



[Art. 1 was intended to follow Chap. xiv. Art. 2.] 

1. A8 the condition of dependence of functions of two 
variables is of fundamental importance in connexion witk the 
tlieory of ordinary differential equations, so the generalized 
condition of dependence of functions of any number of vari- 
ables forms a fundamental part of the theory of partial differ- 
ential equations. This ia contained in the following proposi- 
tion. 

Prop. L Km,, u„...u„ are functions of x^, a;,,. ..a:,, 
but are as such so related that some one of them is expressi- 
ble as a function of the others, or more generally that there 
exists among them some identical equation of the form 

F(u„u,,...u.)=0, (I), 

80 that as functions of a;,, a;,, ... a;, they are not mutually 
independent, then, adopting the notation of determinants, the 
condition 



du, 


du, 


du. 


d., 


d^r 


-d.. 






du. 


dx. 


da:,' 


ix. 


d". 


du. 


du. 


d., 


■^r 


'"dx. 



is identically satisfied. Conversely, if the ^bove condition be 
identically satisfied, the functions n^, u^, ...u, are not mata- 
ally independent in the sense above explained. 



ART. 1.] 



ADDITIONS TO CHAPTER XIV. 



57 



First let it be noticed that the Proposition is but a general- 
ization of that of Chap. ii. Supposing U and u to be two 
functions of a? and y, the condition of their dependence is 
affirmed to be 

dU dU 



= 0, 



dx ' dy 

du du 

\ dx^ dy 

i. e. it is the result of eliminating da?, Jy, from the equations 

dU. dUj ^ 
-d^^-dy^O, 



du J du ^ 



and therefore it is 



dU du dUdu_^ 
dx dy dy dx^ ^ 

as expressed in Chap. ii. 

We proceed to the general demonstration. 

Let the first member of (1), considered as a ftmction of 
UyU^f...u^ be represented for brevity by F; then differen- 
tiating, we have 

dF^ dF. ■ dF ^ ^ 

-T- du^ + --7— au^ + ... + -T— du^ = 0, 

from which it follows that if du^^ du^y ... du^^^^ are equal to 0, 
then is du^ equal to ; or, since u^jU^^^u^ are fiinctions of 

du. y . du. , . du. , 



da?j * dx^ "'" * dx^ 



dx. 



dx. 



du, 



«»-t 



€lx^ 



da>, + ^dx,...+^dx,^0 



dx^ 



dx. 



(3), 
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then is 

s."^+a?,*'-+sr. '*'■=» '<! 

Thus the last n ec[iiationa, linear with respect to 

dx^, dx^, ... dx^, 

are not independent, and therefore by the theory of linear 
equations the determinant of the system vanishes identically. 
Now this is expressed by the condition (2), 

It remains to prove the converse, viz. that if the condition 
(2) be identically satisfied, the functions w,, m,, .■■ m„ will not 
be mutually independent. 

First, the n — 1 functions m , w,, ... w^ are either mutually 
independent or not mutually independent. 

If not, then the n functions w,, « , ... m„ are not mutually 
independent, and the Proposition to oe proved is granted. 

If M,, «j, ... !*,^j are raatually independent as functions of 
a;,, x^, ... ar,, they may be made to take the place of w — 1 of 
these quantities, e.g. x^, a;,,.., x^^ in the expressing of u,, 
i.e. we may, by means of the expressions for «,, «,, ... w,^,, 
eliminate from that of «„ the quantities x,, a;,, ... ic^,, and ao 
express u„ as a function of m,, m,, ... m„_j and a;,. Suppose 
this done, then the system (3), (4) will be converted into 

(?Mj = 0, dui = 0, (iu^, = 0, 

(?M, , _ du, , du, , du, , 

-i— dw, + J— flw. . .. + -j— du., + J— aa 
du^ ' du^ * aWn-i dx„ 

Now, the determinant (2) vanishing, the equations of the 
linear system (3), (4) are not independent ; therefore those of 
the transformed system, as written above, are not independ- 
ent; therefore the last equation of that system must be a 
consequence of the others which manifestly are independent. 
Bat from the form of that last equation we see that such can- 
not be the case unless we have 
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■which implies that u, is a function of «,, u,, ... k^, merely. 
Hence the fuuctions u,, m,, ... w, are not independent, as was 

to be shewn, 

The firat memher of the equation of condition (2) is com- 
lonly called the functional determinant of u,, m^,, ... «, with 
respect to ar^, a;,, ...a^,. The proposition may therefore be 
expressed as follows. 

The condition of dependence or independence of any sys- 
tem of functions of as many variables is the vanishing or 
non-vanishing of the functional determinant of the system. 

On account of the great importance of this proposition it is 
desirable to illustrate it by an example. 

£z. Are the functions 

x + 2y + s, x — 2i/-r3z, 'ixy -xz-i- iyz - 2z' 
mutually independent or not? 
The equation of condition is 

I, 2, 11 



Sy-s, 2a: + 4s, -31-1-4^-43 | 
that is, 

4 (- ic + 4y - 42) + 8 (2y - s) - 2 (2a: -1- 4a) = 0, 

"which is identically satisfied. Hence the functions are de- 
pendent. In fact, representing them by w, v, «t, we have 

iw = u' — n", 
[Art. 2 was intended to follow Chap. xiv. Art. 4.] 
2. As it has been shewn that a primitive 

«-*H (1) 

leads to a linear partial differential equation of the form 

Pp+.Qq^B (2), 
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provided that u = a, w = 6, are integrals of the syatem of ordi- 
nary differential equations 

dx dy dz ,.■. 

T-'Q^R ^^'' 

it is evident that we shall obtain a aolution of the partial dif- 
ferential equation (2) by constructing the system of ordinaiy 
differential equations (3), deducing their general integrals 



and then constructing from these the primitive (1). 

But the question arises, Will this be the most general solu- 
tion of the partial differential equation given ? 

That it will be so, may be shewn by means of the general 
proposition. See Art. 1. 

For let M = represent any aolution whatever of the given 
partial differential equation. Differentiating this with respect 
to X and y, we have 

dw dw dv) dw _ 

dx dz" ' dy dz ' 

substituting the values of p and q formed from this in the 
given equation, we have 

dx dy dz 

which must be identically satisfied. 

In like manner, u=a, v = h being solutions of the same 
equation, we find 



which must be identically satisfied. 



rial 
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Eliminating P, Q, S from these three equations, it results 
that the fiiTtctional determinant ofw, u, v, with respect to x, y, z, 
will identically vanish. Hence k; is a function of v. and v, 
and the equation w = is a particular case of 

-F{.,.)-0, 

which ia thus shewn to be the general Integral of the given 
equation. 

We aie thus led to the following general Bule. 

BuLE, To integrate the equation I^+Qq=Il we must 
Jbrm the system of ordinary differential equations 

dx dy _6x 

deduce their general integrals in the form 

u — a, v=b, 

and construct t?ie equation 

F{u, v) = 0. 

Thia will be the general solution sought. 

[Art. 3 was intended to follow Chap. xiv. Art. 5.] 

3. The above theory may be extended to linear partial dif- 
ferential equations of the first order, without regard to the 
number of the variables. 

First, the theory of the genesis of such equations ia ex- 
pressed in the following proposition. 

PeoP. a primitive equation of the form 

F{u^,u„...u„)=-Q (1), 

in which «,,«,,...?/, ate any given functions of the vari- 
ables z, dependent, and x , a:,, . , . a;„ independent, will satisfy 
the linear partial differential equation obtained by eliminating 
dz, dXi, ax„...dx„hoia 
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expressed as total differential eqnationa witli respect to the 
primitive variables, and the equation 

dz —Pidx^ —pjdx^ ... —p^dx^ = Q. 

Of this important proposition I propose to pve two distinct 
proofs. 

Ist proof. Forming the total differential of the given 
equation we have, on representing its first member by F, 

dF, , dF, . dF, „ 

flu, ' AM, ' au„ 

Xow this cannot be tme for all forma of the function unless 
we have the separate conditions 

dM, = 0, dw, = 0, dw, = 0, 

Strictly to prove this, suppose F,, -f,, . . . F, to be any n 
distinct and independent tunctior.g of u^, u,,...u„, and as 
such, distinct and independent forma of F. Then the above 
equation gives 

dF^ , dF,, dF, , 

J—* du, + T— ' du, ... + T— i du. = 0, 

(IM, * dMj ^ du, " ' 

dF,, dF,, dF,^ 

du ' du ' du 



1 



du, ' dtt. 



= 0. 



L 



Now F,, F^,...F^ being independent, their functional 
detemiinant with respect to «,,«,,...«,, does not vanish. 
This again is the condition necessary and sufficient that the 
above system of linear equations may be independent ; and 
this lastly being the case, their only possible solution will be 



du, = 0, d«, = 0, d«„: 

as was to be shewn. 



mt 



• 
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These equations in their developed expression 
du. y du. y du. , du. , 

pdx, + pdz = 0. 

dx^ * dz 

^dx, +^& = 0, 

ctej * dz 

enable us to determine the ratios of dx^^ dx^^ • . . dx^, dz in 
the form 

dx^_dx^ _dx„_dz , . 

X, " ^ '" "" X, ~ jB ^"^^ 

where X^, X,, ...X„, i2 are functions of the original vari- 
ables. And now, forming the equation 

Pidx^ '\-p^dx^ ... +pndx^ — £& = 0, 

and eliminating the' differentials, we find 

X^p^+X^p^...+X^p^ = B 
for the partial differential equation sought* 

2nd proof. Differentiating the given primitive with respect 
to a?j, as contained explicitly in the functions w, , t^,, . . . w„, and 
also implicitly in the same through Zy we have, on represent- 
ing the first member of the equation by F^ 

dF 

dF /du 
du^ \dx 
or 



l/^^_l ^^ dF (du^ duA 

h\^i ^' dz J '^ dir,[dx, ^^' dz) - 



+^i^;+i..^)=<>> 



Since 



dFdM^dF^da^ .dF^dF ^^ 
du^ dx^ du^ dx^ '" du^ dx^ dz ^^ ' 

dFdu^dFda^ .dFdu^^dF 
du^ dz du^ dz *** du^ dz dz 
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du^ , dF ^ 



Differentiating thus with respect to the remsdning inde* 
pendent variables, we obtain finally the system 

du^ dx^ du^ dx^'" du^ dx^ ' dz 

dl^da dF^du^ .§^^,§E. =0 
du^ dx^ du^ dx^'" du^ dx^ dz** ' 

^dFdu^dFdu^ dF ^ dF 
du^ dx^ du^ dx^'** du^ dx^ dz^"^ ' 

from which,^ in combination with the equation 

dF du^ dF du^ dF du^ dF ^ 

we can eliminate 

dF dF 
du^^ du^^ 

The result will be 



•••••• 



4F dF 

du^ ' dz * 



du^ du^ 
dx^ ' dx^*' 



du 
dx^ 



>P1 



du^ du^ 
dx^' dxj 

du^ du^ 
dz ' ^" 



dx' 



P* 



du^ 
dz 



,-1 



=0, 



or, conrertiDg rows into colnmns, 



du^ 
dx^' 


du^ 

dx;'" 


du^ 
-dx,' 


du^ 
dz 


du^ 
dx^' 


du^ 
dx,'"' 


du^ 


du^ 
dz 


Vx^ 


Pj,... 


— jPn> 


-1 



= 0. 
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which is the determinant form of the result affirmed in the 
proposition. 

The second of the above forms of demonstration seems to 
be preferable to the first, in that it rests only upon the consi- 
deration of the one general form of the function F. I have^ 
however, given the two proofs, chiefly in order to illustrate 
an important remark, viz. that, in nearly all general re- 
searches connected with partial differential equations of the 
first order, two modes of procedure, the one involving the 
use of differentials, the other that of differential coefficients, 
may be employed, and that between the forms to which these 
respective modes give rise, a certain law of reciprocity will be 
found to exist. 

The theory of the solution of the partial differential equa- 
tion 

follows immediately from that of its genesis. If we repre- 
sent by 

the integrals of the system of ordinary differential equations 
(2) a solution of the given partial differential equation will 
be represented by (1). That this will be also the most gene- 
ral solution may be shewn by the argument of Art; 1. For 
if w? = represent any solution, then since 

dw dw _ ^^ , ^ «- A 

we find 

from which, in combination with the corresponding equations, 
y^du y. du^ J. ir ^^1 J- 7?^^i -n 



Y^^n^y^ dUn i Y ^^n , -pdu^ 



B. D.E. II. 
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eliminating X^, A'„ ... -V., E we obtain a result which ex- 
presses that the functional determinant of m, Wj^, ... w, with 
respect to the original variables is virtually 0. Whence v> ia 
a function of «,, «,, ... «,, and the proposed solution is in- 
cluded in the one to which the above method of solntion 
leads. 

That method may therefore be stated in the following Knle. 

KtJLE. To integrate the linear partial differential equaliim 





X,|l+X.|!.. 


+ X 


dz 
dxj^ 


= It 


rmthe 


system ofordiwi 
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equations 
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dz 
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«. = «.. 


..«. 


= «., 






•F(« 


„«„■.■ 
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[The general observations were intended to follow Chap. 
XIV. Art. 6.] 

General obaervationa. 

4. The relation which exists between a proposed linear 
partial differential equation and its auxiliary system of ordi- 
nary differential equations should he carefully studied. While 
it ia proper to say aa above that the general integral of the 
one requires the knowledge of all the integrals of the other, 
it is also proper to describe that general integral simply as 
the most general form under which an integral of the aaxi- 
liitry system can appear. If 

«, = «,, w, = a,, ... M„ = a, 
are integi-ala of that system, then 

jF(«„W„...U„)=^ 
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is the one general form of an integral of that system, and 
due regard being had to the arbitrariness of F, this is equi- 
valent to 

^^(Wj, Wj, ... wj = 0. 

5. The form which the auxiliary system assumes when 
the given partial differential equation is deficient in any of its 
terms should be noticed. 

If ^ ss 0, the auxiliary equation 

becomes, on clearing of fractions, 

(&?! = 0. 
And thus, if X^, -Xj, ... X^ vanisli, the given equation being 

the auxiliary system will be 

diCjSsO, dir, = 0,..*d!a?,= 0, 

dXf^^ __ dXf^ _ dx^ _ dz 
X^4i Xy+a X^ X 

and the integrals of this system being of the form 
the general solution of the given equation will be 

F{X^, aJo Wr+i, Wn) =0. 

This conclusion would follow also from the principle laid 
down in Chap. xiv. Art. 2. 

Linear partial differential equations in which the absolute 
term is wanting, and which are therefore of the form 

■rr dz -XT dz -rr dz 

5— a 
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may be termed homogeneous. As in this case one of the 
auxiliary equations is 

rf^? = 0, 

the general integral will be 

Wj, Wj, w^i being found by the integration of the remain- 
ing auxiliary equations 

"X'x — X* 

When X,, X,, X), do not contain z, the solution is Iwst 

exhibited in the form 

6. Every linear partial differential equation can be converted 
into a homogeneous one containing one additional variable. 
For it is shewn in Art. 3, that if w = be any integral of 

•y dz yr dz i T- ^^ ^ V 

then is 

X— 4-X— -4-X ^4-T— =0 

a homogeneous equation with a new variable. 

From the general integral of this equation, that of the 
former one may be deduced by making u = 0. 

7. The solution of partial differential equations is some- 
times facilitated by introducing a new system of independent 
variables. The actual transformation is greatly facilitated 
by the following symbolical theorem. 

Theobem. If the partial differential equation 

Y dz y. dz , xr ^-^ v 
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be expressed symbolically in the form 
in whicli 

then, if y,, ^a, ^n^^ ^ ^^^ system of independent vari- 
ables given in expression as functions of the old ones, the 
transformed equation will be 

(%,)|+(AyJ,|...H-(Ayj|-=X 

For, regarding ^s as a function of y,, y., y,, we have 

dz _ dz^ ^ dz^ dy^ dz dy^ 

dXy^ dy^ dx^ dy^ dx^ "* dy^ dx^ ' 



dz _ dz dy^ dz dy^ dz dy^^ 

^n, " dy^ dx^ dy^ cte^ *" dy^ dx^ ' 

whence, substituting in the given equation we find, as the 
total coefficient of -r- , the expression 

or symbolically, Ly^\ and so on for the other coefficients^ 
The result then is 

It remains only after calculation of Ay^, Ay,, ...... Ay^, as 

functions of a?i, a:,, a:», to express these functions and X 

in terms of y^,y„ y,. 

[It appears from the manuscript that an example was to 
have been supplied here.] 
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[The nest Article may be considered supplementary to 
Chap. siv. Art. 10.] 

Singular Solutions of partial Differential Equations. 



^ idre's theory developed in Chap, xxill, for ordi- 
nan', may l)e applied also without essential change to partial, 
diflerential equations. Eegarding the independent variable 
s as receiving an infinitesimal change Zz through infinitesimal 
change, not in the values of the independent variables 



but in the values of the arbitraiy constants of the complete 
or in the forms of the arbitrary functions of the general mte* 
gral, and performing upon the given equation the operation 
denoted by S, we shall obtain a linear partial differential 
equation for determining the general value of Sz corresponding 
to any particular given value of z. If that linear equation be 
of a lower order than the differential equation given, then the 
equation expressing the value of z + 8z will be a limiting 
form of a solution less complete or less general than the com- 
plete or general solution of the differential equation given, 
and the given solution, formed by making the infinitesimal 
constants in the limiting form actually 0, will be singular. 

Conversely, to deduce singular solutions without the know- 
ledge of the complete or the general integral, we ought to 
construct the equations of condition for the reduction of the 
equation determining Sz to a lower order than the eqaatiou 
given, and the most general solution of the differential equa- 
tions of condition so formed, will be the most general expres- 
sion for the singular solutions of the differential 
given, 

Ex. (j)x — qy)' q + imx' (a — xp) = 0, 

in which 

^~ dx' " dy' 
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Kepresenting the first member of the equation by F^ we 
have, on operating by S, 

dp dx dq dy dz ~ ' 



and the conditions 



dF ^ dF ^ 



necessary to reduce the equation for 8z to a lower order give 

{px - qy) q — 2ma? = 0, 
{px-qy) ipx-dqy) = 0. 

From these we find 

p = 3m*ic*y*, q = m^ary'^f 

definite and simultaneous values of ^ and q^ which being sub- 
stituted in the given equation lead to 

z = 2m^x^y^f 

and this, as it gives the same values of p and q as those 
obtained before, will necessarily satisfy the given equation. 
It is therefore a solution, and firom the nature of the analysis, 
a singular one. 

Legendre shews that this singular solution is also dedu- 
cible firom the general integral of the given partial differen- 
tial equation. That integral is the result of the elimination 
of a fix>m the two equations 

{^ (a)}' — 2axip {a)'\-az^ mxy = 0, 

{if> (a) - ax] if>{a) - 2x<t> (a) + ;s = 0. 

To deduce the sing^olar solution he supposes ^ (a) to be not 
simply a function of a, but a function of a and of one or 
loth of the independent variables. He expresses the varia- 
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tion of ^(a) derived from this new source by S, and operating 
on the nrst equation with 8, finds 

{2^ (a) - 2ax} S^ (a) = ; 

therefore . ^ (a) = ax. 

Substituting this in the equations of the general integral, and 
eliminating a, we find 

z = 2m*a5*y* 
as before. 

Legendre states his theory of the derivation of the singular 
solutions of partial differential equations from the equations 
themselves with great brevity, but still as a general theory. 
And there is nothing in the statement that carries^ with it any 
apparent restriction upon either the order or the degree of the 
equations given. Until however we are in possession of a 
perfect theory of the genesis of partial differential equations 
we shall not be entitled to say that Legendre's theory of 
their singular solutions is a perfect one; for until then we 
cannot even define, in a perfectly general way, the nature, of 
the operation denoted by 5. 




( 13 ) 



[The nest three Chapters all relate to the subject of partial 
differential equations 01 the first order. The manuacripta do 
not appear to have received their final revision from Professor 
Boole. It is certain that he intended tlie contents of Chapter 

XXV. to form a part of the new edition; and it is highly 
probable, althougii not certain, that the contents of Chapter 

XXVI. and Chapter XXVII. were also to be included. 

The three Chaptera are mainly derived from two memoira 
Ijy Professor Boole, published in the Fhiloaopkical Traita- 



ft The 

fcliy Pr 

M-iietiona 



The firat memoir is entitled On Simultaneous Differential 
Equations of the First Order in which the Number of the 
Variables exceeds hy vwre than one the Number of the Equa- 
'Ions: it occupies pages 437. ..454 of the PkUoaophical Trans- 

Mons for 1862. 



The second memoir is entitled On the Difkrentiat Equa- 
ma of Dynamics. A sequel to a Paper on Simultaneous 

Differential Equuliims: it occupies pages 485. ..501 of the 

Philosophical Transactions for 1863. 

The first memoir was finished before Professor Boole had 

seen Jacobi's researches, which are cited at the beginning 

of Chapter xxvr; these reaearclies indeed could only just 

have been published. In his second memoir Professor Boole 

'ibes Jacobi's methods, refers to his own already pub- 

ihed, and points out the nature of the connexion betweea 
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CHAPTER XXV. 

ON SYSTEMS OP SIMULTANEOUS LINEAR PARTIAL DIFFEREN- 
TIAL EQUATIONS OF THE FIRST ORDER, AND ON ASSO- 
CIATED SYSTEMS OP ORDINARY DIFFERENTL4L EQUATIONS. 

1. The term simultaneous is here applied to a system of 
partial difTereDtial equations, to siguiiy that in that system 
tliere is but one dependent yariable, the general expression 
of which, as a function of the independent variables satisfy- 
ing all the equations at onee, is the object of search. All 
linear partial differential equations of the first order being re- 
ducible to the homogeneous form, we shall presuppose this 
reduction here. Under this form indeed the problem actually 
presents itself in Geometry, in the theory of partial differential 
equatioua of the second order, and in Theoretical Dynamics. 

We are sometimes led, in connexion with the same class 
of inquiries, to systems of ordinary differential equations 
marked by the peculiarity that the number of the variables 
exceeds by more than one the number of the equations. Such 
systems are intimately connected with the former — stand 
to them indeed in a similar relation to that which the 
Lagrangean auxiliary system bears to the sivffle partial dif- 
I'erential equation from which it arises. The theory' which 
explains this connexion, and grounds upon it the method of 
solution of both systems will form the subject of the present 
Ciiapter, 



Connexion of tite Systems. 

2. Prop. I. The solution of a system of simultaneous 
linear partial differential equations of the first order may be 



^ 
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made to depend upon that of a system of ordinary differential 
equations of tlie first order in which the number of the vari- 
ables exceeds by more than one the number of the equations. 

The system, of partial differential equations being reduced 
to the homogeneous form, Chap, sxiv. Art. 6, let n be the 

number of the equations, a:,, x^, a:,,, the independent 

variables, and P the dependent variable. 



Then from the n given equatioi 

dP dP 

dx'-- 



i determining 

dP 

' dx' 



B obtain an equivalent system of equations which, by trana- 



j position of its terms to o 



side, 



the redueed i'oriu 



dp 

' di„' 

" dx 



"'<i',„ 



..(1), 



Multiplying theee equations by the arbitrary constants 



ind adding the results, we have 
dP ,, dP 



+ (M„ + M„ + M.,)7 



+ (M„+V- +MJj 



■(2), 
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a single partial differential equation which, on account of the 

arbitrariness of Xi,^^, \,, is equivalent to the system 

from which it was formed. 

Of this equation the Lagrangean auxiliary system will be 



\ \ 



<^-4-t 



X^A^^ ■+ \A,^ ...... + \A^^ 



<^«+r 



M.r+M»-V + M. 



(3), 



whence, eliminating X, ,Xa, \j we have the system of 

ordinary differential equations 



dx^^ — A^^dxj^ — A^^dx^ ... .--A^jdx^ = 0" 
dx^a — -^igdbj — Aj^dx^ ...."A^^dx^ = 



(4). 



dx^ — ^ir^l — ^ar^a • • • •-^nr^n = 

These equations being included in the previous system (3), 
any integrals 

w = a, t? = J, t(? = c, &c« 

of them will be integrals of it. Therefore w, v, w,... will be 
values of P satisfying the partial differential equation (2). 
For they will be the only values which can satisfy it inde- 
pendently of \j, \, \. Hence they will satisfy the 

equivalent system (1), and the general integral of that system 
will be 

F{u, V, w, ...) =0 (5), 

the form of F being arbitrary. 

Thus the relation of the system (4) to the system (1) is the 
same as the relation of the auxiliary system of a single linear 
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partial differential equation to that equation. And the ground 
of this relation is seen to be the same in both cases. The 
one form necessitates the other. 

3. Instead of employing the above mode of deducing 
the- auxiliary system, we might employ the following which 
is practically more convenient. 

• 

Since any value P which satisfies the partial difierential 
equations determines P=c as an integral of the ordinary sys- 
tem, the latter must be consistent with dP—0 in its de- 
veloped form 

dP. dP . _^dP . ^ 



'a n^r 



T,,. . ^. dP dP dP 
Ehmmatmg ^, ^, ^_ 

by means of the n given equations (1), we have 

dP 
^- {dx^^- A^^dx^" A^dx^ - Ax^») 

dP 
-^ -IZ- ^dx^" AJ.x^- A^dx^ -Aa^n) 



dP 

Whence, equating to the respective coefficients of 

dP dP dP 

^n+l ^11+2 ^i^i^ 

we have the system (4). 

In the same way we can pass from the system of ordi- 
nary to that of partial differential equations. From the equa-- 
tion dP= 0, in tts developed form, we must eliminate a number. 
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of differentials dx,, dx^,... equal to that of tlie given equa- 
tions, and then equate to the coefficients of the remaining 



4 Lastly, the forma! connexion of tlie two systems should 
T>e noticed. The partial differential equations being given in the 
reduced form (1), the ordinary system may be constructed aa 

follows : For any differential coefficient, as -j — , in any 

column after the first, write the corresponding difFerential 
dx„^^, subtract from this the sum of dx^, dx^, dr.^, mul- 
tiplied respectively by the descending coefficients of that 
column, and equate the result to 0. The system of equations 
thus successively formed will be the auxiliary system sought. 



The transition from the ordinary to the partial system may 
be effected by the same rule, substituting only difierentials 
for differential coefficients. 

[It appears from the manuscript that an example was to 
have been supplied here.] 

Up to this point the theory of systems of partial differen- 
tial equations is in analogy with that of single equations, 
But here a difference arises. We do not know beforehand 
what number of integi'als a system of ordinary differential 
equations, in wliicli the number of variables exceeds by more 
tlian one the number of the equations, admits. 

The theorj' which removes this difficulty will be developed 
in the foUowmg sections. It will be shewn that a system of 
linear partial differential equations which admits of aolution 
by the assigning to the dependent variable a value which 
satisfies all the equations in common, must either itself satisfy 
a certain condition, or be capable of being developed into a 
new but equivalent system which will satisfy that condition. 
It will be shewn that when that condition is satisfied, the 
auxiliary system of ordinary, is capable of expression as a 
system of exact differential equations determinmg the inte- 
grals Bought. 
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It will be found convenient to express by a single symbol 
the aggregate of the operations to which the dependent vari- 
able 18 subject in the expression of a partial differential equa- 
tion. Thus the equation 

dz dz dz . 
may be expressed in the form 

A«=:0 

if we assume 

Under this convention the following proposition is to be 
understood. 

5. Prop. II. If AP = 0, A'P=0 represent any two 
homogeneous linear partial differential equations of the first 
order, then will 

(AA'-A'A)P=0 

also be a homogeneous linear partial differential equation of 
the first order, and it will be satisfied by all the common 
integrals of the equations from which it is derived. 

First, the equation will be linear. For, let a:, y represent 
any two variables whatever, or the same variable repeated, 
out of the set fiCj .... a?„, and let -4, B represent any functions 
of the variables a?^ .... x^. Then A may be represented by a 

series of terms of the form A -7- , and A' by a series of terms 

of the form B -7- . Hence (AA' — A' A) P can be expressed 
by a series of terms of the form 



dx\ dy) dy\ dxj ' 
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which, on effecting the differentiations, becomes 

dx dy dy dx^ 

tlie terms containing the second differential coefficients of P 
niutuallj destroying each other. Hence the equation 

(AA'-A'A)P=0 

will be a homogeneous linear partial differential equation of 
the first order. 

The constitution of the coefficients of this equation is easily 
determined. For suppose the given equations to be 

SO that 

then the equation 

(AA' - A'A) P= 

may be written in the form 

-■«-^.S -^-D- 

and, since terms involving second differential coefficients of 
P will disappear, this becomes 



+ (AP.-A'^J^=0. 



dx. 
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We see frura tliU that the form of the result ia the same as 
^ the A or A' from either equation operated only on the coeffi- 
deitts in the other equation. 

Secondly, the above equation will be satisfied by all the 
common integrals of the ei^uations from which it ia derired. 

For, let ^ = c be a common integral of 
AP^O and AT=0, 
then 

A^ = 0, A'^ = 0. 

Performing 'on these the teapective operationa A' and A, 
operations which involve only differentiation together with 
ugebraic processes, we have 

A'A0 = O, AA'^=0, 
I whence, by subtraction, 

AA'0 - A'A^ = 0, 

(AA' - A'A) = 0, 

rom which it appears that is also an integral of the 

(AA'-A'A)P=0, 
B was to be shewn. 



Prop. III. If by the above processes of reduction and 
Brivation we convert a system of partial differential equa- 
tona into a new system, such that if expressed in the form 

AjP= 0, A^= 0, A„,P= 0, 



tne condition 



(AA-AA)-P=o 



pAiall for each pair of equations be identically satisfied, then 
" em of ordinary differential equations corresponding 

new system will admit of reduction to the form of 
act differential equations, the integration of which will 
taablo us to construct the general value of P satisfying the 
yatera given. 
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Ist. Suppose the given aystem of n equations reduced to 
the form (1), marked by the peculiarity that n of the differen- 
tial coeificienta appear only in successive equations and with 
the coefficient unity. Then talcing any two of those equar 
tions (we select the first two), we have 

d , d . d 
A,= T-+^„-^ +A„-s , 

from the forma of which we see that the derived equatu 

{AA-A.A,)P = 
cannot contain either 

dP^ dp 

dx^ dXa' 

It can only, as appears from Art. 5, contain the diSet 
coefficients ~ 

dP dP 

dx^^' dx,^' 

and must be of the form 



'dx^, 'dx,^ ""dx,^ 

It cannot therefore be an algebraic consequence of any of the 
equations of the system (1) from which it was derived. It is, 
unless by the vanishing of 5,, ...,B, it present itself as an 
identity, a new equation algebraically independent. Com- 
bining this with the former ones, we have a system of « + l 
equations admitting of the same reduction as to form fol- 
lowed by the same subsequent process of derivation. And 
the result of each of these completed steps is to convert the 
system into one containing one equation more than before; 
but containing in each of its equations one term fewer than 
before. The process must then end either in the genesis of a 
system of partial differential equations such that the further 
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application of the process of derivation of Prop. II. Bhall only 
lead to identities, or in the emerging of the Bystem 

^-0, ^-0 ;^.0. 

The latter Buppoaition would imply that P is a constant. 
The consequeacea of the former we proceed to examine. 



The final system of linear partial differential equations 
will be of the same type (1) as the original system, but will 
differ from that system in that n will be increased, and r 
diminished by the same amount. We shall therefore simply 
state the form (l), only under the condition 

(A.A,-AA)-P=0, 
and with the altered values of n and r. 

First, then, the common inte^als of the new system will 
be the same as those of the origmal system. Thia is evident 
from Prop. II. 

Secondly. If we write 

m = n + r, 
the first equation of the system (1) will be 



and the auxiliary Lagrangean system of this will have j» — 1 
independent integrals 



among which the k— 1 known integrals (Chap, XXIT. Art. 5) 

x^ = c^, a;, = Cj, x^ = c^ 

G— 2 



L 
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are included. And the general value of P satisfying the 
atove first equation will be 

-p-i't",," «„)• 

The assumption P= a;, would not satisfy the eaid equation, 
for it would lead, on substitution, to 1 = 0. Hence we infer 

that while the functions m,, w,, u are independent with 

respect to each other, they are also independent witli respect 
toic,, so that the wi functions «,, w^, m^,, a;,, are mu- 
tually independent in the sense explained in Chap, xxiv. 

Let u3 now transform the equations of the system (1) after 

the first by introducing m,, w,, W™-,) ic, as independent 

variables. Those equations being 

A.P=0, A,P = 0, 

the result of the transformation will be (Chap. xxiv. Art 7) 



{i.«,)£ + (A.-j£ + (A.O: 



But P = x, being an integral of each of the equations of 
the system (1) except the hret, as appears from their forms, 
we have 

A^, = 0, A„aT, = 0, 

thus the last terms in the transformed system vanish. Further, 
the coefficients of the remaining terms reduce to functions rf 

«!, M w^, mereiy. Tor, considering the coefficient A,u„ 

we have 



{AA-AA)wi = 
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which, since A,u, = reduces to 

A,A,«, = 0. 

Hence A,u must be & solution of A,P= 0, and therefore a 

function of «,, «,, m^,. And so for the others. It 

results therefore that the transformed system is 



,. .dp ,, . dP ,, . dP ^ 



.. B,^j heing the actual independent variables of the 



system. 

But the transformation having involved no loss of gene- 
rality, for a new system of vi independent variables was 
simply substituted for an old one, the condition 

(A.A,-AA)-P=0, 

aatisfied before, will continue to be satisfied in the new sys- 
tem represented symbolically in the form 

A,P=0, A,P=0 A.P=0. 

Any common integrals of this system will also be common 
integrals of the previous system. For as functions of 



I they will satisfy the first equation of that system, and they J 

Irill satisfy the other equations, because the present system is I 

but a transformation of those. The converse is equally mani- 1 

Thus a system of n partial differential equations contain- 
bj^^mdepm^i^TO^Ue^^^^iflfrins^^^tw^cOT^^^^^^J 
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tioii, Iiaa in Tirtue of that condition been converted into a 
system of n — 1 equations between m — 1 independent vari- 
ables, and satisfying the same condition. This then ia con- 
vertible into a similarly constituted system of « — 2 equations 
containing m — 2 independent variables, and so on till we 
arrive at a final single partial differential equation containing 
m—n + 1 independent variables. This equation has m—ti, 
that ia, r integrals, and these are the common integrals of the 
system (1). 

But the system of ordinary differential equations corre- 
sponding to (I) is in number r, and is satisfied by all the 
common integrals of that system. Hence these differential 
equations must admit of reduction to the exact form. 

7. We may deduce from the above investigation the fol- 
lowing Rule. 

To integrate a aystera of simultaneous linear partial diffe- 
rential equations of the first order, 

EuLE. Reduce the equations to the homogeneous fbmi 
(l), express the result symbolically by 

A.p = 0, A,P = 0, A„P = 0, 

and examine whether -the condition 

is identically satisfied for every pair of equations of the sy»- 
tera. If it be so, the equations of the auxiliary system, 
Prop. I., will be reducible to the exact form, and their inte- 
grals being 

u = a, v = h, w = c, 

the complete value of P will be F{u, v, w, ...), the form off 
being arbitrary. 

If the condition be not identically satisfied, its application 
will give rise to one or more new partial differential equa- 
tions. Combine any one of these vrith the previous reduced 



I 
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syateiD, and again reduce in the same 'way. With the new 
reduced syatem proceed as before, and continue tliis method 
of reduction and derivation until either a system of partial 
differential equations arises between every two of which the 
'above condition is identically satisfied, or, which is the only 
poaaible alternative, the syatem 



f-O, f = 0,. 

dx ax. 



appears. In the former case tlie syatem of ordinary equations 
corresponding to the final system of partial difi'erential equa- 
tions will admit of reduction to the exact form, and the gene- 
ral value of F will emerge from their integrals as above. la 
the latter case the given system can only be satisfied by sup- 
posing P a constant. 

Ultimately then the determination of P depends on the 
solution of a syatem of ordinary differential equations reduci- 
ble to the exact form. Thia does not mean that each equation 
of the system ia reducible to the exact form, but that the 
equations may be combined together so aa to form an equal 
number of equivalent equations of the exact form. Generally 
when we know this combination to be possible it is easy 
to effect it, and best to endenvonr to do so. We might how- 
ever employ the method of the variation of parameters as fol- 
lows. Supposing _p the number of differential equationa make 
all but ^ -i- 1 of the variables constant, integrate the reduced 
Bystem, and then seek to satisfy the unreduced syatem by the 
same seriea of integrals with tlie arbitrary constants as new 
variables. The succeaaive integrations and transformations 
of this method would amount to the same thing as those 
upon which the second part of the demonstration ot Prop. iji. 
Teats*. 

Lastly, given a system of ordinary differential equations 
containing a superfluous number of variables without know- 
ing how many integrals they admit, we must, supposing 
P=cto be any integral, construct the corresponding system 

* It WM thai bdeed that tbe >athor wu fint led to that thfloiy. 



1 



J 
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of tomogeneous partial differential eqtuitions satisfied by P, 
and apply to them the foregoing Rule. 

8. Ex. Eeqcired the integrala of tlie simnltaneouB par- 
tial differential equations 



Bepresenting these in the fonn A,P=0, A^ = 0, it will 
be found that the equation 

(A.A,-A,A,)J'=0 

becomes, after rejecting an algebraic factor, 

dP dP ^ 

''^■^dr^' 

and the tbree equations prepared in the manner explained in 
the Bule will be found to be 



No other equations are derivable from these. We conclade 
that there is but one final integral. 

To obtain it, eliminate 

dF dP d^ 
dx' di/ ' dt 
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from the above Bjstem combined mtb 

dP, dP, dP, ^dP. - 
-dx^^dy + ^dz^-j^dt^Q, 

and equate to the coefficient of ^ in the result. We 

find 

& — (< + 3aj^ dx'^ydy — xdt = 0, 

the integral of which is 

« — ip^ — a?'— ^ = c. 

An arbitrary function of the first member of this equa- 
tion is the general value of P. 

[It appears fix)m the manuscript that another example was 
to have been added here.] 
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CHAPTER 5XVI. 



[CH. XXVI. 



HOMOGENEOUS SYSTEMS OP LINEAR PAETUL DIFFEEESTUL 
EQUATIONS. 

1, The theory of homogeneona syatema of linear partial 
differential equations in which when expressed in the synr- 
hoUc form 



A,P=0, A,P = 0, A„P = o. 



..(1), 



the condition 



(AAi-^Ail'-O.. 



..(3) 



ia for all combinations represented by i and J satisfied in 
■virtue of the constitution of the aymbola A;, Aj, forms the 
subject of important researches by Jacobi {Nova Methodua... 
Crelle's Journal, Vol. Lx. p. 1). The following are the most 
important of his results. 

Ist. An integral of any one equation of the system being 
found, other integrals of the same system may be obtained 
without integration, by a process of derivation founded upon 
the condition (2). 

Let be an integral of the first equation of the system. 
Then is the equation 

A,0 = O 
identically satisfied. 

Also the condition (2) being satisfied in virtue of the con- 
stitution of the symbols, we have 

(Aiaj-A,Aj)<^ = 0; 
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and in particular, making i = l, and separating the terms, 

wtich reduces by a prior equation to 
A,Aj^ = 0. 

It appears from tliis that A,i^, if it do not rcdnce to a con- 
stant, la an integral of the tirst equation A,i^ = 0, and, if it 
prove to be not a mere function of tp, a new integral. 

This process may be repeated upon the new integral with 
a similar alternation of results. It will he evident from this 
that if we confine our attention to the two equations 

A,P=0, A,P=0, 

and suppose, as before, tfi to be an integral of the first, then 
will 

A,0, A,(A,.^), A,{A.(A,.^)],... 

or, as these may be expressed, 

Aj^, A)Y. ^iVi ■■• 

■'be also integrals of the first equation ; and this process of 
derivation may he continued until we arrive at an integral 
A/"^ which is not independent, but is expressible as a func- 
tion of prior integrals 



and, sooner or later, such a result must present itself, since 
the number of independent integrals is finite. 

It is further seen tliat the most general symbolic form of an 
integral derivable from the root integral ^ is 

A,-A,B A,"^, 

, o, ff, ft being positive integers. 

The above remarkable theorem was in some degree antici- 
i pated by the researches of Poisson. 
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2nd!y. Jacobi shews how by the aid of such derived in- 
tegrals of the first equation of the system a common integral 
of the first and second equation may be found, and how from 
this integral and its derived series a common integral of the 
first three equations of the system may be found, and bo on, 
until a common integral of the entire system haa been as it 
were huilt up out of previons integrals of less general appli- 
cation. 

Let A, ^', ^", 0'"-^) represent a series of independent 

integrals of the equation A,P = 0, of wliich ^ is the root in- 
tegral, and the rest are derived from it by successive applica- 
tions of the operation denoted by A,, so that 

0' = A^, ^f"-i' = A/-'^; 

also let Aj*^^ he not a new integral but a function of 

^, ^', ^t^-'i. 

Now ^, (^', (^(«-i* being particular integrals of AjP=0, 

the function ii'(^, ^', ^if-'i) will also be an integral of 

the same equation irrespectively of its form. Let us inqoire 
wliether the form of the function can be so determined, as 
to render it also an integral of the second equation A^P = 0. 

We have then to satisfy the equation 

A,i?'{<^, 0', <^<^-'l)=0. 

By the principles of the Differential Calcalus this equatioa 
assumes the form 

But A,i^ = ^', A,i^' = ^", A,^'"-** = ^f^-« ; 

lastly, A,^'*'-" may bv liypothesis be expressed in the form 
/C^, 9^' ip>'^~"). Thus the equation to be satisfied is 
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„dF, ,..iF . ,„„„ dF 



+/».*■ ♦"-'')J^^-». 

a linear partial differential equation of which the auxiliary 
system Is 

Jitll m 



/».< 



Now the integration of this system may be made to depend 
upon that of an ordinary differential equation of the (ji — 1)"' 
degree between the two variables <^'''~^' and if>. 



d<l> ^' ' dtf> 0' * 

d4><^-» _ /{(!>, i>', 0"-^') 

Differentiating the laat equation with tespect to if>, and attend- 
ing to the former ones, we shall be able to express ^ ia 

terms of the variables 0, 0', ^l"-''. Proceeding with 

thia in the same way and continuing the process we ahall be 
able to express the series of differential coefficients 

in terms of 0, <^', ^*'~^'. From these ^ — 1 equations, 

eliminating i^', 0", 0^"-^, we shall have a final equation 

between 



that is, a differential equation of the Iji — 1)"" order between ^ 
and ^1^1. 



J 
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The complete integral of this equation will be of the form 

<^''-!'=/(f <:..«. c^-i)- 

Differentiating this fi.-2 times in succession with respect to 
A, and continually substituting for the differential coefficients 
of (^(^-'i their values aa before assigned in terms of 

1^, tf,\ ^("-11, 

we shall have a system of ^ — 1 equations connecting the 

above variables with the constants c,, c^, Ca-i- Finally, 

solving these equations with respect to the constants, we shall 
possess the integrals required in the form 

F^{<P,<J3, ^'"-") = c.. 

F^.i{<f>,'P' •^'-"} = c^-i, 

and each of these will be a common integral of the first two 
equations of the given system (1). 

[On the back of a page of the manuscript the following 
paragraph occurs, which seems to have been intended as a 
simplification of the preceding argument which begins with 
"The complete integral."] 

Suppose that a first integral of the equation can be found. 
Its form will be 

^[^'^"^ '-df- -d$^)='- 

Substitute in this for the differential coefficients of (^"") 
their values before assigned in terms of <^, (j>', 0",.-.0'"~", and 
we have an integral of the system (3), and therefore a com- 
mon integral of the first two equations of the system (l). 

[We now return to the place at wliich we inserted a para- 
graph.] 

Just in the same way Jacob! deduces a common integral of 
■ the first three equations of the system (1). For representing 



\ 
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any one of the first members of the above Hjatem by if-, and 
deriving thence the new independent integritls A^-\jr, A,'i^,.,. 
he subatitutea an aibitrary function of tuese for P in the 
equation 

A,P=0. 

It is evident that the solution of the partial differential 
equation so found will again be reducible to that of an 
ordinary differential equation between two variables. And 
80 the process ia carried on till all the equations are satia- 
ted. 

2. The above remarkable process was developed by Jacobi 
in connexion with the theory of non-linear partial differential 

auations of the first order. In that particular connexion it 
mits of certain reductions tending to diminish the order of 
the differential equations to be integrated. But these do not 
affect the general principle of the method. It was in this 
special form that the theory of the solution of simultaneous 
linear partial differential equations originated. Jacobi does 
not consider the theory of equations in which the condition 
(2) is not satisfied ; but the language in which he refers to 
the condition shews that he had speculated upon the general 
problem — and it is difficult to conceive that he should have 
meditated upon it and not arrived at its complete solution. 

[The manuscript here gives the first two words of 
the passage from Jacobi'a memoir which is quoted in the 
Philosophical Transactions for 18S3, page 486.] 



J 




[CH. XXVII. 



CHAPTER XSVn. 



OP KOy-LIHEAR PARTIAL DIFFERENTIAL EQUATIONS OF THB 
FIR8T ORDER. 

t. In treating the present subject we shall first consider 
that class of non-linear partial differential equations of the 
first order which involves two independent variablea, and 
then proceed to the general theory. The reason fot this 
proceaure is that the particular theory, though of course in- 
cluded in the general one, rests upon a somewhat simpler 
basis, and it was in fact developed by the labours of 
Lagrange and Charpit long before the general theory was 
known. The latter we owe to the independent researches 
of Cauchy and Jacobi, 

[Here the mannsci'ipt refers to the matter contiuned in 
Chap. siv. Arts. 7 to 12 incluBive; and thea^passes on to 
the general theory.] 



General Theory. 
2. Given an equation of the form 

z = ^ (ic^, a;,, . . . a:„ a^, Cj, 



.«.), 



the number of arbitrary constants a,, a„ ... a^ involved being 
equal to the number of the independent variables x^t a^ •••'^ 
we obtain by differentiation and elimination of the constant 
a partial differential equation of the first order. Of this the 
proposed equation is said tc constitute a complete primitiv«>:i 
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The form of tlie above process which it seems best, as 
throwing light upon the inverse problem of deducing the 
•complete primitive from the partial differential equation, to 
employ, is the following. Let the given primitive, solved 
with respect to one of the arbitrary constants a^, be presented 
in the form 

Differentiating with respect to each of the independent vari- 
ables we have a system of n equations of the forms 

^(a?j, ... a?„, z,jp^, a^, ... a„)=0'j 

[ '(2). 

f, (iTj, ... x^, «,^„, a,, ... an) =0 ^ 

These n equations enable us first to eliminate the n — 1 

constants a^, a^, and so deduce the partial differential 

equation sought in the form 

F^{x^, ... x^, z, p^, ... jPn) =0 (3); 

secondly to determine the w — 1 constants as functions of 
iCj, ... x^j ^>i^i> »" Pn ^^ the forms 



-^2 (*^1' ••• ^«' ^^ Pv ••• Pn) ~~ ^S 



•^n ("^i > • • • ^n > ^i Piy • • • Pn) — ^n ^ 



(4). 



As the system formed of these n — 1 equations, together 
with the previous one, is merely another form of the system 
(2) obtained by directly differentiating the primitive, it follows 
that if from these equations we deduce the values of^j, ... Pn 
as functions of 0?^, ... a?„, a^, ... a^, and substitute them in the 
equation 

dz=pjdx^+p^dx^+ ...+Pndxn (5), 

they will render that equation integrable, and its integral 
will be the complete primitive (1), the constant a^ being re- 
gained by integration. 

B.D.E. II. 7 
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Exainining tlie eyalem (3), (4) we see that the first mem- 
Lera of all the equations which it contains are fimctions of 
a;,, ... x„, z, p^, ... p„, while the second membera are con- 
stants. The question then ariscB, What mutual connexion 
exists among these lunctions in virtue of which they yield 
values of p,, .-•jo,, which render the equation (5) inte- 
grable? 

The answer to this question must involve the entire theory 
of the solution of partial ditJerential equations of the first 
order, so far as relates to the determination of a complete 
primitive. Given a partial differential equation of the form 
(3) it is evident that if we can construct a system of associated 
equations (4) possessing the character above described, the 
final value of 3 obtained by integration of (5) will botli 
satisfy the given equation and contain the requisite number 
of arbitrary constants. It does not follow tirom this that 
it will be the only complete primitivej but it will he a 
complete primitive. 



3. The 

Proposition 


relation 


sought 


is 


expressed in 


the 


followin 


Pkoposition. If 














•f(»., 


•■ a^ni e 


?. 


...yj-a, 








■SK, 


.. ic., z 


Pi 


-i-J-S 







represent any two out of a ayalem of n mdependent equaticm 
guch that the values ofp^, ...p^, thence d-etennined would make 
the equation 

dz =p^dx^ + p,dx^ + ... -i-p^dx^ 

integrahh, then the first members of these equation* being 
represented Jhr simplicity by F and C>, the condition 

^ UdF , dF\ d^ dF/d^ d^W „ 
^Kdx^P^d-zjd^-dpX-d^.+P^-dzJr''' 

the summation extending to all values of i, from \ lo n inclatiw, 
will he satisfied identically. 
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RedprocaUyy if the above condition he satisfied identically 
fir each binary combination of junctions in the proposed system 
of eqiuUions, and if these fiinctions be independent^ then the 
values ofp^ ,. . . jt?^, ew functions ofx^^ . . . a?,» , «, which they yield, 
wiU make the equation 



dz =Pi^i +i?2^a + ••• +P*^i. 



integrable* 



It will be convenient to begin with the particular case in 
which the proposed equations do not explicitly contain z^ the 
particular pair to be considered being represented by 

F{x^, ...x^,p^, ..:p^) = a, 



Differentiating with respect to Xi, and regarding ^j, ,.^p^ as 
functions of the independent variables, we have 



dF dF^ _,dFdpn 

dxi dp^ dxi *" dp^ dxi 

d^ d^ dp^, 4.^ ^ 

dXi dp^ dXi '** c^n dXi 



=0, 



= o,J 



(6), 



to which we may give the form 

^ = - 1 — ^n 

dXi ^ dp J dXi 

d^ ^ d^ dpj 

dxi " ^ dpi dxi^ 



(7), 



the summation with respect toj" extending from y= 1 io j-n 
inclusive. 

d<b 
From the first of equations (7) multiplied by -j- subtract 

7—2 
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the second multiplied by -r- , and sum the result with respect 
to i from ^ = 1 to t = 71 inclusive. We have 



t 



dF d^ dF d^\ 



\dXi dpi 






dpi dxj 



- _ :S ^ (^ ^ dpf^dFd^ ^\ , . 

* ^\dpj dpi dxi dpi dpj dxJ 

The expression under the double sign of summation in the 
second member vanishes when ^=y; we may therefore re- 
strict the summation to unequal values of i and j. Now 
as for any particular combination of values, e.g. 2, 3, there 
would exist in the completed member both the terms cor- 
responding to i = 2,^* = 3, and those corresponding to ?* = 2, 
i= 3, it is evident that if we employ the symbol 2y to aenote 
summation with respect to different combinations of { and /, 
the second member of the last equation may be expressed in 
the form 

^ fdF d^ dpjdF d^ dp^ 
^^\dpj dpi dXi dpi dpj dx^ 



dFd^^^dFd^dpA 
dpi dpj dxj dpj dpi dxJ ' 

>\fdpi_^] 
J \dxj dxJ) * 



or 2 (fdFd^^dFd^ 
'^ \\dpi dp J dp J dp, 

so that the equation (8) becomes 
^dF d^ dF d^\ 



^' \dXi 



dpi dpi dxJ 



^[\jdpi dpi dPi dpJ\(iXj dxJ) ^ '' 

The number of terms of which the second member ex- 
presses the sum is thus — ^- — - , and it will be observed fhat 
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as to any particular term it makes no difference in what order 
the numerical values of i and j are assigned to these quan- 
tities; e.g. whether for the combination 2, 3 we make i= 2, 
j = 3, or i=3,j=2; but we must confine ourselves to one 
order. 

1 JTow when the equation 

H dz=p,dx^+p^dx,+ ...+p„dx„ 

V integrable in the manner here supposed, we have for all 
combinations of I'andj", 

dpi dpj 

dxj dx^ ' 



I therefore vanish. 



All the terms in the second member of {i 
and we have 

^ /dF^ d^_dF ^■\ 
* \dxi dpi dpi dxJ 

This is the direct form of the Proposition under the parti- 
cular limitation supposed. 

As F, ^ represent, under the same limitation, any two of 
the first members of the n equations (3), (4), which determine 

pi, ...p., there wilt exist — - equations like the above. 

It is usual to employ for brevity the notation 

^ ^ fdF df dF d^-^ 



"* \dxi dpt dpi dxi/ 



'IF<S,], 



sod this being done the above system of equations expresses 
— functions of the form {F^F^ a 



the- 



sUn 



functions of the - — - quantities of the form ~i;~^'- 

It is hence that the vanishing of the latter series of quantities 
secures the vanishing of the former. 



The converse truth will thi 
— --■ quantiti 



fore be established by shewing 
of the form -f^ — -P- are, when 
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F^, -Fj, ... i^„ are independent with respect to J^u p^j •••!>•> 

fi (t^ ~- ll 
expressible as linear homogeneous functions of the — ^-^ — - 

functions [FiFj]. 

To avoid complexity of expression I shall establish this for 
the particular case of n = 3, and shall shew that the reasoning 
is general. 

The functions F^, F^, F^, being independent with respect 
*^ PvP^y i^8> ^^ determinant 



dF, 


dF, 


dF, 


dp,' 


dp,' 


dp. 


dF, 


dF, 


dF, 


dp,' 


dp.' 


dp. 


dF, 


dF, 


dF, 



dpi ' dp, ' dp, 
does not vanish. This determinant we shall denote by A. 

In (9) writing for J?^and ^ first F^ and jPg, secondly F^ and 
jFj, thirdly F^ and F^^ we have on changing signs the system 



^"^ \dp^ dp^ dp^ dpj\dx^ dxj 
L » ^ \dp^ dp^ dp^ dpj\dx^ dxJ 
L^^^^ " [dp, dp, dp, dpj [dx, dxJ ■*" - J 



V (10). 



dF: 



dF^ 



Multiply the first equation by -i— * , the second by -^ , the 



third by -y-^ and add. Then 
dPi 



-f [^-^J - f ™ - f K^J - ^ (t-^. 
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wlience as A doea not vaniali we liave, on dividing by it, the 

- j-^ expressed as a linear homogeneous function 

of [F,F^, [F,F,f and [F,F,]. 

T n 1.- 1 ■ t ,-1. ^P, dF, dF, 

In like manner multiplying the equations by -J—, -^, -j— * 

respectively, and dividing by A, we obtain -P — ,°' as a 

similar linear liomoeeneous function, and lastly, multiplyine 

. dF, dF, dF„ 

t>y -j—^, -J—, _j~ I ^nd proceeding as before, we obtain 

3— — J-* BS a siinilai' linear bomofrcneous function. 

From all whicli it follows that when [F^F^, K^J, [F^F^ 
vanisb, then 

dp, _ dp^ dp^ dp^ dp^ dp^ 
dx dx, ' dx dx ' dx dx, ' 



The reasoning ia general in its nature. If ^, , F^, ...F^ 
are independent witli regard ^o p^,p^, ...p„, tlie deteniuaaiu 



*,■■■ 


•■■*. 


dF. 
dp,- 


dK 



does not vanish. This determinant is from its constitution 
as a determinant linear and homogeneous, not only with 
respect to any row or column of elements, but also with 
respect to the possible binary combinations which can be 
formed of two rows or columns, ternary out of three rows or 
columns, &c. provided that these combinations are themselves 
of the form of determinanla. In the language of the theorj^ 
auch combinations are called minor determinants. Hence if 
we construct the system of equations represented by (10), and 
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obaerye that the coefficients of any particular terra of the form 

Jli — Zi in the several equations form a system of such 

minora to the general determinant (11), it will be plain that 
the equations can by multiplication and addition be brought 
to a form in whicli the coefficient of that particular term will 
l>e A. At the same time the coetEcienta of all the other 

terms of the form -^ — i-^ will vanish. For a little atten- 

«3-j axi 

tion will shew that they will be what the determinant A 
would become on making two of ita columns or rows of ele- 
ments equal, and therefore will be identically equal to 0. 

Thus the Proposition is generally established for the case 
in which e does not explicitly appear in the functions 

F^,F^ F^. 

When s does appear in those functions the equations [6) 
will be replaced by 

dXi "' lis dp^ dxi dp, dxi ' 



:+Pi 



d^ d<^ dp, (?* dp„ 



dj^i "' dz dj), dxi dp, djCi 

from which it ia seen that the theorem above established will 
only need to be changed into the form employed in the state- 
ment of the general Proposition. 

As the above ia one of the most important propositions in 
the entire theory of Differential Equations, it may be deaire- 
able to illustrate it by examplea. 

[There are no examples in the manuscript.] 

4. We resume the general theory. 

The integration of non-linear partial differential equations 
may be effected by two distinct methods, both resting npou 
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Habt. 

^nbe ground of the above Proposition. The first of these I 
Bimethods, originally eatabliahed bj a different analysia froml 
"that which will here be employed, was discovered by CauchytJ 

(Keercices d'Analyse), and rediscovered by Jacobi (Creifea;! 

Journal). The second method, discovered by Jacobi at a | 

later period, forma the subject of his posthumous memairoj 

Mova Metkodus.. 



Cauchy's Method, 

We will, as before, begin with the case ia which z does J 
not appear explicitly in the proposed partial differential equa-«j 
tion, which we shall represent in the form 1 

F,{x^, x„ p„ iJ„) = (1). 

"We have seen that to find a complete primitive of the I 
equation it is necessary and sufficient to construct a series of 1 
equations 

■^>K. ^n. P. P.)=<:] 

(^)- 

i^,(a!, ar„ p„----.p„) =aj 

ich that not only shall the conditions 

[F.FJ = 0, [F,F,] = (3), 

mnecting the new functions F^, F„ with F^, be identl- 

Jly satisfied, but also the series of conditions 

K^.l = o W. 

f^ and F^ representing any two of the new functions re-' I 
id to. 

The first of the above series of conditions amounts to 1 

lis, that i^j i^ must be integrals of the partial differen- i 

tl equation 

K^-o (5)- 
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It is the peculiar aim of Cauchy'a metliod to determine the 
integrals so as to cause the second series of conditions to be 
satisfied also. And it is shen'u that this will be attained if 
the integi-ala of (5), which foim the first members of (2), are 

auch that the particular values which », p^^ assume when 

a;» is made to receive any constant value, as 0, are differential 

coefficients with respect to a:;,, a; of any single function 

of those variables, the form of which may be arbitrarily 



1 



The necessity of this condition is obvious. If the general 

values of ^,, p^ are differential coefficients of a function b 

with respect to a;,, a:,, then the particulai- forma which 

p , p assume when x^ receives any constant value are 

simply differential coefficients with respect to a; a;^ of 

what s becomes under the same circumatances. To prove its 
sufficiency we must shew that when it is satisfied the condi- 
tions represented by {4} will be satisfied also. 

Since F^ and F^ are integrals of {F^F\ = 0, 

[F,F„] = 0, [F^F^ = Q (G). 



Also, since if in (l) and (2) we give to x, a particular con- 
stant value, as 0, and then in (2) regard^, as a function of 

^v a^M. i'l- P-i 

determined by (1), the system (2) will virtually contain only 

ai,, a;^,, p^, Pn-i^ 

of which p^ p^^ are differential coefficients of a single 

function with respect to a:,, a;._,, it follows from the pro- 
position of Art. 3, that any two functions F^ and J^ will 
satisfy mutually the condition 



s; 



=„_, fdF, aF\ _ dF„ ^A ^ ^ 
"^ ydxt dpi dpt dxj ' 



the differentiations having reference to 

Xi, a;^i, p,, p^^. 



ART. 4.] EQUATIONS OP THE FIEST ORDER. 107 

explicitly as they appear in F^ and i^j, and implicitly as 
involved in^,^. Thus the developed form of the above equa- 
tion is 

^'^ \\dXi dp^dxj\dj>i dp^dpj 

\dpi dp^ dp J VdaJ, dp^ dxj) ' 
the forms of -^ and -^ being determined from (1). 

Perfonning the multiplications, the above equations wiU be' 
reduced to the form 

'^ \dx^ dpi dpi dxJ 

^ dK 2<=«-i fdpn dF,_dp^ dF\ 
dp^ **^ \dxi dpi dpi dxJ 

_ dF, 2<=n-i (dpn dF, _ ^« ^a") ^ .... (1). 

dp^ ^^^ \dxi dpi dp^ dxi) "**^ '* 

But from the form of the total differential of (1) we see that 

dFy dl\ 

dPn^_dXi dp^ _ dp^ 



dx^ dP^' dp^ dF\' 
dpn dp^ 

Hence 



2i=«-i (dpn ^a _ dp^ dF\ 
*"^ \dxi dpi dpi dxJ 



\dpj *'^ \dXi dpi dpi da\J 
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Now since by Art. 3 

*~^\cirj dpi dpi dxj ' 
we have 



^^'^ \dxi dpi dpi dxJ 



■" \dx^ dp^ dp^ dxJ' 



therefore 



^*=^ \dXi dpi dpi dxJ 



[dp J \dx^ dp^ dp^ dxJ 



In the same way 



*"^ \dxi dpi dpi dxJ 



\dpj \dx^ dp^ dp^ dxj ' 



The substitution of these values in (7) gives 



^*^ \dxi dpi dpi dxJ 



\dpj \dp^ \dx^ dpn dpn dxJ 
dPn \d^n dp„ dp^ dxJ) ' 



or 
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^'"^ \dxi dpi dpi dxj 

d^n dpn dpn dx^ ' 

<i^n(^ Hit ^^a ^^ 
*=^ \fl^a;< dpi dpi dxj ' 

which is precisely the equation 

[FaF,] = 0. 

We see therefore that to solve the partial diflferential 
equation 

F^{X, X^, p^y i?n)=0, 

it is only necessary to construct the linear partial differential 
equation 



^'^\dxi dpi dpi dxJ ' 



and to obtain n — 1 independent integrals of this 

F^ix^y a^n, Pi,....p„)=«2» 

Fn{x^j X^, P,, P^=^a^y 

such that if we determine from these conjoined with the given 

equation the values of p,,....^^, then those of j!?^, p^. 

shall, when x^ is made constant, be the partial differential 

coefficients of one and the same function ot aj^, x^^ with 

respect to these variables in succession. 

Now provided that we can find all the integrals of the 
above partial differential equation the particular determination 
required may be effected in the following manner. 

The Lagrangean auxiliary system consists of 2w — 1 ordi- 
nary differential equations 

dx^ dx^ dp, dp^ ,. 

~7f^^ ~r^~^" "^ 

^ ^1 dPn d^i doCn 
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These admit of 3n — 1 integrals, one of which will be i^, = c, ; 
and this will agree with the given equation if we make c, = 0. 
We have therefore, besides the particular integral J^ = 0, 
2/i — 2 integrals of the form 

^, (a;,, ...iF„, ^j, ,..p,) = Cj 



0s»-i (a:. , ■ ■ ■ if, , ^. , . . .^ J = c,^ 

Now suppose it required to find a value of s as a function 
of ar,, ...Xn, which shall satisfy the given partial differential 
equation, and shall reduce when a;, = (or any numerical con- 
stant) to a particular given function of x^, ...a:^,, which we 
will represent by 1|r{x^,,..x^^), Then on the assumption 
that a:, = 0, we nave first the given equation 

z = -^(a;,,...a:^,} {10)j_ 

secondly the derived equations 

„ _ tl^(^- T,.,) , 



de. 



10)^ 



(II). 



Pn-l-- 



J^(t„...x,J 



Make in the 2n — 1 integrals a;, = 0, and suppose at the 
same time x„ ... a:^,, j),, ..■^„_, to receive therem the same 
values as in the above derived equations. Then from the 
3n — 2 particular equations which we thus possess in the two 
systems united (particular because under the assumption that 
"'n^'*)! ^s can eliminate the 2n — 1 particular values of 
a;,, ...x^j, p ...p„_^, and so obtain n — 1 equations among the 
constants. These express the conditions which are necessary 
andsufficient in order tliat the values of j;,, ■■■p^.i thnsderivM 
from the integral equations may, when a:„ = 0, agree with the 
values assigned in (11). Accordingly if we substitute in these 
equations of condition for o,, . . .c^_, the general values 0, , -- .^n,_i 
we shall obtain n — 1 equations between a;,, ... x„, p,, ...p„, 
which will at once be particular integrals of the system (8), 
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and will possess the property that the values of p|,...p,^ which 

they in conjunction with F^ = give will when x„ = reduce 

to the values given in (11). Hence these values with that of 

derived from the same equation will make 

t^s— P, t^, — —p„dx„ = 

an exact differential equation. In the integral of this it will I 
only remain to determine tlie constant so as to make the value 1 
of s agree with that given in (10). Alt the conditions wUlil 
then be satisfied. 

We may collect the results of the above investigation into J 
the following Rule : 

To obtain an expression for s as a function of the ind&i 
pendent variables x^, ...x^, which shall satisfy the f 
aifiTereatial equation 

F{x^, ... »„, p^,...p,) =0, 

*iid shall when x„ is made equal to (or to any numerical 
Talne) reduce to a given function of a:,, ...x^^, which we will 
Represent by ^ (ar, , ... a;^,) . 

HtlLE. Construct the linear partial differential equatioi 

'"'W-Ci lip* <ipi f^^iJ ' 
and forming its auxiliary Lagrangean system deduce its in- 
tegrals 

<^, = e,,... 0fc 1 = Chi 1. 
in addition to the known particular integral F=0. 

Between the above integrals and the equations 

d^ix„...x._,) . _ '^tfe--0 

^' rf.c, ' ■" ^""' dx.^^ 

eliminate, after making x^ = 0, the quantities 
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In the resulting n — 1 equations replace 

^2 by <^3, ... c^^ by ^^^y 

and we shall have a system of equations which with i^= 
will determine values of^^, ...pn> which will render 

dz —jp^dx^ — ... ^pndx^ 

an exact differential. The integration of this will give th6 
integral sought. 

In the case in which the given partial differential equation 
is of the form 

F{x^,...x^, z, p,, ...^J = 0, 
z being contained explicitly, the linear equation to be solved is 
^ ._, (fdF . dF\ dP dFfdP . dF\l ^ 

and the argument by which it is shewn that the integrals of 
this to be employed in conjunction with F= for the deter- 
mination of p^, ... pn need only be so conditioned as to make 
p., ...pn.^ differential coefficients of one and the same function 
ot a?j,...cp^j when a?„=0 is in character the same as that 
already developed in the present Article. It is only necessary 

to substitute in its exposition j— + A3- for -j— , and so 

dXi dz axi 

for the other functions. 

But as the auxiliary system 

dxj^ _ _ dxn dz 

dp, dp^ ^'dp, •;• ^Vjp« 

^A dp^ , . 

dx^'^P'dz . dx^'^^'^dz 
virtually includes the equation 

dz —p^dx, — ... --pndXn^O, 
the ultimate expression of the Eule will be as follows: 
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To obtiun an expresaion for ^ aa a function of a 
which shall satisfy t. 



F{x^, ... x„ z, p,, ...p,) = 0, 

I Mid shall when x. is made equal to (or to any particular 
L constant yalue) reduce to a given function ■^ {x^, ... x^^ of 
Fthe independent yarialiles a;,, ... x^^. 



\ lie the 2n — 1 integrals of the auxiliary system (12) which ai 
V'ttdditional to the particular integral F= 0. Make in these 2 
T aquations a^ = and forming the further equations 

« = i/r {x^, ... a-^,), 

_ d^ (a;, , ... X, .,) 



_d^{x^ , ... x^ 



eliminate the 2w quantities ic, , ... x„_-, , z, p,, ... p„. AVe thus 
obtain n equations among the constants c^,...Cm- 

Substitute in these equations i^, for c,, ... ^, for c„, and 
wo have n equations connecting x„ ... x^, z, p^,... p„, from 
which with the aid of the given equation p,, ... p^ may be 
eliminated, and tliere will result a single equation connecting 
»,, ... x^ with z. This is the integral sought. 

[It appears from the manuscript that an example was to 
have been supplied here,] 

5. Cauchy's method is evidently a general one. But its 
generality is not of the same kind as that which belongs 
to Lagrange's solution of linear partial differential equations. 
It conducts us, not to a form embracing every possible 



^ 
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solution, but to a system of results from which every possibit 
Bolutiou may be derived, by arbitrarily varyin<; the fonu of 
the thnction which expreases the initial state of tlie dependent 
variable, that is the value of z when a:, = 0, and then per- 
forming certain elirainationa. To obtain a complete primi- 
tive we should only have to assume aa the form of z when 
2r„ = a function of the variables x^ ... x^, involving ii 
independent constants. The form of this function is arbitrary. 
Each distinct determination of it under the conditions leads 
to a distinct complete primitive. The number of such com- 
plete primitives is infinite. 

There are some most important problems in which the 
knowledge of a single complete primitive is all that is re- 
quired. For this purpose the method of Jacobi wliich wc 
shall now give may be employed. 



Jucohi's Last Method, 

6. Supposing 3 to be not explicitly involved in the given 
partial differential equation 

i^, (j;,, ... x,,p,,p„) = 0, 

which we shall as before represent by F, = 0, the problem of 
the discovery of a complete primitive consists in the finding 
of n — 1 equations 

i^,= «„ F^ = a,„ 

such that between any two functions F^F^ the relation 

KJ^] = « (1) 

shall be identically satisfied, The valuer of ^,, ... p„ deducei! 
from the equations, by rendering 

d3—p^dx^ — ...—p,^dx„=0 

integrable lead us to the complole jrlniitjve expressed by its 
integral. 
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Now the idea upon which Jacobi's later methods rest is 
that of directly solving the different systems of linear partial 
differential equations flowing from the general condition (1), 
not of solving, as in Cauchy's method, one of those equations 
and then limiting that solution by conditions which virtually 
involve the satisfaction of the others. 

It is evident that the entire series of — ^^ conditions 

(1) will be satisfied if we determine F^ to satisfy the single 
equation 

then F^ to satisfy the system of two simultaneous partial 
differential equations 

[i^,ig = 0, [i^,i^3] = 0, 

then F^ to satisfy the system of three simultaneous partial 
differential equations 

[f,f;\=o, [F^^=0, [F,F^^0, 

and so on, until finally F^ is determined by the solution of 
the system of n — 1 partial differential equations 

[F,F,] = 0, [F,F,] = 0, [F^,F:\ = 0. 

Now all these are particular cases of the general problem 
of determining a function P which shall satisfy simultaneously 
the equations 

[i^,P] = 0, [i^,P] = 0, [i^„P] = 0. (2) 

F, jP,, ... F„ being given functions between each pair of 
which the equation 

IF,f;\ = 

is identically satisfied. Here F will represent in succession 
the series i^, F^, ... F^. 

The given system is one of homogeneous linear partial 
differential equations. It belongs to the class of systems the 

8—2 
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general theory of which is discussed in Chap. xxvi. But it 
is not necessary to apply the theory in its ffeneral form. We 
need only a single integral ; for a single valuei of each of the 
functions i^j, jFj, ... jF„ suffices in combination with the given 
value of i\ for the determination of a complete primitive. 
Now it may be shewn that the system is oi the class dis- 
cussed in Chapter xxVi. If expressed symbolically in the 
form 

A,P=0, A,P=0,...A,P = 0, 

the condition 

(AA-AA)i'=0, 

will be identic$.lly satisfied. Hence Jacobi's method for the 
treatment of systems of this kind may^be applied. 

That the system ia of the kind asserted is a consequence of 
the following proposition. 

Proposition. If the equations 

[wP]=0, [t?P] = 
are expressed in the symbolic form 

AP=:0, A'P=0, 

then the derived equation 

(AA'-A'A)P=0 (3), 

will be equivalent to 



\[uv'\ P\ = 0. 

" *"^ \dxi dpi dpi dxj * 
*"^ \dxi dpi dpi dxJ * 
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Qce -,— is the coefficient of 5- in A'P, and -j~- 

B coefficient in AP, its co efficient in the derivedeiiuation (3) 
llrill be {Chap. xxv. Art. 6), 



^i^nfdu d'v da d^v dv d'u dv d'u \ 
\dxi dpidxj dpi dxidxj dx^ dpidXj d-pi dxtdxj 

_d^^i=nfdu dv du dv'S 
dxj *"' \cixf dpi dpi dxj ' 



-M. 



TT / » . ' . ' . 1 11 *i=B fd \uv\ dp d\uv\ dP\ 

Hence (AA - A A) P = S^:? ( — ^--^ -. \^ -,- 

■* ^ dXj dpj dpj axjj 

-[w-p], 

whence the Proposition is established- 
Applying this to the system (2) we see that any derived 
equation will be of the form 

But [J^-PJ] = by the conditions .given ; hence the condi- 
tion (A,Aj — AjAj) P = 0, is identically satisfied. 

The results of Chapter XXVi. being thus directly applicable 
to the system under consideration, we see that a common 
integral of the system (2) may be found by a series of alter- 
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nate processes of integration and derivation. We begin by 
seeking an integral of the first partial differential equation. 
By a process of derivation, always possible, followed by the 
integration of a differential equation between two variables, 
we arrive at a common integral of the first two partial di£Fe- 
rential equations. Again, by a process of derivation followed 
by the solution of a differential equation we obtain a common 
integral of the first three partial differential equations. And 
so on, until a common integral of all is obtained. 

7. Another solution of the above problem has recently 
been given. Beginning as in Jacobi's method by finding an 
integral of the first partial differential equation, a process of 
derivation agreeing in principle with Jacobi's, only more 
extended, may lead us without further integration to a point 
at which the discovery of a common integral of the entire 
system will depend only upon the solution of a single diffe- 
rential equation of the first order susceptible of being made 
integrable by a factor. Failing this, it will enable us to 
convert the given system of partial differential equations into 
a new system possessing the same general character, but con- 
taining one equation less. Upon this the same process may 
be tried with a similar final alternative — and so on till the 
required integral is discovered. [On the Differential Equa- 
tions of Dynamics. Philosophical Transactions^ 1863). 
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CHAPTER XXVIII. 



PARTIAL DIFFERENTIAL EQUATIONS OF THE SECOXD OJtUEK. 



( 

[This Chapter is a reconstruction on a larger scale of part 
of Chapter xv. At the end of the Chapter reference will 
1)6 given to other writings of Professor Boole on the subject 
liere discussed.] 

!• The general form of a partial differential equation of 
the second order is 

F{x,y,z,p,q,r,8, «) = (1), 

where 

dz dz d^z d^z d^z 



P^ZTZ^ ? = x.» ^ = :xj«' ^ = T;::7n» ^ = 



dx' ^ dy' da?' dxdy' df' 

It is only in particular cases that the equation admits of 
integration, and the most important is that in which the 
differential coefficients of the second order present them- 
selves only in the first degree ; the equation thus assuming 
the form 

Rr^-Ss^- Tt= V (2), 

in which -B, 8j T, and Fare functions of a?, y, z,p and y. 

The most important part of the theory of the solution of 
this equation is due to Monge, and was extended hy Ampfere 
to the more general equation 

Iir+S8+ Tt^- U{8^-rt)=-V (3). 
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This equation, together iritli tLe particular equation of 
Monge, and the equation 

Iir + Ss + Tt+ D"(s'-r(} = 0, 

botU which though falling under Ampere's general form 
possess peculiarities demanding special notice, I propose to 
consider in this Chapter. 1 shall in conclusion make soroe 
observations on the theory of partial differential equations of 
the second order with more than two indepeudent variahlea. 

Monge'a method, and Anipfere'a in bo far as it is an exten- 
sion of Monge's, consists in a certain procedure for discovering 
either one or two first integrals of the form 

»=/(•) w. 

u and V being determinate functions of x, y, s, p, and g ; and 
f being an arbitrary functional symbol, i rom these first in- 
tegrals, singly or in combination, the second integral involving 
two arbitrary functions is obtained by a subsequent inte- 
gration. 

Now this procedure involves the assumption that the pro- 
posed equation admits of a first integral of the form (4). But 
such is not always the case. There exist primitive equations 
involving two arbitrary functions, from which by proceeding 
to a second differentiation both functions may be eliminated 
and an equation of the form (2) obtained, but from which it 
is impossible to eliminate one function only so as to lead to an 
intermediate equation of the form (4). Especially this haii- 
pens if the primitive involve an arbitrary function and its 
derived function together. Thus the primitive 

leads to the partial differential equation of the second order 



1 



..(G), 



but not through an intermediate equation of the form (4). 

It is necessary therefore, not only to consider the ( 
which the assumed condition is satisfied, but also to 



ase in H 
notice II 

■J 
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what has been done in those cases which do not at present 
fall under the dominion of any known method. 



Oenesis of the Equation* 

2. Prop. i. A partial differential eqtiation of the first 
order of the form u =/(v), or its symmetrical equivalent^ 

F{u,v)=0, 

in which u and v are any functions of a?, y, z^ p^ q, always 
leads to a partial differential equation of the form 

Er+Ss+ Tt+U(/-rt) = V. 

For, differentiating the proposed first integral with respect to 
Xy and with respect to y, we have 



dF/dfi du du ,^u\ 
du \dx dz^ dp dq J 



dFfdv dv dv dv \_ 
dv \dx dz^ dp dq )~^ ^ 



dF (du du du .du\ 
du \dy dz^ dp dq J 



dF/dv dv dv , ^^ .\ /v 
dv \du dz^ dp dq ) 



\dy dz^ dp dq 

For brevity, write 

fdu\ ^ du du , [du\ j, du ^ du 

and then eliminating 

dF dF 

du' dv ^ 
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We have 

'fdu\ du du ) {fdv\ dv dv X 

{fdu\ da . du \ {fdv\ ■ dv dv ) 

which, on effecting the multiplication, gives 

{du (dv\ /du\ dv] 
dp \dyj \dyj dp\ 



+ 



(/^\ ^ _ ^ fdv\ j_ fdi£\ dv du /dv\^ 
[[dx) dp dp \dx) \dy) dq dc[ [dyj. 



s 



{(du\ dv ^ du (dv\\ 
[\(lx) dg dq \dx)) 

fdu dv du dv\ . , . 
\dq dp dp dq) ^ 

_ /du\ /dv\ ^ /du\ /dv\ . . 

^\Ty)\d^) \Jx)\dy) ^"' 

a result which, since u and v are by hypothesis given func- 
tions of a?, y, «, p, J, is seen to be a particular case of the 
general form (3). 

We may hence deduce also the conditions under which 
particular forms included in the general form (3) arise. Thus, 
in order that the equation u =f{v) may give rise to a par- 
tial differential equation of the second order of Monge's form 

it is necessary that the condition 

du dv du <?^_^ 
dq dp dp dq 

should be identically satisfied. This requires, by Chap. n. 



Art. 1, that u and v, considered 
not be independent. 
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functions oi p and q, sliould' 



3. The geometrical relations of the equation (3) are 
Temarkable. It may in particular be shewn that an equation 
of this form will be satisfied by the equation of any surface 
■which constitutes the envelope of any system of surfaces 
formed by the variation of three parameters in subjection 
to two arbitrary conditions. , For let the common equation of 
the enveloped surfaces be 

^=f{x,y,a,b,c) (8), 

the parameters a, h, c varying in subjection to the conditions 1 

^, (o, 5, c) =0, tpj {a, h, c) = 0, 

iditions which, determining h and c as functions of a, may I 
i)C reduced to the form 



I 

° 



H"), 



=+(■>)•■ 



■ (9). 



Now the values of p and g being tlic same for any point.! 
in the envelope as for the same point in the generating surface] ■ 
!ve have for all such points 

_ df(x, y. a. h. c) _ df{x, y, a, i, c) 



(!»)■ 

) enable us 1 



dx 

These two equations in conjunction with 
itermine a, b, c as functions of a;, y, z, p, 
dues be 

o= w, h=v, c = w. 

Then substitullng in (9) we have 

equations which hold for all such points. These are then the 1 
partial differential equations of the first order of the envelope* I 

Now each of these equations is of the genera! form (4) j I 
*hence by Prop. i. the partial differential equation of tlie j 
Second order is of the form (3), as was to be proved. 
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Let US actually construct this equation. 

Differentiating the first of the equations (10) with respect 
to X and to y, and regarding therein a as a function of tuose 
variables, and b and c as functions of a, we have 

^^d'f / dy d\f dh ay dc\da 

dx^ \dadx dbdx da dcdx da) dx ' 

3= ^y (dy d^f dh dy dc\da 
dxdy \dadx dbdx da dcdx da) dy ' 

from which we readily derive 

/ d\f\da f ^y^^^O 
V dix? ) dy \ docdy) dx 

Proceeding in the same way with the second equation of 
the system (10) we have 

/ dy \ da ( dy\ da __ 
\ dxdy) dy \ dy* ) dx 

Hence, eliminating -p and -^ , we have 

(-=S)'-('-i)('-f)- 

df dxdy da? dj? dy* \dxdy) ' 

the equation sought. 

Comparing this with the general form (3) we have the equa- 
tions 

dy* dxdy _ tfas* _ 1^ _ da? djf \dxdyl 

IT" ST^ -T'U' V ' 



r.4.] 
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wc amve at the! 



■hence eliminating ,j, ■, j, and t'j' 
equation, 

S'-i-iiUV-IiT)=0. 

This then is the condition which must be satisfied in order 
that the equation (3) may admit of an integral representing 
the envelope of a system of surfaces in which three parameters 
vary in subjeetiou to two connecting conditions. It ia only 

roved however to be a necessarjf, not to be a su^cient, con- 

.tion. 



Solution of the equation Br + Ssi-Tt+ [/'{s' - )■() = V, when 
a first integral of Oie form F[u, v) =0, exists. 

In the following sections we propose 
t. To shew that when a first integral of the above form 
s, its discovery depends upon the solution of two simul- 
taneous partial differential equations of the first order re- 
solvable into linear equations. 

2nd]y. To shew how fi^m such first integral or integrals 
the second integral is to be obtained. 



Prop. II. If the equation 

Iir+Ss+Tt+U{s''-rt)=V 

admit of a first integral of the firm J^(u, v)=0, in which w 
mtdv are functions of x, y, z,p, q, then will F[u,v) considered 
as a function ofx, y, z, p, q, and represented as such for hrevtty 
hy F satisfy the two partial differential equations of the first 



i 



, !dF\ dF 
\dx/ dq 



--(f)f--(f)('f) 
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\dq/ dq dp \dpj 

, ^ (fdF\ dF . fdF\ dF\ ^ 

(dF\_dF dF (^\^^. ^ 
^dx/^dx -^ dz^ \dy ) dy ^ dz ' 



in which 



Kegarding the function F in the proposed integral F= 
simply as a function of i», y, z, p, q, we have 



^dF\ dF dF 



^ 



dp ' ' dq 
dF dF 



= 



\dy J dp dq 



> 



dp dq 







(11). 



On the other hand, regarding J?' as a function of x, y, z, y, q, 
mediately through u and v, we have the system 



^dx) dp 



dF 
du 

dF 

du \\dy) ' dp 



dq ) dv \\dxj dp dq ) 



= 0, 



'du\ du 
dy) dp 



8+ — t\ + 



du 
dq 



dF (/dv 



dv \\dy 



((dv\dv^^dv)^ 
\\dyj dp dq J 



> (12), 



0, 



and these systems are equivalent 

dF 
Now if from the second of these systems we eliminate -v- 

dF 

and -J- , we obtain (Art. 2), a result which must be equiva- 
lent to the proposed partial differential equation, 

5r + i&+r^+Z7(5«-rO = F (13). 

This equation then considered as a relation between r, *, t, 
must be an algebraical consequence of the relations (12), and 
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therefore of the equations (11). If then we determine alge- 
braically two of the quantities r, s, t, (we select r, t) from the 
system, and substitute their values in (13), that equation 
ought to be satisfied independently of the value of the re- 
maining quantity s. Now supposing^ and q to be both con- 
tained in F, so that neither -7- nor ,— vanish, we have 
from (11), 



r = — 



(dF\ dF 
\dx) dq ' 

dp 



_^S^2/) 



< = 



dF 
dp 



dF 
dq 



substituting which in (13) there results 

R fil\ ^+ t(—\ - + U(—] (— 
\dx J dq \dy J dp \dx J \dy 



dF\ dF dF 

) dp dq 



I \dq/ dq dp \dp J 



+ 



^(f)S^''(f)|}"»-(»)- 



Now as this equation is to be satisfied in virtue of the con- 
stitution of jB, S^ Ty U, F, and the function F, and indepen- 
dently of 8y both the coefficient of s and the absolute term not 
containing 8 must be separately equated to 0. Thus F con- 
sidered as a function of x, y, z^ p, q, and containing p, q, at 
least must satisfy the partial differential equations 

„ /dF\ dF r^ fdF\ dF 
\dxj dq \dy J dp 



+ 



-(SGf)-"^^^- 



\dx J \dy 



\diJ 



dp dq \dp J 



+ T 



dp dq 

dF\' 



i- 



dF\ dF fdF\ dF 




^dx I dp \dy J dq 



}-oJ 



(15). 
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This result may also be established by forming the equa- 
tions of condition which expresa the proportionality of 
R, 8, ... V, to the con'esponding quantities in trie constracted 
equation (7). From these equations of condition it is actually 
possible to eliminate in two distinct ways the quantities 

f-j-J , f — 1 , -J-, J- 1 the result being the formation of two 

partial differential equations for u agreeing in form with those 
above given fOr F, (See the memoir Ueber die partielle Diffe- 
rentiaJgleichung... Crelh's Journal, Vol, 61.) The actual 
transition from the former to the latter rests upon the con- 
sideration that the equation F {u, v) = 0, when F is arbitrary, 
is not really less general than the fonn ^ [F (u, v), v\ = 0, in 
which the <& is arbitrary. And here u has been replaced by 

F(», .). 

The only condition respecting the application of the above 
equations is tliat we do not admit any relations which maiie 

eitiier -;— or -j- to vanish. 
op dq 

5. Phop. Ill, The solution of the system of partial dt^i- 
rential equations established in the last proposition may in a2 
I'.ases he made to depend upon that of simultaiieous linear parHoi 
differential equations of the first order. 

In demonstrating this proposition we shall consider first 
the ease in which U= 0, then the case in which V= 0, laa^j 
the case in which neither of these quantities vanishes. Tm 
ground of this division will appear in the investigation. 

Case 1, Suppose U= 0. Tlie equation then is of Monge's 
form. 



The second equation of the system (15) tecomea 

„ldF\' dFdF „ldF\- 
■"teJ -'^d} d^ + ^Kd^j-"' 
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and therefore breaks up into the equations 

dq ^dp ^ ^i ^dp ' 

i»j and m, being the roots of the quadratic equation • 

5m'-/8^m+ !r=0 (16). 

As each of the above constituent equations is of the form 

dF^ dF 
dq " dp'* 

the system (15) may be reduced to the form 

dFdF 



\dx J dp \dy J dp dp dp 



= 0, 



which breaks up into the equations 



dF ^ ^ fdF\ m(dF\ ^. dF ^ 



The former of these we must reject (Art. 4). There le- 
mains for the determination of F the system of linear partial 
differential equations 



dF_ dF^ -^ 
dq dp " 



MS)-Kf)-''-f=« 



(17), 



and there will exist either one or two systems included under 
this form, according as the roots of the quadratic (16) are 
equal or unequal, 

B.D.E. II. 9 
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Pose II. Let F=0. The system (15) then becomes 

\aqj dq dp \dp J 

]\dxjdp \dy) dq) 
Eliminate Uhj multiplying the first equation by 

\dx/ dp \dy ) dq ' 

the second by 

(dF\ (dF\ 
\dxj\dyj' 

and subtracting; we obtain, after rejection of the common 

. dFdF 
factor -J- -T- J 

' ■^(S)'-«©(f)--(f)-o- 

We shall put this equation in the place of the second equa- 
tion of the system. This we are permitted to do under the 
restriction that in seeking to satisfy the system so changed 
we do not make use of any relations which would cause either 
of the two factors employed -in the process of elimination to 
vanish or become infinite. 

The new equation reduces to one equation, or breaks up 
into two equations of the form 

. ©-^-» (-). 

m being determined by the quadratic equation 
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Uaking (t— ) = *"{ j— ) "i tte first equation of the system 
1(15), we get 

V which breaks up into 



^"■f+'-f^*©-- 



But if we combiae the first of these with (18), we obtain 



(f)=». ©-. 



I and this combination causing botli tte factors empl _ 
[the elimination of U to vanisli must be rejected. There 
3 then tile combination 



_ dF , „dF ^ „ /dF\ 
di dr 



•\dyl 
'(f)" 



..(19), 



and this will represent either one or two systems of equations 
according as the quadratic determining m, has equal or un- 
eqoal roots. 

Caae III. Let neither U= nor V= 0. 

Multiply tlie second equation of the system (15) by an 
indeterminate quantity I, and add to the first ; then we have 



\dxj dp \dy ) dp 

.dPdF 
' dp dq 



Hv-m)'^^=o.. 



..(20). 
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We shall enquire whether It is possible so to determine I 
as to resolve this into linear factors. 

We might investigate thia by resolving the. equation as 
a quadratic with respect to -y- or -7—. But /the form of 
the equation suggests what the forms of the linear factors 
must be if the resolution be possible. For as the squares of j- 

and -v- both appear, and these squares alone, in the fime- 

tion to be resolved, it is clear that -7- and 3— will be the 

dq dp 

only differential coefficients of F. which will appear in both 

linear factors in common. The most general supposition 

possible is then that one factor shall contain j- and -j- with 
( -7- j , the other the same with ( ^ ) . 

Assuming then one factor to be of the form 

I ^+ ni — 4- {—] 
dq dp \dx J ' 

it is seen from the. form of the coefficients of the first three 
terms of (20) that the other factor must be of the form 

jl^.TldJ^U/dF\ 
dq m dp n \dy)^ 

and the resolved form of (20) must be 

\ dq m dp n \dy )] \ dq dp \CLx)) ^ 

Multiplying out and equating coefficients, we obtain the 
conditions 
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Tin 






rp_ Urn 
n 

n 



TP 



F-)SZ = jBw + 



m 



The third and fourth of these conditions are equivalent, 

and give n = l. The first and second are also equivalent, 

T 
and give m = jy. These values reduce the last equation of 

condition to 

80 that I is determined by a quadratic. The resolved form 
of equation (20) now becomes 

To these results we may give a somewhat simpler form by 
making Ul=m; not the m used above. We have then as 
the quadratic for determining m, 

m*-v8m + BT- UV= (21), 

and as the resolved form of (20), 
(^dF dF rrfdFWi dF ^rpdF ^j^(dF\\ 
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Then we liave from the 



Let m, and m, te the values of 
last the two diBtinct equations 






-(£))=«. 



{«f-.f-Kf)} 



■.f--f--©} = 



and it is evident that these will be t 
equations (15) from which they were 



gether equivalent to the 
lerived. 



Now to satisfy these equations simultaneously it is nece»- 
sary that we should equate to one linear factor from each of 
tlieir first memhera. If we equate to the first linear factorfl, 
we have 







-%-' 




<i)- 


■0. 






-f- 


dF 


o- 


= 0; 


whence 


hy 


subtraction, 












K 


-"^f.- 





This combination must therefore be rejected (Art. 4). For 
the same reason must the combination formed by equating lo 
the second linear factors in the left-hand members of the 
above two equations be rejected. There remains then only 
the combinations formed by equating to the first factor of 
one of these members, and the second. of the other. 



Thus we should have the combination 
^dF , dF , 



dq '■dp 



■^m- 



..(22), 
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with the combination which would be obtained from this by 
interchanging m, and m,, 

6. It results from the foregoing investigations that the 
function F is in all caaea to be determined by the solution of 
two simultaneous linear partial differential equations with 
Jive independent variables. Now tlie theory developed in 
Chapter xxv. shews that the number of integrals of such a 
ayatem cannot exceed three. That theory enables us both to 
determine what the number of integrals is, and to construct 
the system of ordinary differential equations, reducible to the 
exact form, upon which their discovery depends. 

We have seen that the knowledge of two integrals u = a, 
w = J of the system enables us to construct a general first 
integral 

F{u,v)=0, 

of the partial differential equation (3). And the solution of 
this first integral would lead us to the second integral which 
is the final object sought. But the direct solving of a partial 
diiferential equation of the first order which is not linear and 
which involves in its actual expression an arbitraiy function 
is difficult, and happily it may be avoided here, Tlie fol- 
lowing propositions will enable us to accomplish the virtual 
solution by a different solution, founded however upon the 
Bame general principles. 



7. Prop. iv. The integrals of the respective systems of 
aimultajteous linear partial differential equations vpon which 
the determinaiion ojF depends are so related that tf from tico 
siich respective integrals the values of p and q are determined, 

they will retider the equation 

dz =pdx + qdg 

integrohle. And in the particular case in which the two systems 
become identical, any two integrals of the system stand in the 
game relation. 



^ 
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For, let ^ be an integral of the system (22), and "S? an 
integral of the associatea system obtained by interchanging 
971 and m^ in the case in which these quantities are different. 
Tnen ^ satisfies the equations 






and ^ satisfies the equations 









But the necessary and suflScient condition in order that the 
values of p and q derived from the equations ^ = 0, "^ = 0, 
may render dz ^pdx — qdy integrable> is 



fd^\ d^__d^ fdT\ 
\dxj dp dp \dxj 



-(f) f-f©- <->• 



See Chap. xiv. Art. 11, Equation (36). 

Now if from the previous equations we determine the values 
of 



'^\ f^\ fd^ fd^ 
- ' \dy)' \dx)' \dy)' 



\dx)' 



and substitute them in the above equation of condition it will 
be identically satisfied. 
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The determination of [-f-] > ••• from the previous systems 

requires that U should not vanish. Hence the proposition is 
established except in the case of Z7=0, which is left doubtful. 

To examine this case let us revert to the system (17) which 
is proper to it. To that system since 

whence Hm^rn^ = T, 

we may give the form 

dF dF ^ [dF\ , (dF\ V dF ^ 

or the form obtained from this by interchanging m^ and w,. 

Substituting in these respective forms ^ and "¥ in succes- 
sion for F, we find 

d^_ ^ fd^\__ fd^\ V d^ 
di'^'^'dp' \dx)^ '^\dy) B dp' 

d9_ d^ fd^ _ fd^ V d^ 

dg^'^^dp ' \dxj^^'''^\dyl^ R dp' 

and these values substituted in (23) reduce it to an identity. 
Thus the proposition is established generally. 

Lastly, as in the case in which the two roots of the quad- 
ratic for determining m are equal, the two systems of partial 
differential equations for determining ^ and "^ become one, it 
follows that if from two integrals of that one system we can 
deduce values of ^ and q these values will render the equation 

dz^pdx-^qdy^Q 
integrable. 
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8. Prop. v. When (he system of simultaneous linear par- 
tial differential equations determining F admits of two integrals 
u = a, v = b, it will admit or teiU not admit of a third inte- 
gral w = c, according as the roots of the quadratic determining 

m are equal or unequal. 

The system in question, (22), IwcomeB when we divide by 
Uaxii write for [~-j~] and (-5-) their full 



dF 
ds 


^»f^ 


m^dF MdF 
V df'^ (J £1° 


= 0, 


dF dF 
d.+fdi* 


T dF m, dF 
U df U di 


= 0, 




or A,F=0, 


A,-F=0, 




wliicli 










A. 


=1- 


d m. 

d^-^V 


d B 
df^U 


d 
d, 


\ 


-i^^ 


d ^T d m, 

dz^Udi^TJ 


d 
Hi 


Hence the equation 










{A,A 


,-A,A 


)F~0 





m,-m,dF /.T. wSdF I, m, . S\dF „ 

In this expression the coefficient of the first term only haa 
heen calculated. 

Now. by the theory developed in Chap. sxv. in order tliat 
the two simultaneouB partial differential equations ahoulil 
have their full complement of integrals (three) it is necessary 
that the above equation should be satisfied identically. This 
involves three conditions, namely, 
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Wj — W, = 0, 

the first of wliich is the one affirmed in the Proposition to be 
necessary. 

Secondly, it is to be shewn that if this condition be satis- 
fied and if the system of given linear equations admit of two 
integrals m — a, t? = J, it will admit of a third. 

Beplacing m^ and m^ by m the system becomes 

dF dF mdF B dF_ 
dy^^dz'^U dp^ U dq^ ' 

dF dF TdF m dF _ 
dx^^ dz^ V dp^ U dq^ 

Now if we construct from this the corresponding system of 
ordinary differential equations, we shall find it to be 

dz ^pdx — qdy = 0, 
dp- ^&J--j^rfy = 0, 

dq-^~dx- ^dy = 0. 

Now it is impossible that the first of these equations should 
b(& integrated without a previous determination of p and </ as 
functions oix^ y, z^ seeing that dx, dy^ dz Ate the three differ* 
entials entering into that equation. Such determination can 
only come firom the int^ration of the second and third equa^ 
tions of the system* But if these equations can be integrated 
in the forms tf = a, t; = &9 then n BBdv being particular valoeft 
of J^ satisfying the partial differential equations, it Mhmn 
firom the last Fropomtiou that the values <^ p and q which 
they will yield will make the first equation ^ inieffrable. 
Hence if the system admits of two integrals it will aamit of 
three; as was to be shewn. On the basis of these Froposi^ 
tions the theory cf the second initgntiim tt^iB* 



i 
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Theory of the Second Integraiion. 

9. First suppose the values of m uneq^ual. 

Then Wi = a, , Di = J, being the two integrals (and we have 
seen that there cannot be more than two) ot one of the systems 
of linear partial differential eqxiationa, and Wj=aj, T!i = i, tliose 
of the other, tlic general lirst integrals of the given system 
will be 

*(«.> fi) = 0. ^(«». v,)=0. 
The values of p and j determined &om these will by 
Proposition IV. render 

dz —pdx- qdij = Q 

integrable, and the integral of this will be the general integral 
of the proposed partial differential equation. For it will in- 
volve explicitly or implicitly two arbitraiy functions derived 
from those in the first integrals. 

It suffices however, following herein Charpit's method, to 
combine one general first integral derived from the one sys- 
tem with a particular first integral derived from the other 
system, e.g. the integrals 

"^ (m,, «J = 0, Wj = a. 

The values of p and q hence derived, and employed as 
before, will lead to a second integral involving one arbitrary 
function and containing two arbitrary constants. This con- 
stitutes a complete primitive from which the general solution 
will be obtained by converting one of the arbitrary constant!) 
into an arbitrary function of the other, and eliminating iht- 
latter between the equation and the one derived from it by 
differentiation with respect to that constant. 

Secondly, suppose the valnes of m equal. 

In this case we have but one system of partial differential 
equations so constituted however that if it admits of two inlc- 
grala it will admit of three. 
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Letu^a, 'B = h, io=c represent these integrala. Then if'l 
these we eliminate p and q we shall obtain a final into- \ 
of the form 

iiicl this constitutes a complete prlmilive from which wo shall 
ieducc the general integral by making i = (^(a), c =i/r (a), 
kuiii eliminating a between the equations 



To prove this let ua combine the general and particular I 
first integrals 



The values of p and q hence obtained make 

ds — pdx — qdi/ = 

□tegrable, and the result can be no other than the remaining 
Ktegral v> = c, or rather what this would become on eliminat- 
ngp and q from it. But since the equations by which this . 
Iitegration are to be effected are equivalent to 

« = o, « = *(a), 

will become a function of x, y, s, a and (a). Also by 1 
Jharpit's method c is to be treated as a function of a, bo that I 
mtimately we have the result above assigned. 

We have here supposed U not to vanish. If it do the | 
theory assumes another but simpler form. Let 



»=/(«), 



'•f-W 



L 



be the two general first integrals. Then, since by the Ci 
dition at the close of Art. 2, if p be eliminated from these ! 
equations q will also disappear, it suffices to eliminate them | 
together in order to obtain the general secoad integral. 
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10. Although the cases in which U= and V= have in 
the foregoing sections been treated for simplicity apart, their 
theory might have been deduced from that of the case in 
which neither P'nor Fvaniahea. 

Thus to deduce the equations for the case of U= elimi- 
nate from the general system (22) -j- and -j- in aucces- 
Bion, and we find 

But from {21} RT- m,m, = UV. 

Sutatitating, and then dividing by Z7we find 
ydF_ /dF\ jffdF\_ 
dp ^\dy) \dx j ' 

the equation determining m„ m, being 

This is equivalent to the results of Art, 5, Case r. 

11. We found it necessary (Art. 3) in order that the gene- 
ral partial differential equation of this Chapter should be satis- 
fied by the envelope of ^ system of surfaces the equations of 
which contain three pai'ameters varying under two conditions 
that the relation 

S'-i-iiUV-BTi^O 

should be satisfied. 

It appears from Art. 8 that this is but one of three condi- 
tions necessary and together sufficient for this purpose, The 
formal conditions for every form of ultimate solution con- 
sistent with the existence of a general first integral i'"(u,D)=0 
can be deduced in the same way. 
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[In tie Sulletin de VAcadSmie Impiriale dee Sciences de 

't Pglerabourg, Vol. IV. 1862, there is an article entitled Gon- 

I nitrations sur la recherche dea integraks premieres des Sqy.a~ 

I Hona diffSrentielhs parliellea du second ordre, par Q, Boldt 

I (La Ic 7 Join 1861). 

The article occupies pages 198 — 215 of the voltune. Al- 
' though the name does not quite correspond, I consider that to 
be a misprint, and I attribute the article to Professor Boole, 
partly from the nature of the contents, and partly hecauae it 
IS known by his friends that he was engaged at a time corre- 
eponding to the date here given in the preparation of a mathe- 
matical article in French. 

The object of the article is to determine the conditions 
Ljaecessary for the existence of a first integral of the equation 



dix'^ 



axdy ay 



where B, S, T, and IF are any functions of x, y, 
Bvid G 



dz , dz 

;, -J- and — ; 

ax ay 

\ also to determine the conditions which must hold in order 
P'that Ampfere's method of integration may be employed. 

In Crelle's Journal, Vol. LXi. there is an article by Pro- 
■ JeBBor Boole, entitled t/eier die partielle Differentialgleichung 
|*wetter Ordnung Hr+Ss+Tt+U {a' - rt) = V. 

The article is dated 1862; it occupies pages 809—333 of the 



I 



Among Professor Boole's mannscripts I found a memoir 
■ closely resembling the article in Crelle'a Journal; it 
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would appear that the memoir was drawn up with a view to 
publication in the Transactions of some English Scientific 
Society, and that this design was afterwards abandoned in 
favour of the article in Crelle's Journal. 

After some hesitation I have resolved to print this memoir. 
Even if the memoir had been identical with the article in 
Crelle's Journal it would have been convenient to the English 
reader to be able to avail himself of the investigations ; and 
the memoir contains remarks which do not occur in the article, 
and which are interesting in connexion with the history of the 
subject. There is some repetition of matter which has already 
been given in Chapter xxviii. ; but I was unwilling to impair 
the completeness of the memoir by abridgment or omission. 
Accordingly the memoir forms the next Chapter of the present 
volume. 

In Article 2 of the next Chapter will be found the pro- 
cess to which there is an allusion towards the end of Article 4 
of Chapter xxviii. 

It is obvious that the subject of partial differential equa- 
tions of the second order was much studied by Professor 
Boole. The chronological order of his writings on the sub- 
ject appears to be as follows : 

1. Chapter xv. of the first edition of his work. 

2. The article in the Bulletin of St Petersburg. 

3. The memoir which forms Chapter xxix. of the pre- 
sent volume. 

4. The article in Crelle's Journal. 

5. The Chapter xxviii. of the present volume.] 




( 1*5 ) 



CHAPTEE XXIX. 



iON THE SOLUrrON OF TUE PARTIAL DIFFEEENTIAL EQUATION 
Sr+8s+Tt+Uis''rt) = r, in which ii, 8, T, U, V 
ABE GIVEN FUNCTIONS OF iC, y, Z, p, q. 

1. The equation, the theoiy of the solution of -which 
I propose to consider in this paper, is remarkable from its 
connexion with Geometry. If the equation of a surface 
contain three constants which vary as parametera in sub- 
jection to any two conditions connecting them, the gene- 
rated envelope will satiaiy a partial differential equation of 
the above form. In other words any envelope of the surface 

F{x, y, z, a, h, c)=0 

formed by the variation of a, h, c in subjection to two con- 
necting conditions 

0,(u, J, c)=0, fja, I, c) = 

is necessarily an integral of a partial differential equation of 
the form given above, 

Now this theorem is the more important, because it is 
only when three parameters in the equation of a sm-face 
vary in subjection to two relations that the envelope poa- 
sesaes, in-espectively of the foi-m of the connecting relations, 
any definite character. If there be but one connecting rela- 
tion it ia possible to determine that relation so aa to make 
the envelope assume the form of any surface whatever, and 
therefore tne possible system of envelopes is in such case 
iS.D.E, 11. 10 
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unlimited. If there be three connecting relations the para- 
meters become absolutely constant and no envelope exists. 

The partial diSerential equation 

Sr + Ss+ Tt+ U{s'~rt) = V 

is remarkable also as inclading all the cases in which a 
partial differential equation of the second order admits a 
first integral of the form 

u =f{v), 

n and v being definite fimctiona of x, y, z,p, q, and/(r) 
arbitrary in form. 

Neither of these statements is sufficiently general to con- 
stitute a theory of the genesis of the partial differential equa- 
tion under consideration, but the second one is more general 
than the first, and is indeed sufficiently so to serve as the 
ground of an investigation which connects tlie solution of 
the equation in all cases with the satisfaction of a system 
of simuttaueous ordinary differential equations of the first 
order and degree. And this is the ground upon which the 
method of the paper will rest. I propose to shew, 1st that 
the solution of the given equation on the assumption that 
a first integral of the form u=f[v) exists requires the satis- 
faction of a system of two partial differential equations <rf 
t!ie first order and second degree ; 2ndly tliat this system may 
be resolved into four systems, each consisting of two partial 
differential equations of the first order and fii-st degree, two 
of which systems are irrelevant and the other two relevant; 
Srdly that the solution of the two relevant systems ulti- 
mately depends on the solution of a system of ordinary 
differential equations of the first order, and that from these 
ordinary differential equations the given equation of the 
second order may be deduced independently of the assump- 
tion above mentioned. 1 shall also discuss the theoiy of the 
second integration. And I shall exemplify another method 
of solution connected by a remarkable kw of reciprocity with 
the above method. 



^ 
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First Investigation. 

2. Prop. I. lfu=f[v) he a first integral of the equation 

Br+8s+2t'{-U{s'^rt)=V. (1), 

then wiU u and v, considered as Junctions of a:, y, z^p, j, each 
satisfy two partial differential equations of the form 

fduV 



^©'-«(i)(l)^<l) 

P /du\ du rp(du\ ^w 
\dx/ dq \dy) dp 

rrfdu\fdu\ . ^y.du du 



(2), 



tn 



, . , (du\ J fdu\ ^ J ^ du ^ du j du , du 
m witch y andi^^) stard ^ -^+p^^,and ^^ + q^^ 

respectively. 

To demonstrate this proposition we shall form directly the 
partial differential equation of the second order of which 
u=:f{v) is an integral and, comparing that equation with (1), 
deduce the conditions for the determination of u and v. 

Differentiating u =/(v), first with respect to x and secondly 
with respect to y, we have 

du du dz du dp du dq 
dx dz dx dp dx dq dx 



•^ ^ ^ \dx dz dx dp dx dq dx) ' 



du du dz du dp du dq 
dy dz dy dp dy dq dy 



^•f* ( \ I^^ ^^^ dz dv ^,dh dq) 
~ [dy dz dy dp dy dq dy) ' 



10— ^ 
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or, if we 
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^ rfw . du ^ fdv\ du du y fdv\ 
represent ^+p^ by ^ , ^+2^ by ^ , 



gbyr,gandgby«, and^by*, 



dy 



/du\ . rfw . du „, . {fdv\ , dv . dv} 






du 
op 






Eliminating jr(i?) we arrive at the partial differential equa- 
tion of the second order, 

(du (dv\ ^ dv fduVi 
\dp \dy) dp \dy)) 

(du fdv\ __ do /du\ dv /du\ ^du fdv\\ 
\dq \dy/ " dq \dy) dp \dx) dp \dxJ) 

(dv /du\ _ du fdvVi (du dv^du dv\ , , ^ . 
\dq \dx) dq \dx/) \dq dp dp dq) ^ 



-©(l)-Q(g) » 

It is seen that as respects the mode in which the qnan- 
tities r, s, t are involved this equation is of the same form 
as the given equation (1). That it may be equivalent, its 
coefficients must stand to those of (1) in a common ratio /i. 
This gives 

du fdv\ dv (du\ ^ , . 

^ i^j -^ \dy) =f^ W' 

du /dv\ ^ do fdu\ dv fdu\ du fdv\ _ -, „. 
dq [dyj dq [dyj ^^[d^)~ 'dp Wx) "^^ ^^^' 



dv /du\ du /dv 
dq \dx/ dq \dx, 



:) =f^T. 



(c), 
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du dv du dv ^ ,j 

-q ^-5 di'i"^- W' . 

^)©-©©-''^-- w- 

As we have here five equations which are homogeneous with 
respect to the four differential coeflScients of v and to /i, it is 
clear that we can, by the elimination of these quantities, 
obtain a relation connecting the differential coefficients of u 
with Rj 8, T, &c. But the peculiar cyclical form of the 
functions in the first members of the above system enables 
us to effect this elimination so as to lead to two final equa- 
tions independent of v and fi. 

Th» „d.iplyi„g (a) by (*) |. W b, (I) |, 

^^ ^y (^) (^) ' ^^ ^*^ ^^ ^ ■£' *°^ *^^'°^' ^® ^°^' ^^ 

rejecting the common factor /i, 

\dxj dq \dyj dp \dxj \dyj 

*^%%'' • (*)• 

Agan, mtaiiplTuig (a) ij (gj , Q,) by (g) Q, (o) bjr 

^j, and (e) by (^ j ^- + (^ j ^ , addmg, and agam re- 
jecting the common factor /^, we have 

-''{©i+(i)|}-« (»)• 
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Hence, u considered as a function of a?, y, », p^ q satisfies 
tlie two partial differential equations (4), (5), botn which are 
of the first order and second degree. 

As u and v enter symmetrically into the system (a), (J), 
&c., V will also satisfy two partial differential equations of 
the same form, viz. the equations 

\dxj aq \dyl dp 

\dxl \dy) dp do 



dy) ' ' dp dq 



^ ... (6). 



J 



Further, these two systems of equations constitute the 
complete system of equations resulting from the elimination of 
/A from the five equations (a), (5), (c), &c. ; for in their deter- 
mination, no factor involving either the differential coefficients 
of u and v, or the quantities -B, S^ T^ i&c. has been rejected 
directly or indirectly. 

I am not aware that the above results of elimination hare 
been noticed before. 

3. Prop. II. The system of partial differential eqtuttim 
above obtained for the determination ofu, viz. 



J. /du\ clu ^ fdu\ du 
\dx) dq \dy) dp 



+ 



KSd)-"-- 



du du 
dp dq 



= 



p. \(du\ du fdu\ du \ _ 
'^ ^ l^j ^'^\dy) dq] ~ 
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mita of resolution, into four systems, each consisting of two 
war partial differential equations of ihe first order. Of these 
m^stems two only are relevant to the solution of the problem. 

For, raultipljing tlie second by an indeterminate quantity ' 
I, and adding the result to the first, we have 

\ax/ cy dp aq \ayj dp 



Now let ua see if it is possible to determine X so aa 
make the first member of the equation resolvable into line. 
factors. We cannot say d. priori that such resolutiou is pos- 
sible as we should be able to do if that member were homo- 
geneous and of the second degree with respect to three, instead 1 
of with respect to the/owr subject variables 






ldii\ /du\ f^^\ fdu\ 
[Z:)' \d^)' [^)' [d^}' 

, du , 
-7- and -V- ■ 
dp ay 



ti 

L 0f( 



Observing that the squares of -r- and -v- are wanting ] 

In the first member of (8) while those of |t-J and I j-1 

appear, we are led to assume as the proposed equivalent of J 
that member an expression of the form 



^(^ii^-miyo-'t}- 



Multiplying the factors of this expression together and I 
then equating the coefficients with those of the first member I 
of (8) we have 
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Bm-\-\m- U+ 8\ (a), 

mm: = T\ (J), 

Xn = XF=7wn' (c), 

Bn! = B {d), 

wm'=r... (e), 

W=F (/), 

From (J), (c), (e?), we find 

T 
w=F, w' = l, w = XF, m=-y.y 

values which will be found to satisfy (e) and (/) also, and 
which reduce (a) to the form 

Px«-fifFx+i?r-crF=o. 

Supposing \ thus determined, the equation (8) becomes 

The result is a little simplified if we retain m in place of X. 
We thus find as the resolved form of the given equation 

m being determined by the quadratic 

If Wj, Wj be the values of m thus found, we have 
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and these two equations are manifestly together equal to the 
system (7). 

Now these equations can only be simultaneously satisfied 
by equating to 0, one factor in the first member of each ; and 
the different combinations which are thus possible give rise 
to four binary syBtems of linear equations. Let us examine 
these systems separately. 

If we simultaneously equate to the two first factors of 
the left-hand members of the last two equations, we have the 
systems 

a system which, when m^ and m^ are different, is reducible to 
the system 

It is clear that this cannot lead to a value of u satisfying 
the given differential equation (1), because it takes no 
account of the forms of 8, J7, and T. Indeed if we actually 
eliminate 



(du\ fdu\ 



^du\ du du 
dy) ^ dp^ dq 



from the above equations by means of the system 



/dii\ du du __ 
\dx/ dp dq " 



dp <k L (10), 



(du\ du du , 
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m 



vn. 



(derived from the assumed first integral u =f(v) by making 
/(i;) = c, and differentiating the result first with respect to a, 
then with respect to y), we find as the result 

Again, if we equate to the two last factors of the right- 
hand members of (10), we have 

which, if m^ and m^ are different, reduce to 

^O-^^'O- ©=»■ 

And it is evident that. neither are these equations consistent 
with the given equation (1), because they take no account of 
8j Z7, and B. The equation of the second degree to which 
they actually lead is 

Fr+r(5'-r«) = 0. 

There remain then the two systems formed by combining 
the first factor of each one of the first members with the 
second factor of the other, viz. 



^{ih- 






(11). 



B 



du 






•^ 



Y 



.(12). 
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That these systems are relevant to the solution of the pro- 
blem under consideration maj be shewn Toy eliminating n:om 
either of them by means of (10) the quantities 

/du\ /du\ du du 

The actual result will be 

V{Br + 88+Tt+U{8''-rt)-V}^0 (13), 

which, except in the particular case of F= 0, reduces to the 
given equation. 

More generally, if in the equation 

u =f{v) 

u and V are any distinct solutions of the system (11), the 
same result of elimination may be deduced. For v by hypo- 
thesis satisfies the equations 






Subtract these equations multiplied by/'(v) from the corre- 
sponding equations of (11), and representing u^f{v) by W, 
we have 

^(f)-.(f)-''f-. 



m. 



(f)-^(f)^^f=«. 



which being of the same form as (11) it follows that 

W=0 or w-/(v) = 

also leads to the partial differential equation of the second 
order (13). 
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4. Prop. III. To reduce the determination of the first inte- 
grals of {1) to the solution of a system of ordinary differential 
equations. 

^ EacH of the systems (11), (12) presents u as satlsMng 
simultaneouslj two linear partial differential equations of the 
first order. 

To deduce the value of u thus conditioned it will obviously 
suffice to multiply in each system one of the partial differen- 
tial equations by an indeterminate multiplier X, to add the 
result to the other equation so as to form a. new equation 
which will, like those from which it is formed, be linear and 
of the first order, and which on account of the indeterminate 
character of \ will be equivalent to the two. From the 
auxiliary equations which we obtain in the process of solu- 
tion, \ must be eliminated. 

If in this way we combine the equations of the system (11), 
we have, on arranging the resulting equation according to the 
differential coefficients of w, 



Rp-\-m^ + \{Tq +»»,;>)!• ^ 



dp dq 

Hence we have the auxiliary equations 

dx ^ dy dp dq dz 

5 + Xm,~Wj + \r~T'""xF""5p + w,5r + X(rj + w,p)' 

du = 0, 

and it is to be remembered that m„ m^ are the roots of the 
equation 
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Eliminating X from the first four of the above equations we 
have 

Udgi + m^dx — Rdy = 0] 

Udp -{- m^dy - Tdx = o\ (I). 

dz — pdx — qdy = ol 

This then is the system of ordinary differential equations 
deduced from (11) upon the integration of which the determi- 
nation of u will depend. 

A similar system, differing from the above only in th^ 
mutual transposition of m^ and w,, is given by (12), viz. 

Udq + m^dx — Rdy = 0] 

Udp-\-m,dy--Tdx^(yt (II). 

dz — pdx — qdy = 0) 

If from either of these systems we can deduce two inte- 
grals of the forms 

it is obvious, from what precedes, that 

u =f{v) 

will constitute a first integral of the proposed (1), and there 
being two systems in question, two such first integrals, each 
involving an arbitrary constant may coexist. 

6. Pbop. IV. To deduce ike second integral of (1). 

It will be necessary to consider separately the cases in 
which m^ and m^ are equal and unequal. 

First let m^ and m^ be equal. 

Both the systems (I), (II) reduce to a single system which 
may be expressed in the form 



djp = -^dx - 'jjdy 
dz = pdx -f qdy 



\ (14). 
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Now, since the condition -f = ^ ia here satisfied, it is 

ay dx 

manifest that if from any two mtegrals of the above system 
of the forms u = a, v = h, simultaneous values of p and a T)e 
determined, these values will render the third equation otthe 
system integrable, and the eflfect of its integration will be 
virtually to determine » as a function of a?, v, and three arbi- 
trary constants, viz. a, b, and a constant c introduced in the 
last integration. Let us represent the result in the form 

z = (f>{x,y,a,b,c) (15). 

Now what relation will this result bear to the general solu- 
tion of the partial difierential equation given, to the solution 
which we should obtain by integrating, not the particular 
eauations w = a, v = J, but the general first integral u =/(t?), 
which includes them both. 

To integrate the equation u =/(v) it suffices to deduce any 
particular equation involving an arbitrary constant h, which, 
m conjunction with u =f{v) will render 

dz —jpdx — qdy = 

integrable, and to integrate the last equation regarding the 
arbitrary constant of integration as an arbitrary function of 
b. The result is a complete primitive in which, by the 
variation of i as a parameter the general integral is implicidy 
involved. 

Now either of the equations w = a, v = b will, in conjunc- 
tion with u =/(c) determine 2> and j so as to make 

dz —pdx — qdy = 
integrable. Take the equation v = J, then u =y*(t7) reduces to 

Thus, in place of the equations u = a, t? = J, of the previous 
section, we have 
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for the determination of p and jr« The constant c introduced 
in the final integration becomes also, according to the above 
theory, a function of J, and the complete primitive is of the 
form 

« = 0{a?,y,J,/(J),^(J)l (16), 

while the general integral is found by eliminating J between 
this equation and its differential with respect to h. 

The general integral therefore represents the envelope of 
the surmce represented by (15), a, J, c being parameters sub- 
ject to any two connecting conditions. 

As w, , m^ are supposed equal, a necessary condition of the 
possibility of this species of integration is that 

ig«-4(5r-CrF)=0 (17), 

the value of m is — , and the system (14) reduces to 



dz — pdx — qdy = 



(18). 



We conclude therefore that if (17) he satisfied and we can 
from (18) deduce a value of z in terms of x, y, and three 
arbitrary constants, the equation expressing that value will be 
a complete primitive, ana the general integral will be found 
by making the constants vary in subjection to two arbitrary 
conditions. 



Ex. Let the given equation be 

xqr + ypt + xy{s^— rt) =^pq» 

Here R=^xq, /S = 0, T—yp. U—xyy V^pq. 
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The condition (17) is satisfied, and (18) becomes 

xydp —ypdx = 0, 

xydq^ — xqdy = 0, 

dz —pdx — qdy = 0. 

From the two first of these we find 

p = ax, q = hy, 

whence from the third, 

cux? ht^ 

z = h —- + c. 

2 2 

This is the complete primitive, and the general primitive 
consists of all possible equations derived firom this by making 
a, J, c vary in subjection to two conditions. 

Ex. 2. Given 

(14-g^)r-2pg/+(l4-p')^-- /7^^ . 

=-(i+/+2')*. 

Here the equation for m reduces to 

IV? + 2p2'ffi +^V = 0, 

whence m = —^2, and the system (18) gives 

— -^— rj + PJC&! + (1 + 2«) % = 0, 
(1 ^jp' + grr 

^ +Mdy + (1 +^') dx^O. 



(1+/+2*)* 
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Subtracting the upper equation multiplied by pq from the 
ver one multiplied by 1 + ^^, and dividing by l+i?*+2'> 



lower 
we have 



(1+/ + 2')* 



whence 



"' + V(l+/+2')~"' 



In like manner. 



2' + V(H-i>*+2')^*' 
Hence determining p and j, 

''''' vii-(^-«r-(y-5)r 

Therefore {x - a)' + (y - bf + (is - c)' = 1. 

From this form of the complete primitive it is evident that 
the general integral will represent all possible tubular surfaces 
formed by the motion through space of a sphere of constant 
radius unity. 

Secondly, let tw^ and m^ be unequal. 

Then since, in neither of the systems (I) and (II) is the 

condition ;t- = ;>■ satisfied, from neither system separately 

can values oip and q be obtained which make dz ^pdx + qdy 
integrable. 

But, as will be shewn, any two integrals obtained, the 
one from the one system and the other from the other, will 
give values of p and q which will render dz = pdx + qdi/ iu- 

B.D.E. II. 11 
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tegrable, and the general solution will consist of all possible 
integrals of the latter equation thus obtained. 

Or if the complete first integral of either system be com- 
bined with any particular integral involving an arbitrary 
constant obtained from the other, the two will furnish values 
of 2> and q which render dz=pdx'\-qdy integrable, and its 
integral will be a complete primitive involving one arbitrary 
function in its expressed form, another in the connexion of 
its two constants ; the general primitive being found in the 
usual way by making the constants vary as parameters in 
subjection to a single arbitrary connecting condition. 

In fact it may be shewn that if we attempt by the process 
of Charpit or Lagrange to integrate the partial difi*erential 
equation of the first order u ^f{v), deduced we will here sup- 
pose from the system (i), we virtually construct the system 
(ii) in the auxiliary equations upon which the process of 
solution turns. I have obtained a direct proof of this proposi- 
tion, but I think it preferable and at the same time sufficient, 
to direct attention to the prior ground upon which it rests in 
the relations of the systems of partial difierential equations 
(11), (12) from which the systems of ordinary differential 
equations (i), (ii) are derived. 

Let P=0 represent any integral of the system (11), and 
Q = any integral of the system (12). Then we have 



dp 



dq ^\dx) \dy) 
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Hence we deduce 

\dp \dxj dp \dx) dq^ \dy) dq \dy J) 

={«©-.(f)}(iV{^©-.(f)}© 

The second member of this equation is identically 0. Hence 
dividing by F we have 

dp\dx) dp\dxj dq\dyj dq \dt/ J ^ ^' 

But this is, the known condition under which the values 
of p and q deduced from the equations P = 0, Q = make 
dz=pdx + qdi/mtGgY8i}AQ; see Chap. xiv. Art. 13, Equation 

(36). 



We conclude then that if from the systems (l), (li) we can 
deduce two corresponding systems of integrals 

u^=a^y ^i = ^j 
then will the first integrals of (1) he 

while the second integral will consist of all possible relations oi- 
tained either 1st by specifying the forms off , f and obtaining p 
and q a>s functions of x and y and integrating dz =pdx + qdy, 
or 2ndly, by specifying one of the functions f, f, leaving the 
other arbitrary, determining p, q, integrating dz = pdx + qdy^ 
and regarding the final constant of integration as an arbitrary 
parameter. 

11—2 



^ — ady = 0| 
\^y — cdx = 0) * 

^x — ac?y = 0) 
\^y — cdx = J 
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Ex. Given ar + bs + ct + e (s* — rt) = h, the coefficients 
being <5onstant. 

Here JS=a, S=h, T=c, Z7=e, F=A. 

Hence m^, m^ are the roots of 

m'^ — Jw + ac — eA = 0, 
and the systems (i), (ii) give 

edq^ + m^ — a^y = 0] 

edq^ + m^dx — ac?y = 0] 
edp + Wj 

Whence the first integrals are 

ej + wijOJ - a2^=/ (ei?+ Wj^y -car), 
e^A-m^x - ay =^ (ep + m^y - ca?), 

from which all possible second integrals are to be derived in 
the modes above explained. 

Let us take the second of those modes anfl give to the 
second of the above first integrals the particular form 

ep + m^y — cx^ Of 

C being an arbitrary constant. From this, and from the 
otlier integral, left in its complete form, we have 

cx-m^y+ ay- m^x+f^ {(m, - m,) y + 0} 

whence, substituting in the formula dz=:pdx + qdy, inte- 
grating, replacing the arbitrary form 



\Mt)dthj {m,^m;)<l>(t)f 
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and introducing an arbitraiy function of (7 for the arbitrary- 
constant, we have 

for a complete primitive. The general primitive consists of 
all possible relations obtained by eliminating G between the 
above equation and 

0=i[a: + f{K-m.)y+(7}] + t'(C), 

when the forms of ^ and -^ are specified. 



Second Investigation. 

G. If from the equation 

^r + iSf5+2Tf+Z7(5'-r«) = F. (20), 

we eliminate r and t by means of the equations 

dp = rdx + 8dy^ 

dq = sdx + tdy^ 
the result will be 

{Rdy' - Sdxdy + Tda? - U {dpdx + dq^dy)'] s 

=^Rdpdy+ Tdqdx- Udpdq- Vdxdy (21). 

There are diflferent considerations (all of them however in- 
volving, as I have been led to think, a more or less explicit 
reference to some theory of the genesis of the given partial 
QifTerential equation) which indicate that its solution depends 
Upon that of the equations obtained by equating to the part 
affected and the part not affected by «, viz. upon the solution 
of the equations 

Rd/ " Sdxdy ■{- Tdj? - U{dpdx + dqdy) = (22), 

Rdpiy+ Tdjdx- Udpdq - Vdxdy = (23). 
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Without entering into these considerations let us inquire 
what consequences may be deduced from these equations 
assumed to be true. 

It is seen that these equations are connected by a remark- 
able reciprocity with the partial diflferential equations (7). 
They will in fact be converted into these equations if we 
change 

dx, dy, dp, dq, U,V, 8 (24), 

into 

dv, du du du jr -rj Q /QrN 

'^dy^'^Tx' 1^' Tp' ^' ^'"'^ ^^^^ 

respectively. From this formal connexion it follows that if 
we multiply (22) by \ and add to the result (23), we shall be 
able to determine \ so as to permit the resolution of the equa- 
tion thus formed into linear factors. Ultimately we shall, as 
appears from Art. 3, reduce the system (22), (23) to an equi- 
valent system of the form 

{^-Rdy-m^dx-^t TJdq) {-m^dy^ Tdx 4- Z7^) =0, 
{" Hdy -- m^dx + Udq) {-m^dy- Tdx+ Udp) = 0, 

w, and m^ being determined by the equation 

or, changing the sign of m, 

(- Edy + m^dx 4- Udq) {m^dy - Tdx + Udp) = 
(- Edy -i- m^dx 4- Udq) [m^dy - Tdx + Udp) 

m^ and m^ being as in the former investigation roots of 



■"^L.(26), 

= oJ 



Equating to the corresponding factors of the first mem- 
bers we have 



— Rdy -{-mydx + Udq = 

— Edy + m^dx + Udq 

m^dy - Tdx + Udp = 
m^dy — Tdx -i- Udp 



= 0) 

= oJ 
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The first of these, m^ and m^ being diflferent, is resolvable 
into 

Udq - Rdy = 0, (& = ; 
the second into 

and it is obvious that neither of these can lead to the given 
partial diflferential equation (1). The first of them combined 
with the equations 

« 

djp^rdx+sdy, dq = sdx + tdy (27), 

leads in fact to the partial difierential equation 

JJ- Ut = Q (28), 

the second in like manner leads to 

T- Ur=0 (29). 

But equating to the non-corresponding factors of the first 
members of (26) we have 



::}. 



— Rdy + m^dx + Udq = Q 
m^dy — Tdx + TJdp 

m^dy - Tdx^ Udp = 

— Rdy + m^dx + TJdg^ = 0, 



Now these systems when completed by the equation 
dz = pdx + qdy agree with the systems (l), (ii) deduced in the 
previous investigation. 

It remains to shew that these systems actually lead to the 
given partial dififerential equation (1) directly. Eliminating 
from either of them, combined with the system (27) the dif- 
ferentials dxy dy, dp^ dqy we shall have as the result 

U{Rr + 88+Tt+U{s^-'rt)-'r} = (30), 

which, rejecting the factor ?7, as from (13) we rejected F, is 
the differential equation proposed. 
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Ground of the Reciprocity above noticed* 

7. The reciprocity above noticed is not of a primary cha- 
racter, but is founded upon two prior laws which I shall pro- 
ceed to demonstrate. 

If from the partial differential equations of the system (7) 
we eliminate Fand substitute the resulting equation in the 
place of the first equation of the system we shall obtain the 
equivalent system 



\dqj dq dp \dpj 



^^|(s)r(l)ll 



= 



^(l)'+«(l)(l)+^(l)' 



\\axj dp \dyj dq^ 



= 



-...(31). 



J 



These equations are both symmetrical and it will be ob- 
served that they are convertible the one into the other by 
changing 



du du du du jj ^ ^ 
dq' Tp' ^' dy' ^' ^' ^ 



(32), 



into 



du du du du 
dx' dy' dq' dp 



, V,-8, U 



(33), 



respectively. This is a law of reciprocity which connects 
solely the differential coefficients of u and the coefficients 
U^ 8, F of the original equation. 

Again w =0 is by hypothesis a solution of the given partial 
differential equation. Regarding it however simply as an 
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equation which is trvs and the truth of which is consistent 
with that of the equations 

dp = rdx + sdy\ 

dq = sdx + tdy) 

and differentiating it first with respect to a?, secondly with 
respect to y, we have 

du du d'z du dp du dq __ 
dx dz dx dp dx dq dx ' 

du du dz du dp du dq_^ 
dy dz dy dp dy dq dy ' 

equations to which we may give the form 

fdu\ __ du du 
\dxj dp dq 

du\ du . du 
s -T- +t -y-. 



_fdu\_ 



dyj dp dq 

Now this system is of the same form as the system (27) 
and will agree with it if we change 

_^du du du du , . 

"dx'^Ty' dp' Tq ^ ^' 

into 

dp, dq, dx, dy (36), 

respectively — a change which does not affect the coefficients 
of the given equation, and which is therefore the expression 
of a law of reciprocity distinct from that last noted. The 
combination of these two laws does however lead to the 
law exemplified in the researches of the previous Article ; 
see (24) and (25). 

The question here arises whether it would not have been 
better to employ from the first the symmetrical forms (31) of 
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the partial differential equations of the first order and second 
degree upon which u depends, than the unsymlnetrical forms 
(7). It was indeed from the symmetrical forms that the 
chief results of this paper were originally obtained, but the 
unsymmetrical forms lead to the same end in a simpler way, 
and therefore they have been made use of in the present 
memoir. 

It may be proper to notice, in concluding this sectiort, that 
the symmetrical forms in ordinary differentials would have 
emerged in place of the unsymmetrical ones of (22) and (23), 
if the quantity 5* — rt had been retained instead of 5, The 
equations 

dp = rdx + sdy^ dg^^sdx + tdy^ 

« 

enable us in fact to reduce the given equation (20) to the 
form 

Edp" + Sdpdg^ + Tdg" - V (dpdx + dqdy) 

= {s^ - rt) {Rdf - Sdxdy 4 Td^ - U {dpdx -f dqdy) }. 

Hence arises the symmetrical system 

Mp^'+Sdpdq + Td^ - V{dpdx + dqdy) = 0, 

JSdy'' - Sdxdy + Tda? - U {dpdx + dqdy) = 0, 

which is connected with the system (31) by the single law of 
reciprocity expressed in (35) and (36), 



Postscript. 

8. At the time when the above investigations engaged my 
attention I was totally unaware that the subject of them had 
been discussed by Ampfere {Journal de VEcoU PolytechniquBy 
Tom. XI.) and recently by Professor De Morgan {Cambridge 
Philosophical Transactions, Vol. IX. Pt. iv.). I feel it there- 
fore incumbent upon me to state why after acquainting my- 



OP THE SECOND OEDER. 



self with the results of theii" laboura, I offer this paper for 
publication. 

The method of Professor de Morgan so far resemtles tlie 
first method of this paper, and that of Ampere the second^ 
tiiat -while the former makes the solution of the problem 
md directly upon that of simultaneous partial differential 
equations of the first order, the latter makes it to depend 
■directly upon the solution of simultaneous ordinary differ- 
'Cntial equations of the first order. The formal connexion of 
these methods by the law of reciprocity is, I believe, esta- 
blished for the first time iu this paper. The system of partial 
differential equations of the second degree (7) has not, so far 
as I am aware, been given before. 

But a point which I think of deep importance is. the follow- 
ing. By connecting, aa in this paper, the differential equa-' 
tiona of the second degree, whether ordinary or partial, by 
indefinite multiplier which is afterwards determmed so as 
admit of the resolution of the system into its component linear 
elements, we aasuro ourselves that each step ot^ the solution 
'Offers a complete sequence to that which has gone before, and 
it only remains then to separate the different elements and 
determine whether they are relevant or irrelevant to the end 
HI view. That any such distinction exists has not, so far an 
I aware, been noticed before. And it seems to me tha' 
more important that it should be noticed because the solution.' 
of partial differential equations in cases far more general than 
those above considered seems to depend upon the satisfaction 
of simultaneous differential equations of a degree higher than 
the first. I have in fact by an application of the Calculus of 
Variations arrived at the conclusion that the theory of the 
solution of all partial differential equations of the second 
order, whatever the number of variables may be, is very inti- 
mately connected with the satisfaction of a system of differ- 
ential equations of the type 



1 

d. 



4 



dF 



dy* 



dF, 



^f" 



'= representing the given partial differential equation, ; 
id y (my two of the independent variables, and r 
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second differential coefScients of the dependent variable with 
respect to x and y. 

I may perhaps at some future day resume the subject, to- 
gether with an inquiry into the theory of the solution of the 
partial differential equation of this paper, when the conditions 
under which the auxiliary equations (i), (ii) are supposed to 
be integrable are not satisfied. 

9. Note. It may be desirable to establish directly the 
converse form of one of the results of Proposition iv. For 
this object we shall shew that the equation of the envelope of 

z = (l){x, y, a, J, c) (1), 

where a, b, c are connected by any two conditions of the 
forms 

i|r (a, J, c) = 0, x(a, J, c) = 0, 

will satisfy a partial differential equation of the form 

Br+Ss+Tt+U{s^-rt)=r (2), 

in which also 

m 

Differentiating (1) we have 

d(b dd} da dS db dS dc 
^ dx , da dx db dx dc dx 

dd) dd) da dd) db dd) dc 
^ dy da dy db dy dc dy 

and by the nature of an envelope these reduce to 

d^ dp 
^^ dx' ^- dy '^^)- 

Again differentiating these equations with respect to a; and 
y, and writing for simplicity 

-^=^ i!l=5 ?-^-=c 

dadx ' dbdx ' dcdx ' 

^'<^ __A, d^<i> _ pf d^<f ) __ ^, 
dady ' dbdy'~ ' dcdy ' 
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we have 

— ^^ J. J ^^ j_ 7? ^^ n^ 
da? dx dx da? ' 

dxdy dy dy dy ' 

dxdy dx dx dx^ 

d^<}> J/ ^« , TV rfi . ^t dc 
dy dy dy dy 

Hence we find 

\ dxdy) V dxVK dy") 

\ dy dy dyj\ dx dx dx 

\ dx dx dxj \ dy dy dy 

^ ' \dy dx dx dyj 

db dc dh dc 



{BC'-B'C)(f^-f^) 

\dy dx dx dy) 



' \dy dx dx dyj 



Now since a, 5, c are connected by two conditions, so that 
h and c are functions of x and y onlj as being functions of a, 
we have 

da db da db __ db dc db dc _ 

dy dx dx dy~ ' dy dx dx dy ' 

dc da ^^ ^ ^r\ 
dy dx dx dy 
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Thus the above equation reduces to 



or 



dy^ dxdy da? d^ dy^ \dxdy/ 

This equation is of the general form (2). Its coeflScients 

-j^ , &c. are determinable as fanctions of a;, y, z^ p, q when 

the form of the complete primitive (1) is given. For this 
purpose the complete primitive with the two derived equations 
(3) suffice. 

Again, comparing (4) with (2) we have as the conditions 
of their equivalence 

B 8 T jr y 

dy^ dxdy da? da? dy* \dxdy/ 

conditions which suppose -B, 8^ jT, ?7, V connected by the 
relation 
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CHAPTER XXX. 



ADDITIONS TO CHAPTER XVII. 



[The present Chapter consists of additions to Chapter xvii. 
Art. 1 was intended to follow Chap. xvii. Art. 1.] 

1. The theory of the solution of linear differential equa- 
tions in a series flows very beautifully from their symbolical 
expression. It is usual in treating this subject to assume the 
form of the series, and deduce from the differential equation 
the law of its coefficients ; but the symbolical form of the dif- 
ferential equation determines in reality the form of the solu- 
tion as well as the law of derivation of its successive terms. 

Let us begin with the binomial equation 
Operating on both sides with {/o(i>)}~^ we have 

in which *^^)=iS)- 

Hence {1 - <^ {D) e^]u = {/, (i>)}"'0. 

Now {^(-D)]"^0 will be determined by the solution of a 
linear differential equation with constant coefficients, and will 
be necessarily of the form 

^P+jB^+CB+..., 

in which A^ B, (7, ... are arbitrary constants, and P, Q^R^ ... 
are functions of the independent variable. 
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We have then 

{l-<l){B)e'0]u = AP+BQ+OB+..., 
therefore m = {1 - ^ (Z>) e'"}"' (AP. + BQ + GB +...). 

Now let us represent ^{B) e*"* by p ; then 
u = (1 -p)"' {AP+BQ + CB + ...) 
= {l + p + p^ + p'+ ...) {AP+BQ+CB + ...) 
= A{l + p + p'' + p'' + ...)P 
+ B{l+p + p'' + p'+...)Q 
+ C(l + p + p'' + p' + ...)B 
+ ... 

Eepresent the first line of the above expression by u„ 
then since 

p" = ^ (D) ^mp) e'* . . . j» times 

= €'^'*^{B + mr)j>{JD + mr-r) <i>{B + r), 

we have 

M. = ^ {P+ e*-*^ (2? + r) P+ €*-*^ {D + 2r) </> (i) + r) P 
+ 6"-»<^ (2? + 3r) ^ (2) + 2r) <^ (I> + r) P+ . . .}, 

In which it only remains to perform the operations indicated 
by <^(i> + r), by <^ (D + 2r) ^ (D + r), . . . on the function P. 

Let ns in the first place suppose the symbolic function 
f^ {D) to be of the form (D — a) (x) — i) ... ; then 

Here P= e"^. Hence substituting in the above expression 
for M, and observing that/(i>) €^=f{n) €^^, we find 

or, since €* = a?, 

r 
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and 

tt = u4a;* {I + ^(a + r) «•■+ ^(a +2r) ^(a + r) aj**+ ...} 
+ Bx"" [I + ^{h + r) of + ^(5 + 2r) ^(5 + r) 0;''+ ...} 

the solution sought. 

Consider now the general equation 

UD)u^f,{D)^u+...+f,{D)^u^O. 

Here we have, representing Yf-A hy ^mifi), 

{l + ^,(2>)6«+.... + ^„(i))€'^}t. = {/,(2>)}-^0; 
therefore 

« = {1 + 4>,{D) e» + ui>)^v {/o(^)r 0. 

Here we have first to determine {^(i^)}"^0, then to deter- 
mine the eflfect of the operation represented by 

upon this. 

Now {fQ{D)Y^O is given by the solution of a linear diffe- 
rential equation with constant coefficients, and will therefore 
be of the form 

A J By (7, ... being arbitrary constants, and P, ^, JB, ... func- 
tions of 6. 

Again, since 

{l + ^,(i>)e« + + ^.(2>)6"r 

^ 1 

l + <^,(/>)6^+ + 4>,{D)e'' 

it may be shewn by a process of actual symbolical division, 

B.D.E. II. 12 
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attending to the laws of combination of symbols, that the 
expression may be expanded in the form 

To determine the functions F^{D), F^{D), we may 

proceed as follows. From the equation 



1-1 



{1+^.(2)) 6»+ +MD)^y 

we have 

= F^{D) + {F,{D)+<I>,{D)F,{D^1)}^-]- (1). 

Hence j;(2>) = l, 

F,{D) + <l>,{D)F,{D-l)^0; 
therefore F, (2>) = - ^, {JD) F, (i> - 1) , 

and so on. Hence F^ (D) , F^ {D) , are determined in suc- 
cession. The general law is as follows: the coefficient of 
€"*^ in the second member of (1), when m is greater than 1, is 

F„(D) + il>,{D)F^,iD-l)+<f>,(D)F^{D-2) + ... (2), 

whence 

F^{I)) = -<j>,(I))F^,{D^l)-,l>,{D)F^{D-i)-.... 

By this formula the successive values of F^{D) can be 
deduced from those of F^^ [D) , F^^ (-D) , . . . . 

Combining the above results we obtain thus for u the ex- 
pression 

u=^{l'^F,{D)^ + F,{D)^^-....}{AP^-BQ + ...} 

= A{P-]-F,{D)€'P+F^{D)e^P+....} 

+ B{Q + F,{D)e'Q + F,{I))e^Q+....} 

Let us in applying this expression first suppose that tbe 
factors of f^ [D) are real and unequal, so that j^ {D) is of the 



ART. 1.] ADDITIONS TO CHAPTEE XYll. 179 

form (D — a) (2> — 5) (D — c) . . . • Further, let ns suppose that 
no two of the quantities a, 5, c, .... difier by an integer. 

Then {f,{D)}-^O^A^ + B^'+...., 

whence we may assume 

Thus the expression for u becomes 

+ B{^ + F^{D) €^')' + F^{D) €<^>«+ ....} 

or, since F {D) ^=F (m) e"*, 

tt =A{^+ F^{a + 1) c"^'«'+i?;(a + 2) €<•*"•+ ....} 
+ B{i» + J; (i + 1) €'"•'« + i^, (J + 2) e'"*!* +....} 

Hence, replacing e* by x, 

tt = ^ {x- + i?; (a + 1) a;*" + i?; (a + 2) a!«*» + , . . .} 
+ B{ai' + F^{h + 1) ar' + F,{b + 2)^^ +....] 

In (2) replace in like manner D by a + 1 and we have, 
putting i for m^ 

Fi {a + i) + <}>^{a + t)Fi_Xa + i''l) 

+ ^,(a + i)i^i^(a+i-2) + = 0, 

or, if i^(a + i) be represented by Wa^.<, 

Put m for a + 1, thus 

This agrees with the law established in [there Is no refer- 
ence in the manuscript, but the law intended appears to be 
that given in Chap. xvii. Art. 9.] 

12—2 
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Secondly, suppose that r of the factors of ^(-D) are equal 
and of the form D — a. 

Then {j^(-D)}"^0 contains a term of the fonn 

Hence the corresponding portion of u is of the form 

|l + i^,(D)€* + i^,(D)6»^+'...J6^(c, + c,5+ + c,.,n 

= |e«^ + e(*^^)^i^^(2) + a + l) + €^^i^,(2) + a + 2)4-...U...(3) 

where v stands for c^ + cfi + c^SP +.... + c^^^"*. 
Now F^(P'\-a + %)v 

= k(a + l) + i^/(a + l)l) + i^/'(a + l)j^+....U, 

which on performing the differentiations becomes a polyno- 
mial of the form 

A+A^ + +A^^sr\ 

We see thus that (3) will assume the form of a series of temw 
€"^, 6^""^*^^ .... each multiplied by a polynomial of the (r — 1)* 
degree in 0. Or arranging the terms otherwise it will con- 
sist of a series of terms of the form 

B, + B,0+ +5,.,^-^, 

in which 5., B^, J5^_j are series involving €^, c^, 

€<"^*, Or lastly, changing e* to aj, the portion of » in 

question is of the form 

B,+B,{iogx) + +B^, aoga?r\ 

-Bq, J?j , . . .. 5_j being polynomials in each of which the lowest 
power of X is a;*, and the successive powers increase by unity. 

This establishes the assumption in [there is no reference 
in the manuscript; probably Chap. xvii. Art, 10 is to be 
supplied.] 
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Thirdly, let f^{D) contain r factors D — a^y D — a^^.^.D-^cL^ 
in which a^ya^y...ar dififer from each other by integers, toge- 
ther with other factors. 

The portion of u corresponding to the factor D-^a^ will be 

in which 

F^ (i)) e-» = - 1^, iP) ^F^, {D) 6<-')« 

Thus F^{p) ^ consists of terms of the form 

i being one of the numbers 1, 2, ..*. w. Hence F,^{D)e^'^^ 
will consist of terms of the form 

4>,{P)^F^{D)e'^^», 

j being one of the numbers 1, 2, ... w. Continuing this until 
t +y + &+ ... = w, we see that F^{D) ^ will ultimately con- 
sist of terms of the form 

*»y» ^> ••• receiving arbitrarilv any of the values 1, 2, ... n, 
and t+j+k + ... oeing equal to m. 

Thus the portion of « derived from A^ will consist of 
all possible terms of the form 

A<l,t{D) ^<I>,{D) e*^(Z?) €*» .... {^ 

^ AMD)MD-t)MD-{-j)... ,^, 
f>mf,{D-xUiP-i-J)'" ' 
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Let ^ = a, i-hj = )8, i +J + k+ ..... excluding the last 

term = fi ; and let the symbolical numerator which involTes 
only direct functions be represented by/(i)), and we hyaye 

^fiP) p(«+-«i)» * 

/,(^)/o(i>-«)/,(^-/3)..-/o(^-/*) ' 

in which a^ ff, ....fi,m are integers ascending by differences 
not exceeding w. 

[A few lines of the manuscript here are obscure, and I 
venture to express in other words the idea which seems to be 
involved. 

Let D — a, denote one factor o{f^{D), then the correspond- 
ing factors in the denominator of 

JfiD) eC'^d' 



/.(^)/.(-D-a)/.(i>-^).../,(^-/*) 

are {D — a^) (i>— a,— a) {D^a,-^ fi) ..........(6). 

Now if a, is not greater than a,, then a^ + fi is less than 
a^+m; hence no factor in the expression {5) can be identical 
with D^m — a^. But if a, is greater than a,, then one fac- 
tor in the expression (5) may be identical with J9 — wi — a^. 

Hence it follows that the denominator of the expression (4) 
may contain D^m^a^ to the power r — 1, but not to a 
higher power.] 

And, since 

= «(•»*•.>• |c. + c,d + . . . + cr^e^ + j^} , 

we see that u will contain r sets of terms together of the form 

^ + 5(logaj) + (7Gogaj)'+ + K(^ozxY'\ 

A, Bj C... being polynomials in x. 

This establishes the rule in [there is no reference in the 
manuscript; probably Chap. xyii. Art. 10 is to be supplied.] 



ABT. 2.] ADDITIONS TO CHAPTER XVII. 183 

[There is no hint in the manuscript as to the position 
■which Article 2 was intended to occupy ; and the reasoning 
does not seem fully developed.] 

2. Peop. The solution of the equation 

f>{D)u+f,{D) 6»u+ +/.(Z?) e-«t = 

"being expressed in the form 

it is not necessary to introduce new constants in interpreting 
JPj(i>), .... ; it suffices to interpret particularly if only uni- 
formly and consistently. 

For let 

{f.{D)ro^AF^BQ+ ; 

and in interpreting 

F^{D)^{AP+BQ + ) 

let a new constant be introduced which was not in the inter- 
pretation of 

F^,{D) 6^*"-^)^ {AP+ BQ + ). 

Now F^{D) eJ + <l>,{D) ^F^,{D) e^"-^)^ 

+ <I>,{D)^F^{D)^'^' + =0, 



therefore 

F^{D) €«^ = - <^,(D) e'F^,{D) e^-^)*- .... 

hence the new constant comes from {/^ (-D)}"*0, and the term 
containing it must be A'P, or JS'Q, ...., where -4', JB* .... are 
constants. Suppose it A'P; 

then as derived from this, 
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The law of derivation is exactly the same aa in the deii- 
vation of JJ [D) e« from i^j., {D) ^*''i^, 

[Art. 3 seems intended for a reconstruction on an extended 
scale of part of Chapter XYir. Art. 3.] 

3. We proceed to consider more fully the theory of the 
binomial equation 

Now the poasihility of solving tbe equation depends upon 
the nature of the symbolic function <j>{-D). It is perhaps the 
most general account of the present state of the theory to say 
that there exist certain primary forms of this function which 
render the equation solvable, and that to each of these pri- 
mary forms an infinite number of the forms are reducible by 
general theorems of transformation. As these theorems 
admit of a statement which is independent of the form of the 
function </i(-D), we shall establish them first. 

Prop. II. The function ^(i?) in tAe equation 

can without otherwise changing the first member of that eola- 
tion he ist affected with any constant factor, or 2ndly con- 
verted iiito (f){I> + a), or 3rdly converted into {(f>{~D)}''. 

First. Let U=f{e^, and in the equation 

let €* = o'e*'. Then '^ = ~Ja'i ^^^ ^^^ equation becomes 
w + o0(i))e'*«=/(aV), 
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in which -D = -^ . Thus ^(-D) has been affected hy a con- 
stant factor a. 

Secondly. In the same system let u = ^v. Then 

or ^v + ^<f> {D + a) ^v =/(€^), 

therefore v + ^ (D + a) e^^i? = e'^^fi^ . 

Here ^{D) has been changed into ^(D + a). 

The result of this transformation may be conveniently ex- 
pressed by the following theorem. 

The equation 

wiU he converted into 

v + <l>{D + a)^v=V 

hy the relations 

Thirdly. In the same equation let 5 = — 5' ; then 

d___d^ 
dd^ dff' 

and we have 

in which ^^Tfff* Hence 

u + e-^^{- (Z) - r)) u =/(e-^) ; 
therefore e^u + ^ (r — D) m = e^'/(€~*^, 

whence w + {^(r-2>)r€^w = {^(r«i))r*€*^/(€-»). 
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In this equation let u = ^v. Then by the last theorem, 

« + {^ (- D)r^v = {^ (- 2>)r/(6-^). 

Thus 

is converted into 

v + {4>{- i))r e^t; = {j> (- 2))r/(€-^), 

in which D = -^ , Jy assuming 

The above transformations leave the index r in the first 

ff 

member unchanged. If however we assume ^ = — , whence 

-3^ = a -^ , we should have 

By combining this with the previous results we see 
that it is possible to convert 61B) into 6(aD + b)f and into 
{<l>{aD + b)r. 

But the most important transformation of the function <f>{D) 
is that which is established in the following proposition. 

[The proposition referred to is Prop. iii. of Chap. xvii. 
Art. 3.] 

[Article 4 was intended to follow the words "or subse- 
quently in the derivation of t^" in Chap. xvii. Art. 4.] 

4. It becomes therefore important to establish rules for 
the treatment of the constants which in these different ways 
arise. 
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Now the entire process of solution consists of three stages, 
namely : 

Ist, the determination of V bj the equation 

2ncll7, the solation of the transfonned equation 
Srdlj, the determination of u by the relation 

Let us consider these separately, supposing tf) (D) to con- 
tain a single factor -j^ — j- which is made to disappear in the 

generation of '^(jD), so that a and h differ hy a multiple of r. 
Thus the given equation is of the form 

u-^i.{D)er»u=U (6). 

The transformed equation is of the form 

in which «=P,g±->, r=P,J±-> 

First, suppose a—b = nr, where n is positive. 
Thus 

tt=(2> + a) {D + a--r) ... {D + a-nr+r)v^ 

F=|(i) + a)(i> + a-r) .., {D + a-^nr + A ' U. 
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Hence 

+ (7,e-<^^)*+ + a.€-(-'^4....(7), 

where C^ is a particular value of 

|(D+a)...,(2) + a-wr + r)l U. 

The part containing the constants will consist of terms of 
the form 



= (D + a) (Z) + a - r) ... {D + a-nr + r) ll +1^ (2>) 



^B 



r(«-*1« 



+ f (D) e'^f (2>) €^ + ...I (7€-(*- 

= C{D + a) (2) + a - r) ... (Z)+ a - nr + r) L'^^-*^^ 

+ f (i>) e-t**-*^)^ + f (D) f (D - r) e-^--""-*-)^ + . . .}. 

Now all these terms vanish up to the one containing 
^-ia-nrye . thcrcforc we have to perform the operation 

(7(i? + a) (2> + a-r) ... (2> + a-wr + r) on 

where y = w — i — 1 ; that is, we have to perform the operation 
C{D + a) {D + a-r) ... (D + a-nr + r) on 
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Now ^{D)-^{D-r)...-^{D -jr) e"'""'* 

= T^ [nr ~ a) ^{nr — r — a) ...■<^inr —jr — a) e~(°-"is 

= ££-("■')•; 

therefore we obtain 

BCiD + a) ... (D+a-nr + r)L-i"''^-i-f{D)^e~i''-''^+.. 

= SG{D + a). . .{D + a- fir + r)k+f {0)^0 

Thus thia expression ia the same in form for all values of I 

Therefore all the terms containing an arbitrary constant € 

jn (7) are equivalent to only one term. 

Secondly, suppose a — h = — nr. 

Then u = i{D + b){D + h-r)...iD + h-nr + r)rv, 

V={D + h){D+h-r)...{D + h-nr + r)U. 

Here there are no constants in V. But u contains n arbitrary! 
constants not in v, and as there ia no subsequent process in.l 
the method for destroying these or reducing them to mutnalB 
dependence, it ia necessary that the relations connecting them'J 
should be sought by comparing the aolution with that girenl 
by the method of development in series. 

Note. It would be better to reduce (G) to the form 

llxsfore the demonstration. 

[Article 5 was intended to follow Chap. xvii. Art. 7. 

There is a memoir by Professor Boole on the subject of thian 
lArticle, entitled On the Differential Equations which deter 
Wmitie the form of the Jtoola of Alijebraic Equations, Tia 
■memoir occupies pages 733 — 755 of the Philosophical Trai 
■- -WM for 1864.] 
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5. If we agree to regard as primary tliose forma of tino- 
mial equations which are integrable but not through any 
reduction effected by the Propositions of Art. 3, and to which 
equations through the application of those propositions other 
equations are reducible and so made integrable, it becomes 
very important to enquire what these primary integrable 
forma are. It does not appear at present possible to give a 
general answer to this question, but so far as is known, such 
forms if belonging to differential equations of a degree higher 
than the first stand in a remarkable connexion with the 
theory of algebraical equations. By the study of this theory 
Mr Hariey was led to the conclusion that y defined as an 
implicit function of x by the algebraical equation 

3,"_„y+(„-l)a! = (8), 

n heing greater than 2, satisfies the binomial differential 

equation 

y 27(7rri) ... (lJ-n + -2) '*«'"' 

in which i^ = x. In this expression the factora of the nnme- 
rator are equidifferent, as of the denominator, their common 

difference being , but the equation is not resolvable by 

Propositions ii. and iii. into formSj the integrability which 
bad before been recognised. 

The above result first reached by induction was confirmed 
by Mr Cayley by the aid of Lagrange's theorem, 

To the form (8) all algebraic equations of tlie third, fourth, 
and fifth degrees ai'e known to be reducible. 

Mr Hariey has subsequently found that y considered as a 
function of x defined implicitly by the equation 

r-«!/'-'+(n-l)a: = 
satisfies the symbolical differential equation 

tT' [{n - 1) DT'y -[n-\){nD~n- 1) [nD - 2]"-^^ 
= [«-l]'"'e'' 
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the factorial notatioa according to which 

[w]* = 7n(ni — 1) (w — 2) ... (w — w + 1) 

being here adopted. 

These results are implicitly involved in a more general 
theorem which I shall nOw demonstrate. 

Theorem. If y^, y, y, are the n roots of the algebraic 

equation 

y-ay"-^ + l = 0, 

and if the m*** power of any one of these roots be represented 
by Uf and log a by d, then u as a function of satisfies the dif- 
ferential equation 

\ n n \n n J 

« m^ 



^u = 0. 



And the complete integral of the above differential equation 
will be 

u=c^,''+cj/,^ +0^:^. 

Let y* = z, then the given equation may he expressed in 
the form 

in which 5 = — 1. Hence, by Lagrange's theorem, 

m m n— 1 J m 



db 



+ 



r^lF)"l»1-*=-. 



the general term being 



rt^arp)!'^}. 
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which on effecting the differentiations and adopting the ho 
tonal notation becomes 



971 



m + (n — 1) r 



n 



-1 



n [rj ' 



r-l m-r 



(£ 



and this expression will he found to represent the first term 
as well as tne others of Lagrange's expansion provided that 
we interpret the form 

M»b7l, and|>rT)7i^. 
Further, the above general development includes the n 

m 

particular developments of u or y* arising from the giving 



to J* 
sents 
succession. 



its n particular algebraic values. In this way it repre- 
the m"* power of each of the n roots y^, y„ f....y» in 



Now representing the above general term by w^^*', we shall 
have 



m 



^r = 



m + (?^ — 1) r ^ ' 



n 



m 



^r-n — 



n[r7 
"m + (w — 1) r 



A 



r— 1 mr-r 

b' 



"Ir-ii-J wiK-r 



n 



— n 



w [r— w]*^ 



+1 



Therefore, after reduction and replacing 6 by — 1, 



'^^ -- 


wi + (w — 1) r 


\n n J 


Wr-« 


M" 





(9). 



It follows therefore that the complete series of which the 
general term is n^'' will if represented by u satisfy the diffe- 
rential equation 
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..-L^ — " \4 ^" " ^ "■.=o...(i). 

If we integrate the equation in a series (Chap. xvii. Art. 9), 
the initial terms of the value of u wiU be 

c,+ Cfl + c,a* + a^.o--', 

the succeeding terras teing formed from these "by the law 

(9). Hence, if the arbitrary constants (7,C, C,^ be so 

determined as to make the above initial terms agree with 
tlie first n terms of the Lagrangeau expansion in any of its 
particular forms, the succeeding terms will also agree, and 
the Lagrangean expansion will thus become a particular inte- 
gral of the equation (I). The aggregate of such particular 
integrals, each affected by an arbitrary constant, will therefore 
also be an integral of the differential equation, and will, iq 
fact, constitute its general integral, subject to exception only 
in the case in which for a particular value of tn the integraJa 
Ifi'' y"> yn"" cease to be independent. 

For instance, if m=— 1, and we reduce the equation to 
the form 

it is seen that except when w = 2, we have 

yr'+yr' +y„''=o. 

Here then the solution 

tt=6'^."+C^,'- + C^- (10) 

to be general for it becomes 

»-(o.- (7jj,f+ (0,- c.)s,-"-+(c-.- c.)r'->. 

id virtually involves but n — 1 arbitrary constanta. 
If, however, we give to the integral the form 



- csr+ oj,: ...+ c„j-„+ c."^^"^— 



B. D. E. II. 



13 
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the last term of which becomes a vanishing fraction when 
7?i = — 1, we find for the general value of u in this case 

+ ^n (yi"" logyi + y*"" logy, +y»"' logyj, 

and in this way we may proceed in failing cases generally. 

Lastly, it may be observed that in certain cases the differen- 
tial equation (1) admits of reduction to an order lower by 
unity than its own. And in particular this happens in the 
failing cases above noticed. Thus, if in (I) we make tw = — 1 
the equation will be expressible in the form 

whence, operating on both nMOnbers with (J) — n + 1)"S we 

have 

1 fj^— 1 1 '1«^^ 

\nr-' w - i ^-^ D -. i _ 1 €*^e^ = C€^^'>\ 
"- -^ n[^ n ^ J 

The general integral of this equation will be expressed by 
(10) provided that a proper relation be established between G 

and tne constants 0^, C^, C^. If we choose to determine 

C so as to give to the integral the particular form y""\ we shall 
find on substituting for u its Lagrangean development making 
w = — 1, J = — 1, and calculating the coefficient of a*"* or 
^n~i)$ i^ ^Y^Q £j.g^ member of the differential, 



^^[n-3] 



n 



Hence, if n be greater than 2, we have (7=0. It follows 
therefore that if n be greater than 2, the equation 

[D]^'u - - l"!—^ 2> - 1 - llT'e^u = (II), 
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in Which ^^a has for its general integral 

^ = c'lyr* + ^ay«"' • • • • + ^^tjr^^x^ 

y^yy^^ y«_i being any n — 1 roots of the equation 

It may be useful to notice the forms which the above 
results assume when 6 is changed into — 0^ and therefore D 
into — jD ; see Art. 3. 



It will be found that (I) becomes 

ri> - IT 



M — 



n 



L 71 n^ \n n 



') 



^« = (Ill), 



of which the integral is therefore 

« = Cjy," + C^," + C4^.", 

i/i,y,, y, being; the roots of the equation 

1 



r--y""'+i=o 



(11); 



and log a being denoted by 0; 

while as the equivalent of (II) we have 

[D - 2]"-' 



u — n 



n n 



-i€^t^ = 0. 



(IV), 



of which, supposing n greater than 2, the integral is 

jfify^j y».i being any n — 1 roots of the same algebraic 

equation. ^ 

13—2 
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Mr Harley's results may readily be deduced from the 
above. Thus it will be found that the equation (11) re- 
duces to 

, r - w< + (n - 1) a? = 
if we make 



R-l n—l 



y'=(n-l) -»-«, a = -i ^ 



Hence, making a?= e^ and representing ^, by i)', wfe have 
for the transformation of (IV) 

(n-l)^'^V 

C— C « 

n 



n 



w — 1 



t* = (n - 1) ~ e * U 






Substituting and multiplying the result by €«», we find 



'-(^r- 



w jy 2n — 1 



n— 1 n—l 



«-i 



W] 



s-r 



.6<'^»)<>^=0, 



which is Mr Harley's first equation. 



If in (I) and (III) we make 1 = a, whence ^na n-no, 

and at the same time change a into ah ", and y into yh *, we 
shall obtain the following somewhat more general statement 
of their united import. 
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u 



tnr ^ 



0, 



w — & 



JD - 1]' 



n- 1 
n 



i) + a-l 



fU^ 



r 



£-... 



^w = 0, 



arc both satisfied hy the general integral 

trAen y^, y,, .... y» are <7^e roots of the algebraic equation 

provided that for the first eqttation a = €*, and for the second 
a = €"*. 

• If n = 2, the above equations assume the forms 

f -(«-!)• 

6. [The two principal papers by Mr Harley on the dif- 
ferential equations exhibitea on page 190 are the following : 

(1) On the Theory of the Transcendental Solution of 
Algebraic Equations, Quarterly Journal of Mathematics, 
Vol. V. pages 337... 360. 

(2) On a certain class of Linear' Differential Equations. 
Manchester Memoirs. Third Series. Vol. II. pages 232... 245. 
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In a letter bearing date January 1S| 1864, Professor Boole 
pointed out to Mr Harlej that his second equation might also 
be deduced from the general theorem discussed in Art. 5. 
Emploj^ing the above notation the deduction may be pre- 
sented in the following form. 

The equation (11) will reduce to 

if we make 

2^ = (n-l)"*ar"^^, a = - (w - l)"aj* ; 

and for the transformation of (III) we have 

t«=(n-l)"»€"" w. 

These substitutions being effected we arrive, after some 
slight reductions, at the following equation, 

n•[(n-.l)2y-7^^]•-^JDV-.(n-l)[n2>'-.m-l]Vu' = 0, 

which, making wi = 1 and u = t, gives 

w* [{n -1)2)'-. 1]-^ m - (n - 1) [/lD' - 2] V< = 0, 

an equation which admits of reduction. In fact, operating on 
both members with (2^ — 1)'^, and determining the constant, 
as in the former case, by the aid of the Lagrangean expansion, 
we find 

n*-' [(n - 1) lyy-H ^ (»2)' - n - 1) [nl)' - 2]'^€^« = [w - 1]**^, 

which is Mr Harley^s second equation. 

The references and deduction here given were to have 
been added to the memoir which is cited in page 189, ac- 
cording to Professor Boole's desire; but by some accident 



ART. 6.] ADDITIONS TO CHAPTER XVII. 199 

they were not printed, and the omission was npt discovered 
until after his death. 

Mr Harley has lately succeeded in obtaining the fol- 
lowing extension of Professor Boole'3 theorem. 

The differential equation 

[_ n ax n j [n ax n J 
is satisfied by the m^^ power of any root of the equation 

y* — xy*"'' + a = 0, 
u being considered as a function of x. 

From this he deduces the following ; the differential equa- 
tion 



n 



L r dx r \ L ^^J 






is satisfied by the m^ power of any root of the equation 

y*- w^*"*" + (n- 1) aj = 0. 

For the materials of this Article I am indebted to Mr 
Harley.] 
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CHAPTER XXXI. 



THE JACOBIAN THEORY OF THE LAST MULTIPLIER. 



1. A STSTEM of n differential equations of the first order 
and degree containing n + 1 Tariables admits of n integrals 
of the form 

Wj, Wg, ...t«n being independent functions of the original vari- 
ables. When 71 — 1 of these integrals have been found they 
enable us to eliminate w — 1 variables, with their differentials, 
from the given system of equations, and so to obtain a single 
final differential equation of the first order between the two 
remaining variables. The final equation admits of being 
made integrable by a factor, and its solution so found would 
constitute the n^^ and last integral of the system. We pro- 
pose in this Chapter to develope the theory of the above 
integrating factor as established by Jacobi. The term ' prin- 
ciple of the last multiplier,' which is more usually employed, 
seems objectionable; for the essence of Jacobi's discovery 
consisted not in demonstrating the existence or the nature of 
the last integrating factor, but in the peculiar foi'm of the 
method which he gave for its determination, and in the rela- 
tions which are implied in that form. The discovery may be 
briefly said to consist in this; viz. that instead of forming by 
means of. the w — 1 known integrals the final differential 
equation between two variables and applying methods analo- 
gous to those of Chap, v., to determine its integrating factor, 
we construct antecedently to all integration a linear partial dif- 
ferential equation of the first order, any one integral of which 



i 
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■will enable us to assign an integrating factor of the final differ- 
ential equation, whatever the order of the previous integrations 
may have been. Again, this partiSi differential equation de- 
pending for its construction only upon the form of the system 
given, we can often by examining it affirm beforehand that if 
all the integrals but one of the system be in any way found, 
the final integral will be deducihle by quadratures. This 
happens in the case of the most important of all systems of 
differential equations — tliat of Dynamics. 

Further, an ordinary differential equation of the n'" order 
being reducible to a system of n differential equations of the 
first order, Jacobi'3_ theory may here also enable ua to pre- 
dicate the poBsibility of the last integration when the previous 
integrations have been effected. 

Beginning with a single differential equation of the first 
order reduced to the form 
I dx _dy 

~x~ y 

in which X and T are functions of the two variables x and », 
we know by Chap. v. that the integrating factor /i will be 
given by the solution of the partial differential equation 

^+^Q-0 (1), 

ax ay ^ ' 

the form of which should be carefully noticed. 

Consider next a system of two differential equations of the 
first order expressed in the genera! form 

dx di/ ds , , 

X-f°2 ^^'' 

X, Y, arid Z being functions of the three variables x, y, z, and 
suppose one integral, represented by 

*{^.y.e}=« (3), 

to be known. The function ^ [x, y, s), or, as we shall express 



^ 
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it for brevity, ^, will obviously satisfy the .partial differential 
equation 

^s+^^+^s=»-" w. 

of which indeed the given equations form the Lagrangean 
auxiliary system ; see Chap. xiv. 

If from the given integral we determine ^ as a function of 
X, y and c, and substitute its value in the first of the given 
differential equations, viz. 

dx ^dy 

the latter will be converted into a differential equation be- 
tween X and y. But we may leave to the equation ils prior 
form, provided that we regard X and T as Sanctions ot the 
variables x and y, both explicitly, as they appear therein, and 
implicitly as they are iavolTed in z. And this being so, the 
equation (1) will become 

d (fiX) d ifiX) dz d(fiY) djfiTj dz^ 
dx dz dx dy dz dy ' 

The values of -r- and ^- in this equation must be found 

dx dy 

from the known integral (3) ; they are 

dz _ d(l> d<f) dz ^ d<f> d<f> 
dx^ dx ' dz^ dy dy ' dz^ 

substituting which we have 

d{fiX) d^ d(jiX)d^ d(fiY) d(f> d{iiY)d^_^ .^. 
dx dz dz dx dy dz dz dy " 

This then is the partial differential equation for determin- 
ing /A. But the construction of this equation supposes ^ to 
be known. We propose to shew that fi can be determined by 
a process in which the only partial differential equation to be 
solved can be constructed without the knowledge of <j). 
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Since bj actual differentiation 

AfA^^-^ — fA^^:—*^^ — ^ 
dx\ dz) dz\ dx)'^ dx dz dz dx* 

it follows, writing fiX for A, that 

d(jiX) d<f> d(jAX) d^^d f ^d<f>\ ^ f -y^] 
dx dz dz dx" dx y*^ dz) dz y^ dxj ' 

Similarly 

d(^Y) di> dijiY) d4>^d r Y^]^ — ( Y^] 
df d» dz dy dy\ dz) dz\ dy)^ 

Lastly, we have 

Now adding the last three equations together we see that 
the first member of the result vanishes by (5) : we have thus 

-lKS)-s«)-l('^t)-»- 

The second line of the first member is equal to 
and therefore vanishes by (4). There remains then 
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Hence if we put 






we have 



d(MX) ^djMY) ^djMZ)^^ ,g, 

dx dy dz 



If then hy the solution of this equation a value of M diS" 
tinctfrom he found, the function -r? will he an integrating 

dz 
factdr of that final differential equation which remains when z 
has been eliminated from the system (2) hy means of any 
known integral <f>=:c. 

It will be observed that the equation for M is analogous in 
form to the equation for fi in the previous system. And this 
suggests the form of the general theorem. 

Thus proceeding to the case of a system of three equations 

dx _^dy ^dz _ dt 

we see that if 

be a known integral, -^ therefore satisfying the equation 

X^+Y^ + Z^± + T^^O (7), 

dx dy dz at 

then the system 

dx ^dy ^dz 
'X'^lt^'Z 

will virtually involve only the variables a?, y, «, since t 
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through the known integral becomes a function of a?, y, z. The 
equation (6) now becomes 

djMX) d{MX) dt d{MY) ^ d{MY) dt 
dx dt dx dy dt dy 

djMZ) d{MZ) dt^^ 
dz dt dz ' 

or putting 

dt ^ dyp" ^ dyp" 
dx dx ' dt 

d{MX) dy^ d{MX) d^ d{MY) df d{MY) dy^ 
dx dt dt dx dy dt dt dy 

djMY) d±_d{MY)d±^^ 
dz dt dt dz ^ 



and this is equivalent to 

dx dy dz dt 

and therefore becomes on rejecting the term in the second 
line by (7), and putting 

d(sx) d(irr) djiTZ) diNr) _^ 

dx dy dz dt '" ^ '* 
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If from this equation a value of N distinct from be ob- 

N 
tained, then -3f=^, and therefore 

'dt 

N 



/A = 



dt dz 



This is the final multiplier, i. e. the integrating factor of 
tti^&BAL differential equation between x and y which remains 
when z and" ^ baw^^been eliminated from the given system by 
means of the two knomr integjrals. In calculating fi from 
the above formula we must proceed ai^ follows. The value of 

-^ must be found from any given integral -^ = ^^ biot that of 

1? must be found from another integral from which by means 

of the former one t has been eliminated. Thus the general 
forms of the integrals will be 

<l> (a?, y , z, c) = c'. 

Lastly, the values of -^ , -^ found as above, and that of 

N given by any solution (distinct from 0) of the partial dif- 
ferential equation (8) having been substituted in the expres- 
sion for fly we must eliminate z and t from that expression by 
means of the two known integrals. The resulting function 
of oj, y, c and c' will be the integrating factor sought. 

'■ The reasoning above employed is in its nature quite inde- 
pendent of the number of the equations of the original sys- 
tem. The general theorem to which it leads may be tnus 
stated. 

Theorem. The system of n differential equations 
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■ 1)61115 given, if a system of n — 1 integrals 

be so reduced "by elimination that tlie varialjle y, shall not 
appear in tf>,, the variables y,, i/^ shall not appear in tf>^, and 
80 on, then the integrating factor /* of that hnal differential 
I equation between x and y, will be given by the formula 



in which M represents any integral distinct from of the par- 
tial differential equation 

(ic dy^ "" dy^ 

In applying this theorem the expression for ij. must lie 
freed from all the variables except x and y^ by means of the 
given integrals. 

This isJacobi's theorem. On account of its great importance 
I propose to give another demonstration of it founded upon 
the Calculus of Variations. 

2. Second demonstration founded upon (he Calculus of 
Variations. 

It will be most convenient to present the proposed system' 

of differential equations under the symmetrical form 

dx, _dx^_ _ diP^ 

the independent variables being x , a:,, x^ of which 

X,, Xf, .... X, are any functions. We have thus n — 1 difft- 
rential equations, and we are to seek the integrating factor of 
the differential equation which remains when by means of 
n — 2 known integrals n — 2 of the variables with their diffe- 
rentials have been eliminated. 



i 
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Suppose P= c to be any integral of the system, then P 
satisfies, and it suffices that it satisfies, the partial differential 
equation 



-^1 "izr + -^a TJT • • •• + -Xj, -j^ — 



dx. 



dx^ 



dx. 



(a). 



Now if in place of 0?^, ajg* ••••^« we introduce a new sys- 
tem of independent variables m^, i^a, .... i^^ which are functions 
of the former, then we shall have 

^dx^ ^dx^"" ^ dx^ 

_TjdP jjdP jj dP^ 

^du^ ^ dui"'* ** du^* 

U , Ui, ..*. Un being functions of u^, u^, .... u^. And by the 
theory of the transformation of multiple integrals, 



dP 



fVx^ + X^^ +X^dxdx dx 
J [^'dx,^^'dx,""^^''dxj'^''^''"''^'' 



where 



H= 



du^ 



du^ 
dx.. 



du^ 
dx^ ' 



d(K 
dx^ 



The foregoing equation we may express in the form 
SI -Xi-T- rfa?, cfecj, ... dx^=S I -U^du^du2...,du^. 



Hence, representing by B an operation of differentiation 
which affects only the form of P as a function of o;^, a;^, ... fl*. 
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or of u^j 11,, ....ti^, and not the independent variables them- 
selves, we have 

and therefore integrating by parts and equating the portions 
on each side which remain under the sign of n-fold inte- 
gration, 

21 -j-^BPdx^dxt....dx^ 

Whence again transfqrming the integral in the first member 

^ p dXi -jp <?Uj du^ .... du^ 
J dXi H 

= 2j -^(■^]SPdu^du^....du^, 

and this being true quite irrespectively of the form of P, 
we have 



H^ dx,~^ du\H)' 



In this' equation Jacobi's theorem is virtually contained. 
For let the given equation be multiplied by any factor. Then 
changing in the above Xi into -MZi, and Ui into MU^^ we 
have 

1 ^ d{MX:j ^^ d (MU\ 
H^ dx^ ^ duXH J* 

Hence, if if be determined to satisfy the equation 

aXi 
B.D.E. XL 14 
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we shall have 

si(¥)=« (»)• 

This is wholly independent of the relations connecting 
u^, u^,...u^ with cTj , a?j, , . . . . a;„. Now choose the w — 2 variables 
Wj, w,, .... w,^ so that ««j = Cj, Wj = Cj, .... w,^ = c„_j shall be 
integrals of the given partial diflferential equation (a). Then 
that equation transformed becomes 

of which the auxiliary ordinary equation is 

At the same time the equation (J) becomes 
d /M^ \ d fM 



du^i 



i^^yuP-)-"- 



M . 
Hence -n- is the integrating factor of the preceding diffe- 
rential equation between u^^ and w„. 

Jacobi's theorem in its most general form is thus seen to 
be the following 

Theorem. If the system of diflferential equations 

dx^ ^ dx^ dx^ 



XT — "vr •••• — V" 

be transformed by the introduction of a new system of vari- 
ables Wi, Wg, ....w„, so chosen that 

shall be integrals of the given system, then the final differen- 
tial equation between u^^ and u^ shall have for its integrating 



\ 
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M . 



factor -^, in which M is any function satisfying the partial 
differential equation 

djMX,) djMX,) d ( JfX,) _ 



dx 



dx, 



dx. 



and H stands for the determinant 

du, du. 



dx^ 

du^ 
dx/ 



dx^ 



m 



n 



dx^ 



The form of Jacobi's theorem obtained by the previous 
demonstration may be deduced from the above by choosing 
for u,^,Un two of the original variables, for example a;,^, a?„, 
and transforming the integrals w^, u^, .... Wn_2 so that u^ shall 
contain only x^ ... a?„, u^ shall contain only x^*..x^j and so on. 



Examples. 

3. Jacobi has established by means of the above theorem 
the very remarkable theorem that in any ordinary dynamical 
problem the forces depending not upon the time but upon the 
material constitution of the system, if all the integrals but 
two of the dynamical equations are found, the two remaining 
integrals can be found by quadratures. 

1st. In a dynamical system of free points the forces act- 
ing upon which depend only upon the position of the points, 
we have if we represent the entire system of rectangular co- 
ordinates taken in any order by a;, y, «, ... and the correspond- 
ing resolved forces divided each hj the corresponding mass 
by Xf Yf Zf... the system of equations 



d^x _^ y. d^y _ 



Y 



U— a 
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or putting 

dx , dy _ f 



• • • 



;ja __dx _dy _ dx* __ dy 
X y X Y 

Now as X, F... do not contain t we may consider first the 
system 

dx _dy _ dx' _ e?^ 



and it is evident that if we can find all the integrals of this 
system, t will be given by the equation 

» 

dx 



{dx 



x' having been first converted by means of the supposed in- 
tegrals into a function of x. 

To determine the last multiplier of the system last written 
we have first the equation 

d jMx') d {My') d{MX) d{MY) ^^ 

dx dy "" doi/ dy '" ' 

which since X, F... do not contain a?', ^ ... is satisfied by 
ilf = a constant Giving to the constant the particular value 1, 
we see that if 



^i^^if **j|"~ ^a) •••• ^n-«~" ^1 



fr-% 



are 71 — 2 integrals of the system, and if by means of these 
we eliminate n — 2 of the variables and construct the differen- 
tial equation between the two remaining variables, the inte- 
grating factor of that equation will be -^ , in which fi" is ike 
functional determinant of Wj, Wg, .... w„. . ^^ 



/ 
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2niily. Suppose tlie system subject to a material comies- 
ion wbicli establialiea an equation of condition among some 
or all of tlie co-ordinates. If we represent the co-ordinates 
taken in any order and multiplied each by the square root of 
the corresponding mass by a;, ^, ... the corresponding resolved 
forces by Xj F, ... and the equation of condition expressed by 
means of the above modified co-ordinates by = 0, the diffe- 
rential equations will be 

d^ ' dx' d^ dy ''"' 

the transformation above employed reducing all the equations 
to the same type. [See tbe next Chapter.] 

Making 

dx , dy , 



^ 



the system becomes 



dt = 



^ dy _ dx dy 

"""""^x+xf-'r+xf 



and the Jacobian equation for M becomes 



dx ^ dy ■ 






Now ^ does not contain a:', y' .... Let us inquire whether 
it is possible to determine ,1/" also as a function oix,y.... 
without x, y\... so as to satisfy the above differential equa- 
tion. 



J 
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The equation would become 

dao ^ dy \dx dx dy' dy j ' 

or if we write 

^Tx^y dy ^' 

*^-^(eS-||^ •■•■)-». 

and from this we must eliminate X. 

Now since ^ = 0, we have by differentiating and puttir 
dx ^ , 

and again differentiating 

Md^^d^ dy' _^ 

dx dt dy dt "" ' 



or since 



a,'«^+«'»^4. A.9x'v' ^^ + 
"^ da?^^ dy'^'-'^''"'^ dxdy^'- 
^^dx^^ dy^"" 



AKT. 4.] OF THE LAST MULTIPLIER. 215 

and diflferentiating with respect to a?', 



or if we make 



Similarly 



ahith--^' 



<^>"'- 



<*$0-o- 



Therefore 



dx ax dy ay 

(d\ d(f> dK d<f> \ n a 

or 

« 1. • • a\ dd) a\t ad> 

and now elimmatmg -j-, ZT '^'1~' T^"^ '" 

we obtain QBM- MBQ^O, 

which is satisfied by Jf = Q. 

4 [Among Professor Boole's manuscripts I found five 
pages in German, forming part of a memoir, which was pro- 
Dably intended for Crelle s Mathematical Journal, The 
memoir was to have discussed two applications of the Calcu- 
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lus of Variations; one to the Jacobian Theoiy of the Last 
Multiplier, and the other to the Solution of Pfaff 's equation 

X^dxi + X^dx^+ +X^dx^=^0. 

But there is only a single paragraph relating to the second 
application. 

The manuscript contains the same demonstration of the 
Jacobian Theory of the Last Multiplier as in Art. 2 of the 
present Chapter; after this demonstration some remarks occur 
of which the substance will now be given.] 

It is worthy of notice, that Jacobi in the 36th volume of 
Crelle's Journal, deduced by the aid of the Calculus of Varia- 
tions the result on which the preceding demonstration of the 
Theory of the Last Multiplier depends. In fact, he shewed 
that if V denotes any function of 



dz 



dz 



a?,, a?,, ... ic„, z, ^^, •.. ^^^, 



and Fbe transformed by the introduction of a new system of 
independent variables w^, u^, ... w«, then the following rela- 
tion holds, 



'dV d dV 



dz dx^ J dz 
dx^ 

d[b.V) d d{LV) 

dz 



dz du 



* d 



du. 



d_ dV\ 
dx^ -J dz j 

d (?(AF) 
du^ , dz 
du^ 



where 



A = 



dx^ 
du^ ' 



dx^ 
du^ ' 



dx^ 
du^ 



dx^ 

du. 
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Jacobi applies this result to the transformation of the ex- 
pression 

d'V cPV cPV 

But neither Jacobi himself, nor any other person, so far as I 
know, has drawn attention to the application of the result 
which I have given here. 

[The substance of the single para^aph relating to the 
second application of the Calculus of Variations will now be 
given.] 

Clebsch has earned the thanks of all who are interested in 
the higher parts of the Theory of Differential Equations, since 
he has performed the same service for Pfaff s problem as 
Jacobi did for the Theory of Partial Differential Equations of 
the first order, and thereby for the equations of Dynamics. 
But while I recognise the great importance of the results, I 
consider it desirable to give a simpler deduction of the system 
of partial differential equations therein involved, and on which 
the other results depend. 
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THE DIFFERENTIAL EQUATI0N:S OF DYNAMICS. 



[It will be seen that this is only a fragment of the Chapter 
which was to have appeared under this title.] 

I do not propose in this Chapter to discuss the origin and 
interpretation of the diiferential equations of motion or to enter 
into those details of their application which are found in all or- 
dinary treatises on Dynamics. But they constitute a system 
analytically so remarkable from the forms in which it is 
capable of being expressed, and from the general methods of 
integration which emerge out of those forms, that they are 
well deserving of a special attention. 

Referred to rectangular co-ordinates the differential equa- 
tions for the motion of a system of points free or connected 
are 

"* 'd^-^^^d^''^''dx'-- 



THE DIFFERENTIAL EQUATIONS OP DYNAMICS. 

Here m ia the mass at the point [x,y,z), m' that at {x',y',z'), 
X, Y, Z the resolved forces at (a;, y, £) tending severally to 
increase those co-ordinates, and so on. Lastly 



are the equations of condition each of which may involve all I 
the co-ordinates, and X, ^... are indeterminate multipliers. 

The above is usually termed the first Lagrangean form of the 
differential equations. In applying it we must either elimi- 
nate \, jji,... from the given equations, and then by the equa- 
tions of condition just so many of the co-ordinates with their 
differentials, or ne must retain X, ^,. . . as variables so conditioned 

"fliat the values of -^ , -^ , ... in the system shall satisfy iden- 
ically the differential equations involving ~- , -^ de- 
rived &om = 0, 1^ = 0,... viz. the equations 



to 

4 



-.^ = 0, ^=0, 






The first Lagrangean system may by a slight ti-ansforraa- 
jlion be reduced to a form in which all the equations arc of 
.one type, viz, of the type which they would have if all the 
masses were equal to unity. 

For taking the first equation of the system and dividing^ 
1^ m* we may express the result in the tbrm 



from which we see that if a;, y ... had been taken to represent 
Ihe entire system of co-ordinates taken in any order and mul- 
"rolied each by the square root of the corresponding mass, and 
C, Y ... the corresponding resolved forces taken in the same J 
order and divided each by the square root of the correspond- 1 
Dg mass, the system of equations would have been 
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all being of one type. In general investigations this form 
is to be preferred. 

From the first Lagrangean form another known as the 
second Lagrangean, and from this again a third known as 
the Hamiltonian are derived. The second Lagrangean form 
is properly speaking an expression for the effect of a trans- 
formation of co-ordinates in the most general sense upon the 
original system, i.e. of a transformation which in place of 
a?, y, ... the entire system of given co-ordinates substitutes 
a new system of variables f, ^, ... the expressions of which as 
functions of a;, y, ... are known. It is not necessary that this 
new system of variables should be co-ordinates in the proper 
sense of that term, determining three by three the positions 
of the several masses; it suffices that they should in their en- 
tirety determine and be determined by the co-ordinates given. 

The second Lagrangean form may be established as 
follows : 

Differentiating the equations ^ = 0, '^ = 0,... with respect to 
any one of the new variables f we have , 

r 

d(f> dx d<f) dy 

dx d^ dy d^ '" ~ ' 

d-yjr dx d^ dy _ 
dx d^ dy d^"* ' 

whence if we multiply the equations of the given system by 

dx dy 111 1 

-757, -ji, ,... and add, we have 
di d^' 

dac^ d^x dy^ d^y _^ Y— V^ 
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CHAPTER XXXIII. 



OK THE PEOJECTION OP A SUHPACE ON A PLANE. 



[The following memoir -was found among Professor Boole's 
manuscripts; a Title and IntrodEctory Eemarks were to have 
been prehxed, but with thia exception the memoir appears to 
be finished for publication. It 13 sufficiently connected with 
the subject of Differential Kr[uation3 to find a place in the J 
present volume. 

The memoir by Sir John Hersche! to wliich allusion is ] 
made is entitled, On a new Projection of the Sphere; this was 
read before the Eoyal Geographical Hociety of London on 
the 11th of April, 1859, and was printed as part of the Jonmal 
of the Society, Vol. xxx. 1860, pages 100. ..106. A chart o" 
the World on Sir John flerachel's projection has been pub* 
lished by A. and C. Black of Edinburgh. 

The history of the subject will be found in Chapter xxin 
of tlie Coup d'a^il htstorique sur la Projection des Cartes i 
Oiographie... Par M. D Avezac, Paris, 18G3. 

For the materials of thia introductoiy notice I am indebted 
to Sir John Herschel.] 

1. Let 3!, y, 3 be the rectangular co-ordinates of any point ' 
on the given surface; ^',y' the co-ordinates of the correspond- 
ing point on the plane of projection. Let the equation of the 
given surface be 



F{x,y,s)=0; 



or, for simplicity, 
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The condition of projection upon which Sir John Herschel's 
investigations are founded, and which we shall adopt here, is 
that of the similarity of corresponding infinitesimal areas on 
the surface and on the plane. The object of the problem then 
in general is the discovery of the mode in which a;', y depend 
upon a;, y, and z in accordance with the above condition; its 
object in any particular case is the determination of a?', y' as 
functions of a;, y, 2, 

Eegarding then a?', y' as ultimately functions of a?, y, z we 
have 

dx =-r- <ix + -r- dy-\-'-r dz. 
ax ay ^ dz 

in which da?, 6?y, dz are not independent, but are connected 
by the condition 

dF 7 dF , dF , 

J— oa: + -- ay 4—7- dz = 0. 

dx dy ^ dz 

Now for brevity write 

da;' __ dx* _ - dx' __ 
dx" ^ dy" ^ dz" '* 

M.^a! ^=&' ^^e 

dx ^ dy ^ dz ' 

dF . dF ^ dF ^ 
d^-^> ^=^> ^=^' 

then 

dx' = adx + hdy + cdz (1), 

dy =^adx + b'dy + c'dz (2), 

0=^Adx + Bdy-^ Cdz (3). 
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Now the condition of the similarity of infinitesimal cor- 
responding areas may be resolved into the two following 
conditions, viz.: 

1st. The equality of their corresponding angles. 

2ndly. The proportionality of their corresponding sides. 

And these conditions we shall introduce separately. 

1st. Assuming any point x\ y* on the plane of projection, 
let x' alone vary, and the infinitesimal line generated is da?', 
while (since dy' = 0) (2) and (3) become 

a!dx 4- Vdy + c'dz = 0, 

Adx + Bdy+ Cdz = 0, 

whence, if we write 

L-=^Bd'-Ch\ M=-Ca-Ac', N=Ah'-Ba\ 

. dx dy dz ,. 

we have T^M^N (^^' 

so that the direction cosines of the infinitesimal line on the 
surface F corresponding to the line dx' on the plane {x\ y') 
will be 

L M N 

(2,2 + jkP+JV^)4' [U + M'+N^)^' {U-^M^+N^)^ ^ ^' 

In like manner, if y' alone vary, we shall find for the 
direction cosines of the infinitesimal line on the surface F 
which corresponds to dy' on the plane 

L' ^ N' ... 

{L"+M''+N'')^' (Z'^+if'+i^*)*' {L'''\-M''+N")^ ^ ^' 

where L'^^Bc- Cb, W ^ Ca-Ac, IfT ^Ab -Ba. 

By the first of the conditions of similarity the ungle be- 
tween these lines on the surface must be a right angle since 
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dx and dy' are at right angles. Hence we have, from 

(5) and (6), 

LL' + MM' + NN'^O (7). 

2ndly. The ratio of the length of the element dx' to the 
corresponding element on the surface is 

dx' 
HJdx' + dif + dz^ ' 

or, ty (1), 

adx + hdy + cdz 
'4^+df-\-dz'' 

and therefore ty (4) 

aL + bM+cN 

equating which to the corresponding expression for the ratio 
of the length of dy to that of its projection on the suxface, 
we have 

Now if we substitute for i, My N, L\ M\ N' their values, 
we shall find 

aL + hM+cN=-A{Vc-hc')+B{c'a-ca') + G{a'b - aJ'), 

a'r + VM'+ c'N' = A [hc'^ h'c) + B {ca!- c'a) + G {aV-^a'b), 

and the second members of these equations differ only in sign. 

Thus (8) may be expressed in the form 

U {Vc - be') + B (c'a -ca')-\'G (aJ- a J)l 

X I T -rh = 0...(9). 
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But the first factor of the first member of this equation 
being the determinant of the system 

adx + hdy + cdz = 0, 

ddx + Vdy + cdz = 0, 

Adx-\rBdy + Clfe = 0, 

expresses when equated to ^ero the condition that if in the 
system (1), (2), (3) dy* vanishes da?' shall also vanish; and 
ax and dy being independent, this condition cannot be satis- 
fied, so that (9) reduces to 

1 L__=o 

whence 

i'* + -¥^ + i^-i*-Jf-2^«0 (10), 

and this, with (7), will fully express the conditions of simi- 
larity* 

2. If we multiply (7) by 2 V— 1, and add and subtract 
the result firom (10), we obtain the equivalent system 



-1)« = 0] 
-1" = 0^ 



(r + z V- 1)» + {M + M^j^xY + (i^ + iv^ v-i)« = 

(II - L V^)* + (ilf - if V=n)» + (J^- N^\) 
XT Tt . -r r"^ dF dsd dF dx* 

Now -^ ± -^ v-1 = j;: -3r - -jt :s; 

ay €U8 az ay 



\dy dz dz dyj 



dF djx'tyW'-l) dF djx'ty'*/^) 
dy dz dz dy * 

an. E. II. V& 
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Writing then 

a;' + y ' V^ = M, a:' - y V ^ = t^, 



we have 



f--^ dFdu dF du 
i'-.iV^ = — — - — — 



In the same way 



JJf'-f JfV-l^^^-T- :3~ 



dFdu^dFdu 
dz dx dx dz^ 



' *" dz dx dx dz^ 

'^ " dx dy dy dx^ 

dx dy dy dx" 



Substituting which in the system (11) there result 



fdFdi^ 
\dy dz 



dz dy> 



+ 



fdF du 



fdFdv 
\dy dz 



dFduV /dFdu 
dw \dz dx 

/dF 
\dx dy 

dF dv\* fdF dv 
dz dy) \dz dx 

/dF^ 
\dx dy 




dJFdi^ 
dy dx) 



= 




^...•(12), 



dy dx) 



2.] 
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Hto whicli we may give the Bomewliat more convenient form 

_/dFda dFdu 
\dx dx dy dy 



©'- 



I")} 



dz dz) 



..(I). 



m^m-mm^m-m 



(dF do^ dF dv dF dv\' 
\dx dx dy dy dz dzj 



■■(11). 



These are partial dL£Ferential equationB of the first order, 
■serving to determine u and w as ftmctiona of x, y, z. 



I 

V> But it is not necessary to solve tlie equations in their 
W general form. For, x, y, and z being connected by the equa- 
tion of the surface, the above equations may always be so 
reduoed as to involve only two ludependent variablea. As 
latitude and longitude determine the position of a point on 
the earth, so two co-ordinates of any given species wdl deter- J 
mine the position of a point on the given surface, and these 1 
co-ordinates, when fixed upon, become the independent varia- 1 

tbies of the problem. 
Let s and ( represent such co-ordinates, and let their ex- I 
pressions in teims of x, y, z give 
ffhi 
reci 
1st 



{x,y,z), t = 4>^{x,ij,z), 
'hich equations combined with that of the given surface will I 
reciprocally determine x, y, s as functions of s and t. Then " 
Ist the differential coefficients of F which in the equations 
(I), (II), arc functions of a;, y,s maybe transformed into func- 



tidly, we have 

du _ du ds du dt 

dx da dx dt dx' 

du du ds du dt 

dy ds dy dt dy ' 

du du da du dt 

de da dz dt dz' 
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and as T- > ^... are known functions of a?, y, «, they also 

are expressible in terms of 8 and t. The result of these sub- 
stitutions will then be to convert (I) into, a partial differential 
equation in which u is the dependent and s and t the inde- 
pendent variables, and this equation being, like (I), of the 
first order and second degree in the differential coe^cients of 
w, will be of the form 

For V we shall have an exactly similar equation with the 
same coefficients. 

The above equation is, by the solution of a quadratic, 
resolvable into two equations of the form 

du du du du 

To these correspond the respective auxiliary equations 

dt+\d8=0y dt + \ds = (13). 

If the integrals of these are 
respectively, then we have 

Now V being determinable by an equation of the same 
form as u, it follows that of the above two values of u one 
must be assigned to v, so that the solution of the problem will 
be contained in the system 

u = <f>{S), v = yjr{T), 
or in the system 
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The particular forma of the arbikary functions i^ and -^ 
will depend solely upon tbe nature of the problem under con- 
sideration. 

One other point remains to be noticed. ■ Tlie first mem- 
bers of (12) are essentially positiTCj being composed of 
squares ; so are then the first members of (I), (II) ; and so, 
if the intermediate transformations are real, is the first 
member of the equation wlioae coefficients are" P, Q, It. 
Hence the quadratic determining \, \ will have imaginary 
roots of the form a + /3^/-l. Ultimately therefore it will suf- 
fice to integrate one equation of the system (13) and then to 
deduce the solution of the other by changing J— 1 into 
- V- 1. 

3, Application of the above formulce when the given sur- 
face ia an Mate spheroid, such as the earth. 

Let the plane of the equator be that of projection, the 
centre being the origin. Let the co-ordinates x, y pass 
througli the meridians of and of 90° respectively, and z 
through the poles. The equation of the surface will be 

"^^4-^ (1*). 

where a is the earth's equatorial, b its polar radius. Let also 
the latitude of the point x, y, z be represented by s, the 
longitude by t. We have 







dF_2x ^F_2^ dF_ 

dx a* ' dy a" dz~ 


2z 


and 


[ substituting in (I), 






(^ 


^-i^{(S^®'^(£)} 








rx du ^y *_, 


z d„ 

ff a. 



^1 
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<.,if.erep««.|,1>y*^ 






(15). 



Now as 0?, y are rectangular co-ordinates in the plane of 
the equator^ and x passes through the first meridian^ we have 

^ = tan «. 

X 

Again, representing in the annexed figure the meridian of 
the point P, or (a?, y, z) 
touched by the straight line 
QR in the same plane, we 

have (7Jf=V ^+y, MP=z. 

Therefore if V^*+p=r, the 
equation of the meridian is 

— 4- -1 = 1 




that of the tangent 






r', z being current rectangular co-ordinates of the tangent 
Hence 



a^z 



t3inCQR = Tr- = 



h'z 



But CQR = latitude. Therefore finally 



8 = tan 



-1 






(16), 



and we must now transform (15) so as to make 8 and t the 
independent variables. 
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From the above equations combined with (14) we find 

__ ah cost ^ ahsint ^ a tan s . . 

^""VA» + tan»1J' ^""VAV+tan^l' ^" AVA» + tai?^'"' ^ ^' 



and substituting in (15), 



— (co8«-^ + sm^-T- + tan«^J =0 (18). 



Again, 



du _^du ds du dt 
dx" ds dx dt dx^ 

du ^du ds du dt 
dy" d^ dy dt dy^ 

du^^ ds du dt 
dz" ds dz dt dz* 

Now ~ — — hhx — sing cosg costjH 

where jEr= A* + tan*s. In like manner 

ds _ — sing cosg BJatjH 



dy 


ah 


da _ 

c5 


h COS** JH 
a 


dt 


— tantJS 
ah ' 


dt 


costi/H 

ah ' 


dt 
de~ 


0. 
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Hence 

du JS I . du . ^du\ 

-r- = ^^^- — sms COS* cos ^ -7- ^ sm ^ -J- , 
dx ah \ ds dtj 

du SS [ . . ^c?M . \du\ 

-5- = ^^S- -- sin « cos 5 sm ^ -7- + cos ^ -y- , 
dy ah \ ds dtJ 

dz ah \ dsj 

Substituting these values in (18), and dividing by the com- 
mon factor -gTg- we have on teduction 

ah , 

(|)V cos'. {1 + (A' - 1) cos'«r ( jy = 0, 
which is resolvable into 

^ - /Ti C0S5 {1 + (A'- 1) cos»«} ^ =0, 

^ + 7^ cos 8{1 + {h* - 1) cos*5}^ = 0, 

partial diflferential equations of which the integrals are in- 
cluded in the common formula 



^°^ jco3s{H-(A»-l)co8»«} 

J COS* ^ W 1 4 (/i* - 1) cos'^ 
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f da ^fco&sds / . , a* - 5'\ 
= eM ; a ' 2 > since e' = — a---l 



logtan(^ + |)+|logJ 



1 — e sin « 



V4 2/ 2 ^l + esm^ 



='»«{(r57i^)'-'a-l)}- 



Hence 



•'-*hl(^^/'"(i^'^)i*'^^]' 



oip changing ^ (^) into ^ (e*), 



«=*{(li7S-:)''-(!+i)-i. 

"*{(t^)''«»(M)-'-}- 

Let r and 5 be the polar co-ordinates of that point in the 

Elane of projection which corresponds to the point whose 
ititude and longitude on the surface are « and t; and let 

^ /I -6 sin5\l . fir , s\ 

then the complete solution assumes the very simple form 

r€^Vir= ^ (fi^«V=i), r€-«^^= t {Se^^^ (III). 

Of particular deductions the most interesting is that which 
arises from the supposition that the parallels of latitude are 
projected into circles round the pole. This requires that r 

B.D.E. II. 16 
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should be mdependent of t, a condition which is satisfied in 
the most general manner by assuming 

we then find 

whence, on multiplication ^nd division, 

whence, A and B being new arbitrary constants derived from 
(7 and a' 

If we observe that 5 and t should vanish together, we have 
-8 = 0, and the equation 0= ±nt shews that the surface of 
the sphere will be projected into a sector of a circle, the arc 
of which is to the circumference of the circle as n : 1. Thus, 

if 71 = -, the sphere is projected upon a quadrant^ and so 

on. 

The other equation gives 

I \4: 2/j Vl + ^sms/ 

If 5 = we find r = -4, whence A is the distance of the 
ec{uator from the pole in the plane of projection, and if that 
distance, which is arbitrary, be assumed as the unit, we liave 



\\4: 2j) \1 +e Bin sj 



for the distance from the pole of that parallel whose latitude 
is ^. We may give to this expression a better form by 
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assuming jp^-^ +8, and introducing an auxiliary quantity j 
determined by the equation 

eco8p = co8j. 
We have then 

-H)"("'i)" 

The following table gives the values of r for the sphere and 
for the spheroid whose eccentricity is '08 (which is about that 
of the earth), for each ten degrees of polar distance, for the 

values w = 1, and w = j • 



PoUr 


n= 


»1 


n 


1 
■"4 


DiBtanoe. 










Sphere. 


Sphtroid. 


Sphere. 


Spheroid. 


10« 


•0876 


•0880 


-5439 


•6447 


20® 


•1763 


•1774 


•6480 


•6490 


30® 


•2679 


•2694 


•7196 


•7206 


40<> 


•3640 


•3668 


•7767 


•7777 


60<> 


•4663 


•4682 


•8264 


•8272 


60® 


•6774 


•6792 


•8717 


•8724 


70* 


•7002 


•7017 


•9148 


•9163 


80» 


•8391 


•8400 


•9571 


•9674 


900 


1-0000 


1-0000 


1-0000 


1-0000 


100» 


11918 


1-1904 


1-0448 


1-0445 


IIQO 


1-4281 


1-4260 


1-0932 


1-0926 


120'> 


1-7321 


1-7266 


1-1472 


1-1463 


1300 


21446 


2-1367 


1-2101 


1-2089 


1400 


2-7476 


2-7340 


1-2876 


1-2869 


1600 


3-7321 


3-7114 


1-3899 


1-3880 


1600 


6-6713 


6-6372 


1-6432 


1-6409 


1700 


11-4301 


11-3681 


1-8387 


1-8368 



THE END. 
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